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ABSTRACT 


Measurements  of  turbulence  field  data  are  reported  for  recirculat¬ 
ing  utter  flows  under  steady-state  and  transient  flow  conditions  In  an 
Improved  1/15  scale  FFTF  outlet  plenum  test-cell.  A  dual-channel 
Laser  Doppler  Anemometer  Is  used,  and  a  unique  method  for  performing 
transient  measurements  Is  developed  and  analyzed.  The  experimental 
behavior  Is  predicted  with  a  modified  version  of  the  VARR-II  fluid 
dynamics  code,  and  the  suitability  of  using  the  steady-state  turbulence 
model  closure  assumptions  for  analysis  of  transient  flows  Is  evaluated. 

A  relationship  between  K-q  and  K-e, turbulence  models  is  derived,  and  the 
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Steady-state  and  transient  mean  velocity  data  are  well -predicted 
throughout  most  of  the  test-cell.  Calculated  spatial  distributions 
of  the  scalar  turbulence  quantities  are  qualitatively  similar  for 
both  models;  however,  the  predicted  distributions  do  not  match  the  datp 
over  major  portions  of  the  flow  area.  The  K-ovmodel  provides  better 
estimates  of  the  turbulence  quantity  magnitudes  •  Z  '  r  tr>'  ''■v. 
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The  predicted  results  are  found  to  be  highly  sensitive  to  small 
changes  In  the  turbulence  model  free  parameters  and  to  depend  heavily 
on  the  levels  of  Inlet  turbulence.-^  A  parametric  analysis  shows  that 
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cannot  be  significantly  reduced  by^slmple  changes  to  the  model's 
constants. 


The  transient  analysis  Indicates  that,  for  the  coast-down  rate, 
flow-speeds,  and  geometry  used  In  this  experiment,  the  turbulence 
••odel  closure  assumptions  are  adequate.  An  approximate  method  Is 
described  for  determining  the  applicability  of  this  conclusion  to 
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NOMENCLATURE 

constants  appearing  in  transport  equations  of  K-o  model 
ax,dz  donor  cell  coeficients 

8q  relaxation  factor 

CrC2’^y  constants  appearing  in  transport  equations  of  K-e  model 

AL  optical  path  length  difference 

e  turbulent  energy  dissipation  rate 

f  RC  circuit  cut-off  frequency 

f D  Doppler  frequency  (Hz) 

K  turbulent  kinetic  energy 

l  length  scale  or  Prandtl's  mixing  length 

L  laser  coherence  length  or  characteristic  length 

\  wavelength  of  laser  light 

ni  index  of  refraction  of  substance  i 

N2  number  of  fringes  in  the  measuring  volume 

a)  frequency  in  radians 

Up  Doppler  frequency  (radians) 

(jjq  frequency  of  laser 

frequency  of  Doppler  shifted  light 

w  frequency  of  reference  beam  light 

2 

P  pressure 

Re  Reynolds  number,  defined  with  respect  to  inlet  duct 

p  density  of  the  fluid 

a  turbulent  kinematic  viscosity 

av ,  a  effective  turbulent  Prandtl  number  for  transport  of  K  and  e 

t  time 

T.  t  integration  time 


trans. 


turb. 


T 

e 

u 

Uc 

U„ 


ui  uj 

jmr 


u'.v',  w‘ 


V 

W 


xrxj 


transient  time  scale 

turbulence  time  scale 

integration  time  constant  for  RC-circuit 

half  angle  of  the  beam  intersection 

mean  velocity  (vertical  or  X  direction) 

characteristic  velocity 

ensemble  averaged  velocity 

velocity  component  in  i-direction 

time  averaged  velocity 

kinematic  Reynolds  stresses 

kinematic  shear  stress 

fluctuating  velocity  components 

mean  velocity  (horizontal  or  Y  direction) 

mean  velocity  in  Z  direction 

general  Cartesian  coordinate 


CHAPTER  ONE 


INTRODUCTION 


1.0  Overview 

Accurate  design  of  many  types  of  fluid-flow  and  heat- transfer 
equipment  is  aided  by  accurate  knowledge  of  the  fluid  velocities, 
temperatures  and  turbulence  intensities  that  develop  within  the  particular 
mechanical  component  under  study.  Knoweldge  of  these  parameters  permits 
detailed  estimates  of  such  important  quantities  as  local  heat-transfer 
coefficients,  pressure-drops,  fluid/solid  mechanical  interactions  forces, 
and  thermal  shock  effects.  One  important  area  where  this  information 
is  needed  is  in  the  design  of  the  outlet  plenum  region  of  a  liquid-metal- 
cooled  fast  breeder  reactor  (LMFBR). 

The  efficient  heat  transfer  characteristics  of  liquid  metal 
coolant,  combined  with  the  relatively  high  bulk  temperature  change  across 
the  core  (approx.  300°F) ,  present  thermal  shock  problems  for  nozzles  and 
other  above-core  mechanical  components  in  the  event  of  a  reactor  scram 
from  full  power,  with  or  without  pump  coast-down  [II.  In  an  effort  to 
mitigate  these  effects,  current  outlet  plena  designs  provide  for  a  large 
inventory  of  "hot"  coolant  above  the  core,  and  a  "chimney-like"  entrance 
to  the  outlet  plenum.  The  chimney  forces  the  relatively  cool  core-outlet 
flow  to  mix  with  the  existing  high  temperature  fluid  in  the  plenum  so 
that  the  average  temperature  of  the  fluid  contacting  the  high-temperature 
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structure  is  raised  to  some  safe  value. 

One  method  of  predicting  the  success  of  this  design  is  to  solve 
numerically  the  applicable  fluid  dynamics  equations.  A  number  of  computer 
programs  exist  for  this  purpose,  and  newer  versions  and  updates  of  these 
codes  are  currently  evolving.  For  the  turbulent  flow  case  (the  most  likely 
and  hence  most  practical),  all  such  codes  employ  a  turbulence  model  to 
compensate  for  the  fact  that  averaged  versions  of  the  Navier-Stokes 
equations  present  more  unknowns  that  independent  equations.  Turbulence 
models,  then,  are  just  additional  equations  (they  may  be  algebraic  or 
differential)  that  close  the  set  of  equations  which  describe  a  particular 
problem,  thereby  allowing  a  solution  to  be  obtained. 

The  turbulence  models  currently  being  used  vary  widely  in  their 
complexity.  In  general,  they  all  depend  to  some  degree  on  heuristic 
reasoning  for  their  derivation,  and  contain  adjustable  constants  (model 
free  parameters)  that  are  set  based  on  existing  experimental  data. 

Although  the  majority  of  current  models  employ  free  parameters  that  were 
obtained  from  steady-state  flow  experiments,  they  are,  of  necessity, 
being  used  to  solve  transient  problems  like  the  LMFBR  outlet  plenum  flow 
described  above. 

The  purpose  of  this  work  is  to  examine  the  applicability  of  an 
essentially  steady-state  turbulence  mode)  to  the  solution  of  transient 
flow  problems.  The  method  used  is  to  perform  flow-mapping  measurements 
during  a  transient  in  a  simple  geometry  similar  to  the  LMFBR  outlet 
plenum  configuration.  A  comparison  is  then  made  between  the  measurements 
and  the  code  calculations  to  determine  the  accuracy  of  the  model's 


predictions.  In  this  way  it  Is  possible  to  determine  situations  In  which 
the  models  are  fundamentally  Inadequate,  to  identify  vital  input  data 
requirements  for  accurate  calculations,  and  to  guide  improvements  in  the 
models  for  application  to  transient  flows. 

1.1  Background 

This  effort  is  part  of  a  larger  program,  the  overall  goal  of  which 
is  the  validation  and  improvement  of  calculational  methods  for  analysis 
of  LMFBR  outlet  plenum  flow  transients.  The  most  general  situation 
involves  a  transient  in  both  flow-rate  and  fluid  temperature  occurring  in 
three-dimensional  space.  The  approach  taken  for  assessment  of  the  various 
calculational  models  and  methods  has  been  to  separate  the  relatively 
complex  general  transient  into  its  more  readily  studied  parts.  Results 
are  compared  to  the  numerical  predictions  for  that  particular  flow 
situation.  In  this  way,  individual  effects  of  the  various  physical 
parameters  on  the  model's  predictions  can  be  isolated  and  studied.  For 
example,  two  previous  projects  have  dealt  respectively  with  two- 
dimensional  steady-state  isothermal  flows  [2]  and  with  two-dimensional 
thermally  mixing  flows  [3].  This  effort  considers  two-dimensional 
transient  flows  under  isothermal  conditions.  Follow-on  work  has  begun 
on  two-dimensional  thermal  transients  and  three-dimensional  steady 
isothermal  flows. 

1.2  Scope 


In  this  work  a  lab-scale  outlet  plenum  test-cell,  with  internal 
dimensional  ratios  that  simulate  the  Fast  Flux  Test  Facility  (FFTF),  is 


considered  experimentally  and  computationally.  Isothermal  water  flow 
transients  are  produced  in  the  test-cell,  and  Plexiglas  side-wall 
construction  permits  measurement  of  turbulence  field  quantities  throughout 
the  cell  by  a  Laser  Doppler  Anemometer  (LDA). 

The  LDA,  and  its  associated  electronics,  is  arranged  to  measure 
two  orthogonal  velocity  components,  turbulent  kinetic  energy  (TKE),  and 
the  velocity  cross  correlation  (proportional  to  Reynolds  Stress) 
simultaneously  at  any  point  within  the  test-cell.  At  a  number  of 
different  locations,  the  values  of  these  quantities  are  recorded  during  a 
flow  shut-down  transient  that  is  repeated  successively  to  permit  ensemble 
averaging  of  the  recorded  values.  The  result  is  a  record  of  the  mean 
fluid  velocity  components  and  the  turbulence  field  parameters  throughout 
the  test-cell  as  a  function  of  time  during  a  transient.  Details  of  the 
measurement  technique,  data  reduction,  and  the  method  for  achieving  a 
repeatable  flow  transient  are  presented  in  Sec.  3. 

For  the  computational  work,  the  VARR-II  computer  code  is  used  to 
generate  predictions  of  the  same  turbulence  field  parameters  measured 
above.  VARR-II  is  a  two-dimensional  time-dependent  thermal -hydraulic 
code  that  employs  a  unique  two  equation  turbulence  model,  and  was 
developed  for  designer  use  in  the  LMFBR  programs.  Flow  map  "snapshots" 
of  velocity  vectors  and  values  of  the  primary  turbulence  quantities  are 
produced  by  the  code  at  various  times  during  a  numerical  shut-down 
transient  that  is  matched  to  the  experimentally  determined  plenum  inlet 
flows.  These  data  are  then  compared  with  the  experimental  results. 

In  order  to  generate  a  shut-down  transient,  in  either  the 

experimental  or  numerical  case,  an  initial  steady  state  flow  pattern  must 
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be  established.  In  a  previous  work,  Chen  [4]  has  considered  the  steady- 
state  case  In  essentially  the  same  geometry,  with  mixed  results.  He 
suggested  that  departure  from  turbulence  Isotropy  In  the  experimental 
test-cell,  along  with  Insufficient  knowledge  of  detailed  inlet  conditions, 
could  explain  some  of  the  poor  predictions  produced  by  VARR-II. 
Accordingly,  the  test-cell  inlet  has  been  modified  in  this  effort  to 
improve  the  degree  of  isotropy,  flow  synmetry,  and  velocity  profile 
control  in  order  to  determine  their  effect  on  code  performance.  In 
addition  to  the  transient  data,  therefore,  detailed  steady-state 
measurements  were  also  accomplished,  as  were  steady-state  computer  runs. 

A  number  of  numerical  experiments  were  performed  for  the  steady-state 
case  to  investigate  the  effects  of  mesh  size,  boundary  conditions, 
turbulence  model  constants,  and  the  artificial  diffusion  term  that 
accompanies  full  donor  cell  differencing. 


CHAPTER  TOO 


LITERATURE  REVIEW 


2.0  Introduction 

The  basis  for  the  computational  analysis  of  turbulent  flows  is  the 
so-called  "turbulence  model"  that  relates  the  complex  small-scale  motions 
in  the  flow  to  the  large-scale  mean  quantities  that  are  more  readily 
calculated  on  reasonably  sized  machines.  This  review  begins  with  a 
brief  description  of  the  types  and  characteristics  of  available  turbulence 
models.  The  question  of  transient  influence  on  model  performance  is 
subsequently  addressed  as  a  rationale  for  the  experiment  itself.  Finally, 
recent  work  at  Argonne  National  Laboratory  (ANL)  on  a  related  experiment 
and  computational  scheme  is  reviewed  to  put  this  effort  in  context. 

2.1  Turbulence  Models 

Mathematical  predictions  of  isothermal  turbulent  flows  begin  with 
the  well-known  conservation  equations  for  mass  and  momentum.  This  set 
of  partial  differential  equations  (pde)  is  generally  agreed  to  be  "exact" 
in  the  sense  that  its  solution  completely  and  correctly  describes  the  flow 
of  a  single-phase  Newtonian  fluid.  However,  for  turbulent  flows,  this 
equation  set  cannot  be  solved  analytically.  Direct  recourse  to  finite 
difference  formulations  and  digital  computers  proves  impractical  because 
the  exceedingly  small  spatial  scale  of  turbulent  fluctuations  requires 
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a  very  fine  grid  and  the  number  of  mesh  points  needed  for  the  simplest 
problem  exceeds  present  computer  capabilities  by  two  to  three  orders  of 
magnitude  [5  ]. 

On  the  assumption  that  the  important  effects  of  turbulence  may  be 
accounted  for  by  average  or  statistical  measures  of  the  fluctuating  quanti¬ 
ties,  an  approximate  method  for  computing  turbulent  flows  is  commonly 
used.  Following  Reynolds'  classic  time -averaging  of  the  Navier-Stokes 
Equations,  the  instantaneous  velocity  is  decomposed  into  its  mean  and 
fluctuating  parts  and  substituted  into  the  "exact"  equations;  these  equa¬ 
tions  are  ensemble  averaged,  resulting  in  a  new  set  of  partial  differential 
equations  for  the  mean  quantities  [6  ].  Since  all  of  the  dependent  vari¬ 
ables  in  this  set  of  equations  are  averaged  quantities,  they  do  not  vary 
nearly  as  rapidly  over  space  as  the  original  turbulent  fluctuations; 
therefore,  the  problem  of  an  excessively  fine  grid  for  the  finite  differ¬ 
ence  formulation  is  overcome.  However,  the  new  set  of  equations  so 
obtained  is  not  complete,  that  is,  there  are  now  more  unknowns  than 
independent  equations.  The  averaging  process  produces  these  additional 
unknowns  in  the  form  of  correlations  or  mean  products  of  fluctuating 
quantities.  In  the  case  of  isothermal  incompressible  flow,  which  is 
dealt  with  in  this  work,  the  unknown  correlation  that  arises  is  u'.u'.  , 

*  J 

the  so-called  Reynolds  Stress. 

In  effect,  all  of  the  influence  of  the  turbulence  on  the  mean 
motion  is  contained  in  this  term.  In  order  to  close  the  set  of  aver¬ 
aged  equations,  a  relationship  between  u'.u'.  ,  and  one  or  more  of  the 

1  J 

other  mean  quantities  is  required.  The  non-linear  nature  of  the 
applicable  equations  is  such  that  an  exact  mathematical  derivation  of  the 
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required  relation  Is  impossible.  Therefore,  an  equation,  or  set  of  equa¬ 
tions,  must  be  developed  to  specify  the  new  unknowns  in  terms  of  the  other 
variables,  and  thereby  permit  the  solution  of  the  averaged  equation  set. 
These  relational  equations  are  turbulence  models. 

2.1.1  Types  of  Models 

Turbulence  models  are  usually  classified  by  the  number  of  transport 
type  pde's  they  involve.  The  motivation  for  going  to  more  complex  models 
(more  equations)  lies  in  the  hope  of  improved  generality  and  reduced 
reliance  on  empirical  prescriptions  for  solutions  to  engineering  turbulence 
flow  problems.  The  five  major  categories  of  turbulent-flow  pde  models 
are  succinctly  described  by  W.  C.  Reynolds  [7  ]: 

1.  Zero-equation  models  —  models  using  only  the  pde  for  the 
mean  velocity  field,  and  no  turbulence  pde's. 

2.  One-equation  models  --  models  involving  an  additional  pde 
relating  to  the  turbulence  velocity  scale. 

3.  Two-equation  models  --  models  incorporating  an  additional 
pde  related  to  a  turbulence  length  scale. 

4.  Stress-equation  models  —  models  involving  pde's  for  all 
components  of  the  turbulent  stress  tensor. 

5.  Large-eddy  simulations  —  computations  of  the  three- 
dimensional  time-dependent  large-eddy  structure  and  a 
low-level  model  for  the  small-scale  turbulence. 

Launder  and  Spalding  [8  ]  give  a  clear  and  definitive  treatment  of 
these  models,  their  derivation,  advantages,  and  shortcomings.  They  also 
show  that,  for  recirculating  flows,  the  two-equation  model  is  the  simplest 
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that  can  provide  realistic  predictions  of  the  quantities  that  the  designer 
needs  to  know.  Since  the  experiment  in  this  work  involves  recirculating 
flow,  a  two-equation  turbulence  model  is  used  for  the  computer  predictions. 

2.1.2  Two-Equation  Models 

Two-equation  models  attempt  to  determine  turbulent  kinetic  energy 
(TKE  =  K)  and  length  scale  (z)  from  transport  equations.  These  two 
quantities,  in  turn,  specify  the  turbulent  viscosity  (a),  which  relates 
the  unknown  Reynolds  Stress  to  mean  velocity  gradients,  the  basic  goal  of 
all  turbulence  models.  Early  work  in  this  area  has  shown  that  the  turbu¬ 
lence  length  scale,  £,  does  not  diffuse  as  the  gradient  of  z ,  and  so 
transport  equations  for  this  quantity  have  not  produced  successful  models. 
Accordingly,  it  has  become  standard  practice  to  write  the  second  transport 
equation  for  some  other  quantity  that  is  related  to  Z,  usually  either  the 
turbulent  energy  dissipation  rate  (e),  or  directly  in  terms  of  a  itself. 

The  turbulence  model  used  in  VARR-II  takes  the  latter  approach  and  is 
a  so-called  K-a  model  developed  by  J.  H.  Stuhmiller  [9  ].  The  alternative 
K-e  model,  such  as  used  in  the  TEACH-T  code  DO],  was  not  employed 
directly  in  this  work  since  TEACH-T  is  a  steady-state  code.  However,  a 
mathematically  equivalent  model  (fully  described  in  Sec.  4)  was  derived  for 
use  in  VARR-II  and  results  are  reported  for  both  model  types. 

2.2  Transient  Effects 


A  conmon  method  for  generating  the  additional  transport  equations  that 
form  the  turbulence  model  Is  described  by  Bradshaw  [11].  Exact  transport 
equations  for  the  quantities  of  Interest  can  be  formed  by  operating  on  the 
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averaged  momentum  equations.  These  exact  equations  contain  unknown  corre¬ 
lations,  each  of  which  must  be  modelled  In  terms  of  known  mean  variables 
before  the  solution  can  proceed.  This  process,  necessarily  approximate, 
entails  a  degree  of  heuristic  reasoning  based  on  diverse  factors  such  as 
experience,  phenomenological  studies  and  the  mathematical  properties  of 
tensors.  Whichever  closure  assumptions  are  used,  the  resulting  model 
equations  contain  constants,  or  free  parameters,  that  must  be  set  by  experi¬ 
ment.  Empirical  results  used  to  determine  values  for  these  free  parameters 
necessarily  come  from  available  data  largely  derived  from  steady-state 
experiments.  Bradshaw  has  noted  that  “unsteady  flows  present  no  extra 
problems  in  turbulence  modelling  as  long  as  the  total  transport  term  is 
no  larger  than  in  corresponding  steady  flows."  [12]  This  is  clearly  true 
for  the  exact  formulations;  e.g.,  the  transport  equation  for  quantity  5 
is 


3£+  u 
at  ui  axi 


■  production  +  diffusion  -  dissipation 


(2.1) 


total  transport  term 


Strictly  speaking,  the  magnitude  of  the  transport  term,  not  the  relative 
contributions  of  its  components,  determines  the  physics  of  the  situation. 

However,  when  the  unknown  terms  in  the  right-hand  side  of  Eq.  2.1 
are  modelled,  the  free  parameters  are  established  independently  via 
different  experiments.  Generally  the  decay  of  isotropic  turbulence  is 
considered  to  determine  a  portion  of  the  unknown  constants,  and  then  the 
remaining  free  parameters  are  specified  by  homogeneous  shear  flow  experi¬ 
ments  where  the  transport  terms  vanish  completely.  Reynolds  and  Cebeci 
[13]  point  out  that  this  procedure  implies  the  assumption  that  the  same 


empirical  input  suffices  for  both  types  of  flow. 

Of  more  direct  importance  to  this  work  are  Stuhmiller's  closure 
approximations  for  the  K-a  model  used  in  the  VARR  code.  His  unique  develop¬ 
ment  considers  a  single  wavenumber  of  the  turbulence  frequency  spectrum. 
Physical  interpretation  of  the  interactions  with  other  wavenumbers  leads  to 
linearized  forms  of  the  equations  of  motion  for  the  fluctuating  quantities. 
The  concept  of  simulated  fluctuations  at  a  single  wave  number  makes 
explicit  solution  of  the  linearized  equations  possible.  Stuhmiller  suggests 
that  an  advantage  of  this  approach  is  that  it  preserves  the  dynamical 
relation  between  fluctuations  in  addition  to  the  root  mean  square  amplitude 
of  the  correlations.  However,  the  assumption  of  zero,  or  nearly  zero,  time 
derivatives  is  an  essential  feature  of  his  complex  development.  Stuhmiller 
addresses  this  point,  and  partially  motivates  the  current  work,  when  he 
states  that  his  assumption  of  constant  mean  flow  is  an  approximation  that 
must  be  assessed  by  comparison  with  experiment. 

2.3  Related  Outlet  Plenum  Studies 

An  extensive  series  of  experimental  and  analytical  studies  of  thermal 
hydraulic  behavior  in  LMFBR  outlet  plena  is  currently  underway  at  Argonne 
National  Laboratory  (ANL).  The  ANL  reports  on  this  work  are  numerous;  a 
recent  article  by  Howard  and  Carbajo  [14]  references  a  number  of  the 
relevant  papers.  Experiments  performed  in  this  series  were  in  more 
prototypic  (3-D,  cylindrical)  geometries  and  on  a  larger  scale  than  used 
in  the  current  work.  Simultaneous  temperature  and  flow  transients  were 
examined  and  the  effects  of  the  important  dimensionless  groups  studied. 

Based  on  this  work,  calculational  techniques  have  been  developed,  MIX  [15] 
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and  Plenum-3  [16]  for  example,  that  feature  much  shorter  run  times  than 
detailed  fluid  dynamics  codes  such  as  VARR-II. 

MIX  uses  the  vorticity  stream-function  formulation  and  a  zero-equation 
turbulence  model  that  specifies  a  constant  turbulent  viscosity  from 
experiment.  Plenum-3  uses  a  lumped -parameter  approach  that  separates  the 
outlet  plenum  into  three  vertically  stacked  regions  and  solves  the  conserva¬ 
tion  equations  in  each  region.  The  code  predictions  are  generally  assessed 
by  comparing  fluid  temperature  data  at  a  number  of  key  locations  within  the 
plenum.  Agreement  between  codes  of  this  type  and  experiment  has  been  good, 
better  in  fact  than  more  complex,  longer  running  codes,  such  as  VARR-II  [17]. 

The  experiments  performed  in  the  ANL  series  primarily  measured  fluid 
temperatures  and  did  not  provide  detailed  velocity  and/or  turbulence  field 
data.  Overall  evaluation  of  the  more  sophisticated  turbulence  models,  such 
as  that  used  in  VARR,  as  well  as  improvements  in  the  model  free  parameters, 
can  be  more  effectively  accomplished  from  a  data  base  of  specific  turbulence 
field  measurements.  For  example,  a  recent  article  by  N.  Markatos  [18] 
describes  the  application  of  a  K-e  two-equation  model  to  the  transient  flow 
of  liquid  sodium  in  a  geometry  similar  to  that  used  in  the  current  work. 
Turbulence  model  performance  is  not  directly  discussed  in  his  paper. 

However,  the  overall  code  predictions  are  found  to  be  "in  qualitative 
agreement  with  expectations,"  pointing  out  the  sparsity  of  detailed  experi¬ 
mental  data  on  turbulent  recirculating  flows  in  this  geometry. 
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CHAPTER  THREE 
EXPERIMENTAL  WORK 


3.0  Introduction 

Measurements  of  turbulence  field  parameters  (TFP),  :.e.,  two  velocity 

components  (IT  and  V),  turbulent  kinetic  energy  (TKE  *  - ^ — •  *  and 

Reynolds  Stress  (uV),  have  been  made  in  a  plexiglas-walled  test  cell 
that  is  a  two-dimensional  stylized  model  of  an  LMFBR  outlet  plenum.  The 
measurements  were  made  with  a  dual -channel  Laser  Doppler  Anemometer  (LDA) 
which  is  described  in  detail  in  Sec.  3.2.  For  a  typical  measurement 
sequence,  flow  conditions  in  the  test  cell  were  initially  established  at 
a  predetermined  "full"  or  steady-state  flow  rate;  measurement  of  the  TFPs 
were  made  to  determine  the  initial  or  steady-state  values;  then  a  repeat- 
able  flow  transient  was  induced  that  caused  the  inlet  flow  rate  to  decrease 
to  32.5%  of  its  initial  value  in  120  seconds.  During  this  flow  coast-down 
the  values  of  the  TFPs  were  recorded  electronically.  The  same  transient 
measurement  was  repeatedly  performed  at  42  locations  throughout  the  cell  so 
that  ensemble  averages  could  be  calculated  for  all  the  TFPs  during  the  flow 
transient.  At  37  additional  locations,  steady-state  measurements  were 
made  to  provide  a  more  detailed  description  of  the  full-flow  field.  The 
net  result  of  this  process  is  a  series  of  time-sequenced  "snapshot"  plots 
of  velocity  vectors  and  a  record  of  the  values  of  average  TKE  and  Reynolds 
Stress  as  a  function  of  time  during  the  transient.  These  data  are 
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presented  and  discussed  in  Sec.  5. 

This  chapter  describes  the  hydraulic  and  electronic  details  of  the 
experimental  apparatus.  Moreover,  since  turbulence  measurements  in  the 
presence  of  mean  flow  transients  are  not  conmonly  made,  the  specific 
measurement  technique  used  and  its  effect  on  accuracy  are  also  presented. 

3.1  Hydraul ics 

3.1.1  General  Description  of  Loop 

A  schematic  view  of  the  hydraulic  loop  is  shown  in  Fig.  3.1.  The 
loop  is  closed  and  filled  with  distilled  water.  The  filter  and  de-mineral¬ 
izers  are  used  to  maintain  water  purity  and  thereby  to  inhibit  corrosion 
and  provide  a  controlled  concentration  of  scattering  particles  for  use  by 
the  LDA.  Two  pumps  are  mounted  in  parallel  to  provide  a  variety  of  flow 
rates.  The  five  hp  pump  was  used  in  this  experiment  to  obtain  the  maximum 
possible  initial  inlet  Reynolds  No. 

The  heat  exchanger  served  a  dual  purpose.  Initially  it  functioned 
as  a  heat  source  to  raise  the  loop  temperature  to  110°F;  from  that  point  it 
became  a  heat  sink  and  maintained  isothermal  fluid  conditions  by  dissipating 
pump-heat  during  the  prolonged  experimental  runs.  The  parallel -mounted 
rotameters  were  useful  for  providing  gross  indications  of  flow-transient 
repeatability  and  for  monitoring  flow  rates  during  and  between  runs.  A 
more  precise  measure  of  the  repeatability  factor  involved  the  LDA  measure¬ 
ments  themselves  and  is  discussed  in  Sec.  3.1.2.  Each  transient  run  began 
with  both  flow  meters  indicating  100%,  and  at  the  end  of  the  run,  flow  had 
decreased  monotonically  to  32.5%  indicated  on  the  larger  rotameter. 
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3.1.2  Repeatable  Flow  Transient  Generator— (RFTG) 

The  turbulence  parameters  measured  In  this  experiment  were  obtained 
by  constructing  ensemble  averages  (see  Sec.  3.3  for  details)  of  the  time- 
varying  fluid  flow  properties  measured  by  the  LDA.  This  necessarily 
required  a  flow  transient  that  could  be  faithfully  repeated  a  large  number 
of  times.  Further,  the  transient  to  be  measured  had  to  be  initiated  auto¬ 
matically  so  that  the  data  recorder  always  collected  data  from  the  same 
time  base  for  each  run. 

The  technique  selected  is  illustrated  in  Fig.  3.2.  A  two-inch  (air- 
to-open,  spring-to-close)  diaphragm  control  valve  with  a  Moore  Valve 
Positioner/Controller  was  installed  at  the  outlet  of  the  pump  down  stream 
of  the  flow  meters.  During  the  steady-state  (full  flow)  period,  an  air 
signal  is  supplied  via  the  three-way  solenoid  valve  to  the  instrument  port 
of  the  valve  controller.  The  amount  of  valve  travel  is  proportional  to  the 
air  pressure  at  this  port.  By  regulating  this  source  of  air  to  the  valve 
controller  the  steady-state  flow  rate  can  be  repeatably  set  at  any  initial 
value  desired.  This  initial  charge  of  air  was  accurately  checked  for 
repeatability  with  a  Heise  pressure  gauge  throughout  the  experimental 
period.  At  the  start  of  the  transient,  the  solenoid  is  energized,  and  air 
is  permitted  to  slowly  bleed  out  of  the  valve  controller  instrument  port 
(and  parallel -connected  surge-tank)  through  a  micrometer-adjusted  needle 
valve.  The  micrometer  setting  determines  the  rate  of  air  pressure  decrease 
at  the  instrument  port  and  thereby  the  valve  closure  speed,  and  ultimately, 
the  fluid  flow  rate.  This  simple  system  fulfills  the  automatic  requirement, 
and  has  been  tested  for  repeatability  with  excellent  results. 


Two  methods  were  used  to  check  and  maintain  the  flow- transient 
repeatability  throughout  the  data  collection  period.  The  flow-meters  were 
one  source  of  information  for  determining  repeatability  of  initial  and 
final  flow  rates.  A  more  accurate  determination  was  made  at  each  measure¬ 
ment  point  by  using  the  output  voltage  of  the  LDA  Frequency  Trackers  to 
insure  that  the  initial  velocity  at  a  point  was  the  same  (within  1%)  as 
the  previous  run.  This  was  particularly  useful  after  an  overnight  shut¬ 
down.  In  such  a  case,  the  LOA  measuring  volume  was  not  moved  from  the  last 
point  measured  until  the  loop  was  fully  warmed  up  and  the  steady-state 
velocity  at  that  point  had  stabilized  to  the  same  value  (within  1%)  that 
had  been  recorded  on  the  previous  day. 

The  repeatability  of  the  shape  of  the  flow  transient  as  a  function  of 
time  was  ascertained  by  visual  readings  of  the  flow  meters  during 
repetitive  valve  close-downs  with  a  fixed  needle-valve  setting  for  various 
rates  of  closure.  Flow  rates  were  recorded  at  ten-second  intervals  during 
a  number  of  identical  runs;  the  average  value  at  each  sample  time  was  then 
used  to  form  a  plot  of  mean  flowrate  vs.  time.  The  individual  flowrates 
at  each  point  in  time  formed  a  population  whose  standard  deviation  was  less 
than  2%  of  the  mean  at  points  of  worst  agreement,  thereby  indicating  a 
reliable  replication  of  transient  shape.  Further,  during  the  many  hundreds 
of  transients  repeated  during  the  period  of  experimental  data-taking,  key 
flowrate  values  were  always  observed  to  occur  at  the  same  time  during  the 
transient. 

3.1.3  FFTF  Test  Cell 

The  test  cell  used  in  this  experiment  is  a  stylized  two-dimensional 


model  of  an  IMFBR  outlet  plenum  designed  to  exhibit  the  recirculating  flow 
and  turbulent  mixing  characteristics  of  that  geometry.  Measurements  made 
In  this  test  cell  during  a  120  second  coast-down  transient  are  compared  to 
VARR-II  computer  predictions  In  Sec.  5. 

The  FFTF  test  cell  is  a  1/15  scale  model  of  a  two-dimensional  diametral 
section  of  the  Fast  Flux  Test  Facility  (FFTF)  reactor  outlet  plenum  (Fig. 
3.3).  It  Is  constructed  of  a  machined  aluminum  frame,  two  inches  in  depth, 
with  plexiglas  face-plates  to  permit  the  transmission  of  light  for  the 
Laser  Doppler  Anemometer.  This  cell  Is  essentially  the  same  device  used 
by  Chen  for  steady-state  experiments,  and  reported  on  in  reference  [19]. 
However,  discrepancies  in  Chen's  data  (discussed  below)  suggested  a  number 
of  modifications  to  the  cell,  primarily  in  the  inlet  region.  The  following 
describes  those  modifications  and  their  effects  on  the  data. 

Since  measurements  and  calculations  are  carried  out  in  only  half  the 
cell,  flow  symmetry  across  the  cell  center-line  must  be  obtained  initially 
and  maintained  throughout  the  transient.  Lack  of  flow  symmetry  in  the 
unmodified  test  cell  used  by  Chen  is  the  most  likely  reason  that  cell-scale 
mass  balances  obtained  from  his  data  are  not  consistent  .  The  modifications 
needed  to  alleviate  this  problem  were  found  to  be: 

(a)  change  from  one  inlet  channel  to  three  (each  independently 
valved); 

(b)  installation  of  throttling  valves  on  each  outlet; 

(c)  installation  of  four-level  rigid  grids  and  a  shaped 
entrance  at  the  bottom  of  the  chimney. 

These  modifications  are  sketched  and  noted  in  Fig.  3.3.  The  three 
inlet  channels,  each  valved  separately,  markedly  reduce  the  center-line  jet 
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velocity  experienced  in  the  original  cell.  The  outlet  valves  permit 
balancing  of  the  flow  resistances  in  the  two  outlet  channels.  The  shaped 
entrance  region  eliminates  the  sharp-edged  entrance  effects  of  the  original 
cell,  and  the  four  layers  of  stainless  grid  suppress  the  turbulence 
intensity,  further  flatten  the  velocity  distribution,  and  produce  nearly 
isotropic  turbulence  at  the  cell  inlet. 

The  improvement  in  flow  symmetry  obtained  by  these  changes  is  indi¬ 
cated  by  cell-scale  vertical  mass  balances  that  were  made  at  each 
horizontal  plane  of  measurement  in  the  test  cell  at  full  flow  conditions. 
With  perfect  flow  symmetry  and  truly  two-dimensional  flow,  exact 
correspondence  between  upflow  and  downflow  should  obtain  at  each  horizontal 
plane  listed  in  Table  3.1.  Overall  agreement  is  good,  significantly 
imoroved  from  the  unmodified  case,  but  severe  deterioration  of  the  mass 
balance  is  still  observed  near  the  top  of  the  cell.  This  discrepancy  is 
likely  caused  by  three-dimensional  effects  that  become  more  important  in  the 
upper  region  of  the  cell.  Velocity  measurements  off  the  cell  center-plane 
show  that  the  missing  downflow  is  being  pushed  outside  of  the  normal 
measurement  plane  in  this  region  (see  Appendix  E  for  details).  As  a  result, 
velocity  predictions  by  strictly  two-dimensional  fluid  dynamics  codes,  such 
as  VARR-II,  would  not  be  expected  to  conform  to  the  measurements  made  in 
this  part  of  the  test  cell.  Although  these  three-dimensional  effects  appear 
to  be  unavoidable  at  the  top  of  the  cell,  cross-wise  velocity  measurements 
and  the  mass  balance  data  suggest  they  are  far  less  important  over  the  main 
body  of  the  recirculating  region. 

Another  area  of  concern  regarding  Chen's  data  is  that  the  highly- 
peaked  velocity  distribution  at  the  inlet  exhibited  an  inflection  point 
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indicating  that  the  flow  was  not  fully  developed.  This  condition  contri¬ 
buted  to  the  significant  anistropy  of  the  turbulence  in  the  inlet  region, 
and  has  been  remedied  by  the  modifications  performed  on  the  new  cell. 

Since  the  turbulence  model  used  in  VARR  *s  based  on  the  assumption  of 
roughly  isotropic  turbulence,  the  improved  isotropy  of  the  modified  inlet 
can  be  expected  to  yield  measurements  closer  to  the  code  predictions.  This 
was  found  to  be  the  case  and  is  discussed  more  fully  in  Sec.  5.  It  should 
be  noted  that  truly  isotropic  turbulence  can  rarely  be  expected  in  real 
engineering  flow  situations.  As  such,  it  is  usually  preferable  that  a 
hydrocode  of  this  type  not  exhibit  extreme  sensitivity  to  inlet  turbulence 
levels. 


A  dual-channel  Laser  Doppler  Anemometer  (LDA),  operating  in  the 
reference-beam  mode,  was  used  to  make  the  measurements  of  instantaneous 
velocity  components  in  this  experiment.  The  dual-channel  designation 
indicates  that  the  LDA  performed  simultaneous  measurements  of  two 
orthogonal  velocity  components  at  each  point  in  the  test  cell.  The 
reference-beam  mode  refers  to  the  optical  configuration  used  in  the  set-up 
and  is  discussed  in  what  follows.  The  equipment  and  optical  system  used 
in  this  experiment  are  basically  the  same  as  that  used  by  Chen  and  reported 
on  in  detail  in  reference  [20],  However,  some  important  modifications  were 
made  to  the  laser  and  to  the  signal  processing  equipment  that  improved 
overall  system  performance,  and  these  are  described  in  detail  throughout 
Sec.  3.2.3. 


3.2.2  Theory  of  Operation 


The  basic  principle  of  the  Laser  Doppler  Anemometer  is  that  the 
frequency  of  light  observed  emanating  from  a  moving  source  will  be  changed 
depending  on  the  velocity  of  the  light-source  relative  to  the  observer. 

In  the  same  fashion,  the  familiar  "Doppler"  effect  causes  an  apparent 
change  of  pitch  (or  frequency)  in  the  whistle  of  a  fast  train  that  is  in 
relative  motion  to  an  observer.  Yeh  and  Cumnins  [21]  first  demonstrated 
the  feasibility  of  using  the  LDA  to  measure  fluid  velocities  in  the 
laboratory  in  1964.  Since  that  time  the  technique  has  been  studied, 
refined  and  comnerciallzed.  As  a  result,  manufactured  LDA  systems  are 
readily  available  in  modular  form  that  permit  relatively  easy  alignment 
and  straightforward  operation.  A  detailed  discussion  of  the  theory  of 
operation  of  the  LDA  is  available  from  many  good  sources,  e.g.  [22],  [23], 
or  [24],  Chen's  report  [25]  also  includes  a  general  synopsis  of  the  theory. 
Such  details  will  not  be  repeated  here,  but  the  major  advantages  of  the  LDA 
over  conventional  fluid  flow  measuring  devices  such  as  the  hotwire  or 
pitot  tube  bear  repeating: 

(1)  the  measurement  is  performed  with  laser  beams;  no  probes 
disturb  the  flow; 

(2)  the  measuring  volume  is  small  enough  to  provide  a  high 
degree  of  spatial  resolution; 

(3)  no  calibration  is  required; 

(4)  fluid  velocity  Is  linearly  proportional  to  the  Doppler 
shift; 

(5)  the  directional  sensitivity  is  Ideal  for  simultaneous 
measurements  in  two  orthogonal  directions. 


The  major  disadvantage  of  the  LDA  in  an  experiment  of  this  nature  is  the 
possible  failure  to  account  for  a  number  of  sources  of  error  that  can 
modify  the  raw  system  output  signal  significantly.  These  factors  are 


addressed  in  Appendix  B,  "Error  Analysis." 

3.2.3  Experimental  Arrangement 

The  basic  components  of  the  LDA  used  in  this  experiment  are: 

(1)  light  source  (laser) 

(2)  light  scattering  particles 

(3)  system  to  focus  and  collect  the  light  (optics) 

(4)  method  of  monitoring  the  Doppler  shift  or  frequency 
(frequency  tracker) 

These  components  are  described  below  in  general  and  treated  in  detail  only 
insofar  as  they  differ  from  the  set-up  used  by  Chen  [26]. 

3.2.3. 1  Laser 

A  Spectra-Physics  Model  164  Argon  Ion  laser,  capable  of  a  maximum 
power  of  2  watts,  was  used  as  the  light  source.  The  laser  was  tuned  to  the 

O 

green  line  at  a  wavelength  of  5145  A. 

A  Spectra  Physics  Model  589  Air  Spaced  Etalon  was  installed  in  the 
cavity  of  the  Argon  Laser  system.  The  purpose  of  the  etalon  is  to  inhibit 
longitudinal  oscillations  in  all  but  one  laser  frequency.  Operation  in 
the  single  mode  permitted  by  the  etalon  results  in  a  much  improved  coherence 
length  (Lc)  for  the  laser.  A  long  coherence  length  is  required  for  LDA 
operations  when  split-beam  path  lengths  are  not  exactly  equal,  which  is 
the  case  for  the  dual  channel  measurements  being  performed  in  this 
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experiment.  For  the  DISA  optical  unit  configured  for  dual  channel  measure¬ 
ments,  optical  path  length  differences  (AL)  on  the  order  of  5-6  cm.  occur. 
This  length  is  unacceptably  close  to  the  10-20  cm.  coherence  length  of  the 
Argon  laser  when  operated  without  single  frequency  etalon.  Installation  of 
the  etalon  increases  the  coherence  length  to  about  10  meters,  thus 
satisfying  the  requirement  that  Lc  »  AL  for  the  optimal  signal -to-noise 
ratio  [27]. 

3. 2. 3. 2  Light  Scattering  Particles 

The  theory  of  operation  of  the  LDA  states  that  the  Doppler  fre¬ 
quency  measured  is  specifically  associated  with  a  light-scattering 
particle  and  actually  yields  the  particle  velocity  rather  than  the  fluid 
velocity.  It  is  therefore  essential  that  the  scattering  particles  used  in 
a  particular  flow  situation  follow  the- fluid  motion  directly,  or  that  the 
velocity  of  the  scattering  particles  can  be  related  to  the  fluid  velocity 
in  a  reliable  fashion.  In  turbulent  water  flows,  "natural"  seeding  is 
generally  acceptable  for  most  LDA  applications.  That  is,  the  microscopic 
organisms  and  suspended  matter  normally  present  in  ordinary  tap  water  are 
used  as  light  scattering  particles.  As  mentioned  in  Sec.  3.1.1,  tap  water 
was  not  used  in  the  experiment  because  of  corrosion  problems,  and  so 
artificial  seeding  was  required. 

The  properties  of  the  seed  that  determine  its  ability  to  closely 
follow  the  water  velocity  are  its  shape,  size,  and  relative  density  ratio. 
Reference  [28]  reviews  various  methods  of  calculating  the  motion  of  par¬ 
ticles  in  a  turbulent  fluid  flow.  The  authors  demonstrate  that  smooth, 
small  particles  of  density  equal  to  the  fluid  are  best.  In  this  experiment. 
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polystyrene  latex  spheres  (diameter  s  .5  urn,  specific  gravity  =  1.05), 
which  meet  all  three  criteria,  are  used  and  can  be  safely  assumed  to  follow 
the  velocity  fluctuations  exactly. 


Light  scattering  ability  is  a  second  criteria  to  be  considered. 

Although  the  individual  particles  are  quite  small,  reference-beam  operation 
achieves  its  highest  signal -to-noise  ratio  for  high  scattering  particle 
concentration,  which  places  many  particles  in  the  measuring  volume  at  any 
given  time.  This  high  particle  concentration,  combined  with  the  excess 
laser  power  available,  provided  acceptable  signals  throughout  the  experiment. 

3.2. 3. 3  Optical  Arrangement 

The  three  basic  modes  of  LDA  optical  arrangement  refer  to  the 
method  by  which  the  Doppler  frequency  shift  is  amplified  and  detected. 

They  are: 

(1)  reference  beam  mode; 

(2)  dual  beam  or  fringe  mode; 

(3)  single  beam  or  two  scattered  beam  mode  [29]  . 

The  dual  beam  mode  seems  to  be  the  most  generally  used  method  for  water 
flows  with  cormercially  available  LDA  systems.  However,  the  requirement 
to  measure  two  velocity  components  simultaneously,  while  achievable  with 
the  dual  beam  or  fringe  method,  was  most  easily  and  directly  met  with  the 
reference  beam  mode.  Further,  the  excess  power  available  from  the  Argon 
laser  made  it  possible  to  get  good  signals  with  the  high  particle 
concentrations  needed  for  reference  beam  operations. 

The  reference  beam  arrangement  used  in  this  experiment  is  illustrated 
schematically  in  Fig.  3.4.  The  original  laser  beam  is  divided  by  a  90/10 
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beam-splitter  and  the  weak  beam  is  directed  to  the  acousto-optic  Bragg  Cell. 
This  beam  is  shifted  in  frequency  by  40  MHz  in  the  Bragg  Cell  and  then 
passes  through  an  adjustable  neutral  density  filter  that  allows  external 
control  of  its  intensity.  The  frequency-shifted  beam  is  split  again,  and 
the  two  resulting  beams  are  made  parallel  to  the  unshifted  beam  and  arranged 
so  that  an  end  view  of  the  beams  places  the  three  spots  on  three  corners  of 
a  square.  The  two  frequency-shifted  beams,  reduced  in  intensity,  are  on 
opposite  ends  of  the  diagonal;  these  are  the  reference  beams.  The  third, 
high  intensity,  beam  is  referred  to  as  the  scattering  beam.  This  three- 
beam  group  is  then  focused  by  a  lens  (focal  length  *  300  mm.),  and  their 
point  of  mutual  intersection  defines  the  measuring  volume.  The  scattering 
particles  in  the  flow  pass  through  the  measuring  volume,  are  illuminated 
by  the  bright  scattering  beam,  and  the  light  reflected  from  these  moving 
particles  is  shifted  in  frequency  depending  on  the  velocity  of  the 
particle.  This  Doppler  shifted  light  is  of  low  amplitude  and  the  extent 
of  the  frequency  change  is  very  small  when  compared  to  the  original  laser 
frequency.  In  other  words,  the  Doppler  frequency  cannot  be  readily  deduced 
from  this  weak  signal  alone.  A  simplified  description  of  the  process  by 
which  this  signal  is  amplified  and  made  useable  in  the  reference  beam  mode 
follows. 

Recall  that  this  weak  signal  is  produced  by  particles  moving  with 
the  flow  that  are  illuminated  by  the  high-intensity  scatter  beam.  The 
reference  beam  for  each  of  the  two  orthogonal  velocity  components  is  much 
weaker  than  the  scatter  beam  and  light  scattered  by  particles  from  this 
beam  does  not  contribute  to  the  Doppler  shifted  light  reaching  the  photo¬ 
multiplier  tube  (PMT)  detector.  However,  the  reference  beam  is  aligned 
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directly  on  the  axis  of  the  PMT  and  thereby  produces  a  non-Doppler-shifted 
signal  whose  amplitude  can  be  adjusted  by  the  external  neutral  density 
filter  control.  Therefore,  two  light  signals  arrive  at  the  PMT: 

(1)  A  cosajjt  -  Weak  amplitude  Doppler  shifted  light 

reflected  from  particles  passing  through 
the  strong  scattering  beam,  where 

=  w0  +  wd  (sign  indicates  velocity  sense), 

ajg  =  Doppler  frequency, 

to q  *  original  laser  light 

(2)  8  cosu)  2t  -  Strong  amplitude  signal  produced  by  reference 

beam  aimed  directly  at  PMT,  where 

oj2  *  Wq  -  40  MHz,  because  of  Bragg  cell 
induced  frequency  shift. 

The  PMT  is  a  square-law  detector.  The  output  of  an  ideal  square-law 
detector  resulting  from  the  combination  of  the  two  signals  described  above 
is: 

I(t)  3  (A  cosu^t  +  B  cosi^t)2 

?  2  2  2 
=  A  cos  (0 1 1  +  AB[cos(o)^+uL>2)t  +  C0s(a)^-a)2)t]  +  B  COS  U)2t. 

(3.1) 

If  this  signal  is  low-pass  filtered  to  remove  frequencies  above  (u^-o)2),  or 
if  the  electronic  equipment  cannot  respond  to  frequencies  in  the  optical 
range— basically  an  equivalent  phenomenon— the  resulting  PMT  output  becomes: 
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I ( t)  *  A2  +  B2  +  AB  costa^H,  (3.2) 

where 

-  o>2  =  ooq  +_ Up  -  (cl)q  -  40)  =  40  +  ojq 

The  constant,  or  d.c.,  part  of  the  signal  can  be  blocked,  and  only  the 
frequency  (40  +  u^) ,  amplified  by  the  strong  signal  B,  remains. 

This  is  the  heterodyne  principle  that  the  reference  beam  mode  exploits. 
This  relatively  strong  signal  is  then  directed  to  the  frequency  trackers 
for  further  processing  (Sec.  3.3).  It  is  not  necessary  for  the  reference 
beam  to  pass  through  the  measuring  volume,  but  this  procedure  reduces 
vibration  problems  and  makes  for  easier  alignment. 

The  component  of  particle  (or  fluid)  velocity  that  is  measured  is  in 
the  plane  formed  by  the  reference  and  scatter  beams  and  perpendicular  to  the 
bi-sector  of  their  intersection  angle  [30].  As  a  result,  the  three-beam 
arrangement  described  above,  with  two  PMTs,  permits  measurement  of  two 
orthogonal  velocity  components  simultaneously. 

3.2. 3.4  Point  Positioning 

The  measuring  volume  location,  or  X-Y  coordinate  position  in  the 
test  cell,  is  controlled  by  fixing  the  test  cell  and  moving  the  laser  and 
associated  optics  as  an  integral  unit.  This  is  accomplished,  as  shown  in 
Fig.  3.5,  by  anchoring  the  laser,  optics  and  PMTs  to  a  milling-machine 
base  which  provides  precise  3-axis  maneuverability.  While  this  arrangement 
makes  alignment  problems  more  tedious,  fixing  the  location  of  the  test 
cell  is  preferred  because  changing  the  position  of  the  flexible  hoses 
that  are  attached  to  the  test  cell  can  affect  the  flow  pattern  in  the  cell. 
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Measurements  were  made  at  the  various  X-Y  coordinate  locations  on 
the  cell  center-plane.  Because  of  index  of  refraction  effects,  the  center- 
plane  was  located  by  moving  the  measuring  volume  (beam  intersection)  from 
the  inside  surface  of  one  plexiglas  face  to  the  other,  recording  the 
translation  distance,  and  then  moving  back  one-half  this  distance. 

3.3  Signal  Processing 

3.3.1  General  Description 

Frequency  tracking  was  used  as  the  method  for  converting  the  PMT 
signal  into  a  continuous  voltage  that  is  proportional  to  fluid  velocity. 
Simultaneous  measurement  of  two  velocity  components  requires  two  independent 
frequency  trackers.  The  DISA  model  55L20  and  TSI  model  1090  trackers  were 
used,  which  allowed  a  side-by-side  comparison  of  the  two  systems.  Although 
both  units  are  frequency  trackers,  they  differ  in  their  fundamental  method 
of  following  the  varying  Doppler  frequency.  The  DISA  system  uses  a  phase 
discriminator  [31],  while  TSI  incorporates  a  phase-locked  loop  to  perform 
this  function  [32].  Further,  the  two  trackers  differ  somewhat  in  manual 
control  features,  input  filtering,  and  output  signal  characteristics. 

For  this  turbulent  water  flow,  using  the  reference  beam  mode,  the  DISA 
tracker,  with  more  frequency  ranges,  a  finer  Voltage-Control led-Oscillator 
frequency  selection  capability,  and  a  narrow-band  input  filter  was  found 
to  have  more  flexibility  and  a  less  noisy  output  for  most  measuring 
situations  encountered.  However,  the  TSI  tracker  was  able  to  maintain 
track  on  the  poorest  signals,  when  DISA  could  not,  and  with  proper  filter¬ 
ing  of  the  output  signals,  the  electronic  noise  problem  could  be  overcome. 
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A  schematic  diagram  of  the  analog  portion  of  the  electronic  layout  is 
shown  in  Fig.  3.6.  As  indicated  in  the  sketch,  the  two  trackers  used  a 
common  40  MHz  oscillator.  In  Sec.  3. 2. 3. 3  it  was  shown  that  the  Doppler 
frequency  is  contained  in  a  signal  from  the  PMT  whose  frequency  is 
(40  +  oup)MHz.  The  "40"  in  this  expression  arises  from  the  nominal  40  MHz 
frequency  shift  imposed  on  both  reference  beams  by  the  single  Bragg  Cell  in 
the  optical  unit.  Since  the  Bragg  Cell  is  driven  by  a  single  40  MHz 
oscillator,  maximum  accuracy  is  obtained  if  that  same  signal  drives  the 
local  oscillator  which  is  used  in  each  tracking  system  to  extract  ojd  from 
(40  +  uD).  Therefore,  the  two  tracking  systems  were  modified  to  use  the 
identical  40  MHz  driver  for  their  individual  local  oscillators. 

3.3.2  Low  Pass  Filtering 

The  phase-locked  loop  configuration  of  the  TSI  tracker  and  its 
larger  IF  bandwidth  (15%  fixed  on  TSI  vs.  .5%  to  8%  variable  on  DISA)  is 
thought  to  be  the  reason  that  the  TSI  system  could  better  maintain  track  on 
rapidly  varying  signals  of  poor  quality.  However,  although  the  average 
value  of  the  frequency  tracked  in  these  cases  was  the  correct  average 
Doppler  frequency,  the  tendency  to  track  noise  within  the  IF  bandwidth 
range  was  apparent.  As  a  result,  the  output  voltage  of  the  TSI  tracker 
required  low-pass  filtering  to  remove  these  high-frequency  noise  contribu¬ 
tions  prior  to  making  the  mean  square  and  correlation  measurements.  Though 
most  important  when  operating  on  the  middle  frequency  range,  the  low-pass 
filter  was  installed  for  all  measurements. 

The  filter  used  was  designed  and  built  in-house  as  part  of  the 
project.  It  is  an  active  dual -channel  state-variable  low-pass  filter 
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with  precision  adjustable  gain,  damping,  and  cut-off  frequency.  Both 

Q 

channels  are  buffered  to  provide  an  input  impedance  of  10  ohms,  and  the 
filtering  characteristics  are  exactly  matched  so  that  the  outputs  from 
both  trackers  are  treated  equally.  In  water  flows,  very  little  turbulent 
energy  is  contained  in  fluctuations  above  500  Hz;  consequently,  the  low-pass 
cut-off  was  conservatively  set  at  1  KHz.  This  value  effectively  filtered 
out  the  high  frequency  noise  while  preserving  the  real  turbulent  energy 
contribution  to  the  mean  square  measurements. 

The  adjustable  gain  feature  of  the  filter  was  crucial  to  the 
transient  measurements  of  TKE.  Since  it  is  not  possible  to  change  the 
range  switch  on  the  averaging  RMS  voltmeters  in  the  midst  of  a  transient, 
it  is  necessary  that  the  initial  RMS  value  register  near  the  top  of  the 
particular  range  selected  so  that  its  value  can  be  accurately  recorded 
throughout  the  transient.  Once  the  steady-state  RMS  value  of  the  signal 
from  each  channel  had  been  established,  the  filter  gain  controls  were 
manipulated  to  bring  the  values  near  the  top  of  their  respective  ranges  to 
give  maximum  coverage  and  accuracy  during  the  transient. 

3.3.3  Analog  Instrumentation 

The  output  signal  from  each  tracker  is  a  fluctuating  voltage  whose 
average  value  (or  d.c.  component)  is  proportional  to  the  mean  fluid 
velocity  in  the  component  direction  associated  with  that  tracker  (O’  and  V). 
The  mean  square  value  of  the  fluctuating  part  of  this  output  signal, 
corrected  for  broadening  errors  as  discussed  in  Appendix  B,  is  twice  the 
turbulence  kinetic  energy  (TKE)  associated  with  a  particular  direction 
— r  — 7 

( u * 2  and  v'  ).  Finally,  the  averaged  product  of  the  two  fluctuating 


components  is  the  Reynolds  stress  (uV). 

Standard  conmercial  averaging/integrating  devices  were  used  to  produce 
the  five  time-averaged  quantities  discussed  above.  The  hardware  layout 
is  as  shown  in  Fig.  3.6,  and  a  complete  list  of  the  electronic  equipment 
is  given  in  Table  3.2. 

In  addition  to  a  digital  or  meter  read-out  of  the  measured  quantities, 
each  device  produces  an  analog  voltage  that  is  proportional  to  the  time- 
average  being  performed.  These  five  analog  signals  were  digitized  and 
recorded  during  each  transient.  The  characteristics  of  these  signals,  their 
relationship  to  true  mean  values,  and  the  selection  of  the  time  constant 
for  the  required  integration  is  discussed  in  Sec.  3.4. 

3.4  Transient  Data  Collection  and  Analysis  Method 

3.4.1  General 

The  measurement  of  mean  turbulence  parameters  during  a  flow 
transient  presents  the  problem  of  correctly  defining  the  quantities  actually 
measured.  The  basic  assumption  underlying  the  calculational  prediction 
of  turbulent  flows  is  that  equations  for  average  values  of  the  fluctuating 
turbulence  quantities  adequately  describe  the  flow.  The  predictions  of  a 
typical  fluid  dynamics  code,  such  as  VARR,  are  given  in  terms  of  average 
or  mean  values.  Therefore,  for  purposes  of  comparison,  the  quantities  of 
interest  to  be  obtained  from  an  experiment  must  also  be  average  values. 

The  problem  in  the  case  of  a  transient  experiment  is  that  the  very 
existence  of  the  transient  makes  the  normal  method  of  forming  mean  values 
from  turbulence  data  invalid. 


This  difficulty  can  be  clearly  seen  by  considering  the  definitions  of 


the  two  common  types  of  averages  applicable  to  turbulence: 
(1)  Time-Average 

T+t 

Ut  =  Limit  ^  ^  u ( t ' } dt ' 

-T+t 


(2)  Ensemble  or  "Probability"  Average 
N 


l  U-(t) 

Ue  =  -£l— 


+  00 


N 


i. 


U(t)p(U)dU  , 


where  p(U)  is  the  normalized  probability  density  function  which  satisfies: 


J  p(U)dU  *  1 


Also  note,  in  the  theory  of  turbulence,  Ut  =  Ug,  the  so-called  "ergodic 
hypothesis." 

In  the  more  common  experimental  situation  of  steady-state  turbulent 
flow,  the  appropriate  average  values  are  obtained  by  method  (1),  the  time- 
average  technique.  For  the  steady-state  case,  U"t  is  a  function  of  T  and  not 
of  t.  The  integration  time,  T,  is  taken  large  enough  so  that  Ft  is  constant. 
This  method  is  preferred  to  the  ensemble  or  probability  approach  because 
it  avoids  the  very  large  number  of  experimental  repetitions  (N)  and  does  not 
require  that  p(U)  be  known. 

However,  because  of  the  requirement  that  be  constant,  the  time- 
average  method,  strictly  applied,  is  not  appropriate  for  a  transient 
experiment  where  all  mean  turbulence  quantities  are  expected  to  vary  with 
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time.  Also,  since  p(U)  is  not  known,  the  probability  average  is  not 
useful.  The  ensemble  average,  therefore,  is  the  correct  method  for  cases 
where  the  average  value  of  a  fluctuating  quantity  may  itself  be  a  function 
of  time.  Indeed,  though  Reynolds'  classic  time-averaging  of  the  Navier- 
Stokes  equations  recasts  them  into  a  more  computationally  tractable  form, 
strictly  speaking,  it  is  the  ensemble-averaged  forms  of  these  equations  that 
are  solved  when  the  time  derivatives  are  not  identically  zero. 

In  principle,  application  of  the  ensemble  approach  is  straightforward. 
However,  for  transient  measurements  at  fifty  or  so  different  locations, 
the  total  number  (N)  of  transient  runs  required  to  form  statistically  valid 
averages  would  be  impractical ly  large.  Further,  the  speed,  cost  and  data¬ 
storage  capability  of  the  recording  equipment  required  for  such  an  effort 
are  discouragingly  high. 

Accordingly,  a  method  of  approximating  the  true  averages,  similar  to 
that  used  for  the  measurement  of  turbulence  in  atmospheric  and  ocean  flows 
[33],  was  devised.  In  such  cases,  time-averages  are  made  with  averaging 
times  (T)  that  are  large  with  respect  to  the  period  of  turbulent  fluctua¬ 
tions,  but  small  when  compared  to  the  time-scale  of  the  diurnal  or  tidal 
variations  which,  strictly  speaking,  cause  the  mean  values  so  obtained  to 
be  functions  of  time.  A  typical  example  occurs  in  the  measurement  of  ocean 
flows  where  turbulent  components  on  the  order  of  1  Hz  are  observed  super¬ 
imposed  on  tidal  variations  which  exhibit  periods  of  12  to  24  hours.  Appli¬ 
cation  of  time-averaging  in  this  case  gives  a  mean  value  which  is 
essentially  constant  during  the  averaging  period. 

Mathematically,  the  criteria  to  be  met  for  successful  application  of 


the  above  approximation  is: 


45 


Hurb  <<:  ^int  <<  ^trans; 


where 


tturb  3  Time  scale  of  turbulent  fluctuations, 
Tint  *  Integration  time. 


Hrans 


*  Time  scale  of  flow  transient. 


The  objective  in  this  experiment  was  to  minimize  t 


trans 


so  that  time  deriva 


tives  and  their  effects  on  the  flow  would  be  enhanced.  This  constraint  led 


to  a  T1-nt  which,  though  considerably  larger  than  tturfc),  was  not  so  large  as 
to  produce  a  smooth,  pseudo-constant  U"t(t),  and,  naturally,  was  small 
enough  to  follow  the  imposed  transient  to  within  a  known  small  error.  This 
time-averaged  signal  was  recorded  for  a  series  of  identical  transients,  and 
then  a  single  set  of  the  five  ensemble-averaged  functions  of  time  was 
formed  from  this  data. 


The  averaging  method  employed,  therefore,  was  a  combination  of 
ensemble  and  time-averaging  techniques.  The  time-averaging  served  to 
decrease  the  frequency  and  amplitude  of  the  fluctuating  signal  to  be 
digitized;  this  permitted  the  use  of  a  slower  and  less  expensive  data 
collection  system,  as  well  as  a  reduction  in  the  number  of  transient  runs 
needed  to  get  a  statistically  suitable  ensemble  average.  Six  to  seven 
transient  runs  at  each  position  generally  yielded  standard  deviations 
that  did  not  significantly  improve  with  additional  runs. 

Use  of  this  pseudo-ensemble  averaging  technique  entails  a  degree  of 
error  that  depends  on  a  number  of  factors.  These  are  discussed  and 
analyzed  in  detail  in  Sec.  3,4.3.  A  schematic  illustration  of  the  method 
is  shown  in  Fig.  3.7. 


3.4.2  A-to-P  Conversion  and  Data  Recording 


The  five  analog  signals  depicted  as  average  outputs  in  Fig.  3.6  are 
actually  time-varying  analog  voltages  (e.g.  ( t )  in  Fig.  3.7).  In  order  to 

form  the  ensemble  averages  described  above,  these  voltages  must  be  recorded 
for  a  series  of  identical  transient  runs,  and  then  the  values  at  each  point 
in  time  mathematically  averaged  to  give  the  final  output  (e.g.  Ue( t )  in 
Fig.  3.7). 

This  procedure  was  accomplished  by  a  five  channel  A-to-D  converter 
that  sampled  and  digitized  all  five  outputs  at  ten  second  intervals  through¬ 
out  each  transient  run.  The  digital  values  were  transferred  to  magnetic 
tape  at  the  experimental  site  for  later  processing  by  a  previously-available 
Tektronix  4051  mini -computer.  The  Tektronix  4924  cassette  tape  recorder, 
fully  compatible  with  both  the  4051  and  the  in-house  A-to-D,  was  used  for 
data  storage. 

The  A-to-D  was  designed  and  constructed  in-house  as  part  of  the 
project.  This  permitted  tailoring  of  the  five  input  channels  to  the  pre¬ 
cise  requirements  of  the  experiment,  and  provided  a  high  degree  of  accuracy 
and  flexibility  at  a  relatively  low  cost.  General  features  of  the  device 
include: 

(1)  Digitization  of  1  to  5  channels  simultaneously; 

(2)  8-bit  A-to-D  conversion; 

(3)  Variable  sampling  time  (.1  sec.  to  25.5  sec.); 

(4)  Buffered  inputs  capable  of  digitizing  the  voltage  ranges 
0  to  10,  3  to  5,  4  to  6,  0  to  -1,  and  -1  to  +1  with 
maximum  accuracy; 

(5)  Input  buffer  amp  accuracy  better  than  +  0.1%; 


(6)  External  Programming  Capability. 


A  complete  description  of  the  design  and  operation  of  the  A-to-D  converter 
is  included  as  Appendix  0. 


3.4.3  Averaging  Method  Error  Analysis 

The  use  of  the  special  averaging  technique,  described  at  the  end  of 
Sec.  3.4.1,  introduces  an  error  in  the  mean  value  of  the  quantity  measured. 
The  errors  discussed  in  this  section  refer  only  to  those  caused  by  the 
averaging  technique  itself,  and  do  not  include  the  experimental  errors 
covered  in  Appendix  B.  The  method-related  inaccuracies  are  a  function  of: 

(1)  The  nature  of  RC  integration  networks; 

(2)  The  integration  time,  T^nt; 

(3)  The  shape  of  the  average  magnitude  vs.  time  curve  for 
the  quantity  being  measured. 


3.4. 3.1  RC-Integration  Characteristics 


All  five  analog  channels  of  experimental  data  were  averaged  via 
standard  RC-network  integration.  The  solution  to  the  differential  equation 
describing  the  response  of  a  simple  RC-circuit  to  an  impressed  voltage, 
Vj(t),  may  be  written: 


-t/r 

VQ(t)  =  V.(t)-e 


(3.3) 


where 

t  *  RC  (integration  time  constant), 

VQ  *  Output  Voltage  (ideally,  the  average  value  of  V^(t)  at  the  time,  t) 
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Generally,  for  this  experiment,  where  transients  are  begun  from 
initial  steady-state  conditions,  V- (0)  =  VQ(0).  Then  given  any  V.(t), 
Eq.  3.3  defines  the  output  of  the  RC -circuit  as  a  function  of  t  and  t 
Comparison  of  this  VQ( t)  with  the  known  ensemble-average  of  V . ( t )  spe¬ 
cifies  the  error  caused  by  the  special  averaging  technique  employed  in 
this  experiment. 

In  order  to  use  Eq.  3.3,  an  analytic  expression  for  V^(t)  is 
required.  The  stochastic  nature  of  turbulent  velocity  fluctuations  does 
not  permit  straightforward  application  of  Eq.  3.3  to  real  turbulence. 
However,  the  fact  that  any  function  can  be  expressed  as  a  Fourier  series 
suggests  that  a  sinusoidal  approximation  can  be  used  to  bound  the  errors 
involved  and  determine  the  relationship  of  input  signal  frequency,  w, 
and  functional -shape  to  the  error. 

For  example,  consider  the  idealized  case  of  a  mild  coast-down 
transient,  where  the  mean  velocity  and  the  amplitude  of  the  velocity 
fluctuations  are  assumed  to  decay  linearly.  In  such  a  case,  the  instan¬ 
taneous  voltage,  proportional  to  instantaneous  velocity,  may  be 
approximated  by: 

Vi(t)  =  A-Bt  +  C(l-Dt)  sin  cut,  (3.4) 

where,  generally  (A-Bt)  >  10  [C(l-Dt)]  for  all  t. 

Eq.  3.3  gives,  as  the  output  of  an  RC  averaging  circuit  receiving 
Eq.  3.4  as  input, 

V0(t)  =  A-Bt  Correct  ensemble  average 

c  Error  due  to  a  circuit 

+ - [sin  wt-utcosut]  acting  as  a  low  pass 

Hu  t  filter 
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+  — a —  [{l-u^-r^siniot  -  ZiDTCoswt] 

(1-Ho  t)2 

Error  due  to  imposed 
transient 

+  Bt  — —  2  ^  fsintot  -  ajxcoswt] 

1+U)  T 


Exponential  Lag  Term 


(3.5) 


The  error  contribution  to  the  TKE  measurements  is  somewhat  different; 
in  this  case  the  mean  square  amplitude  of  the  fluctuations  is  of  interest. 

A  simplified  representation  during  a  linear  transient  is: 

V^t)  *  (a-8t)  sin2oit.  (3.6) 


Application  of  Eq.  3.3  gives: 


yo  a  \  (o-  st) 


Correct  ensemble  average. 


Error  due  to  circuit 
acting  as  a  low  pass 
filter. 
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imposed  transient 
amplitude. 


(3.7) 


Both  of  the  simplified  input  signals  produce  a  relatively  complicated 
output  which  consists  of  the  "correct"*  mean  value  plus  some  error  terms. 
Fortunately,  most  of  the  error  terms  are  negligibly  small  for  the  frequencies 
and  averaging  time-constants  of  interest;  this  will  be  demonstrated  in  what 
follows.  Moreover,  those  errors  which  are  not  insignificant  are  functions 
of  the  experimental  parameters  and  can  be  minimized  and  bounded  by  proper 
experimental  design.  In  fact,  these  particular  error  effects  ultimately 
determine  the  type  and  rate  of  the  transient  that  can  be  measured  via  the 
pseudo- ensemble  averaging  method. 

*Correct  in  the  sense  that  an  ensemble  average  of  these  input  signals, 
where  phase  angles  are  random,  gives  the  same  mean  function  of  time. 
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3. 4. 3. 2  Optimal  Integration  Time  Constant  (t) 


Inspection  of  Eqs.  3.5  and  3.7  shows  that  selection  of  the  best 
value  of  t  is  a  problem  in  optimization.  A  small  value  oft  minimizes  the 
transient-related  error,  while  a  large  value  of  t  is  needed  to  reduce  the 
frequency-related  error.  Since  transient  effects  are  of  primary  interest, 
the  lowest  value  of  t  consistent  with  expected  turbulent  fluctuation  fre¬ 
quencies,  allows  the  fastest  flow  shut-down  rate  and  represents  the 
optimum  value  for  this  experiment. 

The  smallest  useable  value  of  t  is  limited  by  the  tendency  of  the  RC- 
clrcuit  to  act  as  a  low-pass  filter  and  thereby  fail  to  integrate  that  portion 
of  turbulent  energy  contained  In  the  lower  frequency  fluctuations.  The 
cut-off  frequency,  fc,  of  the  RC -network  determines  the  lowest  frequency  of 
a  periodic  signal  that  is  accurately  integrated.  For  good  integration  of 
periodic  signals  (less  than  1%  error).  Ref.  [34]  shows  that  the  lowest 
frequency  in  the  input  signal,  f^,  must  satisfy 


f,  >  100  fc. 


(3.8) 


where  fc  ■  rar  •  In  terms  of  the  circuit  time  constant,  t=  RC,  proper 
Integration  occurs  for  input  signals  that  satisfy 


(3.9) 


Therefore, T  a  10  sec.,  the  value  used  in  this  experiment,  assures  us  of 
properly  integrating  all  turbulence  energy  contained  in  fluctuations  with 
frequency  greater  than  about  1.6  Hz. 

A  related  effect,  where  random  signals  are  being  processed,  is 
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discussed  in  the  operating  manual  for  the  mean-square  voltmeters  [35]. 

The  expression  for  minimum  t  in  this  case  is 

T  -  • — r  »  (3.10) 

48(4-)  e 

o 

where 

8  =  Bandwidth  in  Hz; 
f  =  Center  frequency; 
e  =  Oesired  error. 

For  e  =  .01,  B  =  500,  and  f  =  250,  conservative  values  for  turbulence  in 
water,  Eq.  3.10  gives  t  >_  10  seconds,  thereby  confirming  our  earlier 
estimate. 

An  order  of  magnitude  study  of  those  error  terms  in  Eqs.  3.5  and  3.7 
that  are  frequency  related  reveals  that,  with  t  =  10  and  turbulent  fluctua- 
tions  of  frequency  >_  2  Hz,  the  worst-case  error  magnitude  is  of  order  10  . 

This  estimate  was  computed  with  values  of  C,  D  and  e  matched  to  the  experi¬ 
ment.  Further,  as  the  minimum  frequency  is  raised  modestly,  the  error 
magnitude  drops  dramatically. 

In  theory,  given  the  nature  of  the  analog  devices  used  to  perform  the 
time-averaging,  the  pseudo-ensemble  method  will  succeed  when  the  fluctua¬ 
tions  are  Gaussian  [36].  However,  since  turbulence  fluctuations  are  only 
approximately  Gaussian,  a  simple  experiment  was  performed  to  determine  the 
magnitude  of  this  effect,  as  well  as  to  evaluate  the  overall  error  intro¬ 
duced  by  use  of  the  optimized  time  constant.  Accordingly,  the  standard 
turbulence  quantities  were  measured  one  inch  above  the  chimney  inlet  at 
full  steady-state  flow,  and  averaged  in  two  different  ways.  The  experiments 


produced  a  negligible  difference  (-1%)  between  the  true  time-average 
(t  =  100  sec.)  and  the  pseudo-ensembl e-average  (t  =  10  sec.)  values  for  both 
mean  velocities  and  turbulence  intensities. 

3. 4. 3. 3  Transient  Shape  Effects 

For  a  transient  experiment,  the  minimum  useable  value  of  t,  with 
respect  to  frequency  related  errors,  is  also  the  optimum  choice.  This  can 
be  seen  by  noticing  that  those  parts  of  the  error  terms  in  Eqs.  3.5  and  3.7 
that  are  related  to  the  imposed  transient  amplitude  are  increased  as  t 
increases.  Once  t  is  established,  the  only  parameters  that  control  the 
error  magnitude  are  those  that  determine  the  transient  speed  (B,  b,  and  D 
in  the  case  of  Eqs.  3.5  and  3.7). 

The  mechanical  characteristics  of  the  control  valve  produced  a 
transient  shape  (i.e.  plot  of  flowrate  vs.  time)  that  was  roughly  congruent 
to  that  produced  by  most  mean  quantities  measured  in  the  experiment 
(e.g.  Fig.  F.i  ).  The  form  of  this  curve  can  be  closely  approximated  by  a 
third-order  polynomial,  and  this  analytic  expression  used  with  Eq.  3.3  to 
estimate  the  size  of  the  error  term.  Since  t  has  been  constrained  to  be 
large  enough  to  keep  frequency  related  errors  negligible,  one  may  simplify 
the  analysis  by  omitting  the  periodic  part  of  V.(t).  This  is  permissible 
because  of  the  linearity  of  the  differential  equation  that  describes  the 
RC  network.  The  frequency- re la ted  error  term  is  simply  superimposed  on  the 
shape-related  error,  and  therefore  can  be  treated  separately.  In  this  case, 
the  analytic  approximation  to  the  shape-related  part  of  the  input  signal 
is 
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Vi(t)  =  at3  +  bt2  +  1  (3.11) 

where , 

a  *  4.672  x  10'7  ; 

b  *  -9.872  x  10"5  ; 

0  <_  t  <_  140  secs. 

Using  Eq.  3.3, 

3  2 

V  (t)  =  at  +  bt  +1  Correct  ensemble  average. 

-  6a -  x2t  +  t3) 

Error  due  to  transient 
,  shape  and  rate. 

-  2b(xt  -  t‘) 

+  e“t7t(6at3  -  2bx2)  Exponential  lag  term.  (3.12) 

The  shape  related  error  is  a  function  of  time.  For  t  =  10  seconds, 
and  a  and  b  matched  to  the  experimental  data,  the  worst-case  error  is 
about  1%  of  the  steady-state  (t=0)  value.  Fortunately,  the  performance 
of  the  averaging  devices  can  be  markedly  improved  by  a  simple  translation 
of  the  time  variable,  i.e.,  let 

t*  3  t-T  , 

then  Eq.  3.12  becomes 

VQ(t*)  =  a(t*)3+  b(t*)2  +  1 

+  aT2[3(t*+x)+T(6e"(t*+T)/T  >5)] 

+  bx2[l-2e_(t*"T)/T].  (3.13) 


In  essence,  the  fixed  time  lag  introduced  by  the  RC-network  has 
been  compensated  for,  and  only  the  curvature-related  errors  are  left. 

The  resulting  worst-case  error,  associated  with  the  shape  effect  for  the 
standard  transient,  is  now  less  than  1%  of  the  steady-state  initial 
magnitude,  an  acceptable  value. 

The  mathematical  time  translation  is  accomplished  in  the  experiment 
and  carried  over  into  the  calculations  as  follows: 

(1)  t  is  fixed  for  all  five  averaging  devices,  for  all 
runs,  and  at  all  locations. 

(2)  The  recorded  inlet  measurements,  which  have  the  time-lag 
built  in,  are  used  to  define  the  time-dependent  inflow 
values  for  the  code. 

(3)  All  comparisons  between  code  and  experiment  are  made 
in  t*-space. 
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Fig.  3.5  TRANSLATION  METHOD 


Fig.  3.6  ELECTRONIC  LAYOUT 


instantaneous 
output  of  LDA 
for  one 
velocity 
component  during 
a  single 
transient  of 
length  tfin. 


Output  of 
time-averaging 
device.  Note 
decreased 
amplitude  and 
frequency; 
also  time 
lag  =  r. 


Final  output 
is  result  of 
ensemble 
averaging  the 
digitized  and 
recorded 


Table  3.1 


Cell-Scale  Mass  Balances 


★ 

I -coordinate 
vertical 
position 

Area  Wei< 
Vertica' 

Upflow 

jhted  Average 
]  Flow  Rate 

Downflow 

%  Mass  Defect  s 
(Downflow-Upflow)  100 

Upflow 

8 

1.22 

1.21 

-0.5 

10 

1.20 

1.17 

-2.5 

12 

1.16 

1.25 

7.7 

14 

1.14 

1.21 

6.1 

16 

1.12 

1.21 

8.0 

18 

1.07 

1.14 

6.5 

20 

1.08 

.95 

-12 

22 

1.12 

.64 

-43 

*  see  Fiq.  C.l  in  Appendix  C 
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Table  3.2 
Equipment  List 


Unit 

Manufacturer 

Model  No. 

Argon  Laser 

Spectra -Physics 

164 

Laser  Exciter 

Spectra-Physics 

265 

Optical  Unit 

DISA 

55L88 

Photomultipliers  (2) 

DISA 

55L12 

High  Voltage  Supplies  (2) 

Baird -Atomic 

31 2A 

Anode  Current  Meters  (2) 

Simpson 

N/A 

LDA  Control  Unit 

DISA 

55L70 

LDA  Frequency  Shifter 

TSI 

985 

Doppler  Signal  Processor 

DISA 

55L22 

Doppler  Signal  Processor 

TSI 

1090 

Turbulence  Processor 

DISA 

52B25 

Averaging  Digital 

Voltmeter 

TSI 

1076 

Integrator 

In-House 

N/A 

Digital  Voltmeter 

Hewlett-Packard 

3456A 

RMS/MSQ  Voltmeters  (2) 

DISA 

55D35 

Oscilloscope 

Hewlett-Packard 

1 41 B 

Low-Pass  Filter 

In-House 

N/A 

Analog-to-Digital 

Converter 

In-House 

DBC  Mark 

Digital  Tape  Recorder 

Tektronix 

4924 

CHAPTER  IV 


COMPUTER  PROGRAM  AND  MODELS  USED  IN  THE  WORK 

4.0  Introduction 

VARR-II  is  a  two-dimensional  time-dependent  thermal -hydraulics  code 
developed  by  Science  Applications,  Inc.  for  the  Westinghouse  Advanced 
Reactor  Division  to  aid  in  the  design  of  the  LMFBR  outlet  plenum.  This 
code  uses  the  two-equation  (K-0)  turbulence  model  developed  by 
Stuhmiller  [37],  A  modified  version  of  VARR-II  was  used  for  all  computer 
predictions  in  this  work,  and  a  number  of  user-set  parameters  in  the  code, 
primarily  affecting  boundary  treatment,  were  exercised.  These  parameters, 
mesh  size  considerations,  and  the  important  modifications  of  VARR-II  and 
their  rationale  are  discussed  in  this  chapter. 

4.1  VARR-II 

The  details  of  the  VARR-II  structure  and  operational  information  are 
found  in  the  user's  manual  [38].  Chen  [39]  also  provides  a  description  of 
VARR-II,  a  complete  list  of  the  partial  differential  equations  involved, 
and  a  concise  explanation  of  the  solution  algorithm.  Accordingly,  this 
section  discusses  only  the  basic  physical  principles  and  special  numerical 
features  of  the  code. 


1 
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4.1.1  General  Principles 

VARR-II  solves  six  simultaneous  partial  differential  equations: 
one  continuity,  two  momentum,  one  energy,  and  two  turbulence  transport 
equations.  Although  strictly  two-dimensional,  the  code  can  handle,  as  an 
option,  a  cylindrical  axisymnetric  geometry  in  r-z  coordinates.  The 
momentum  equations,  differenced  in  their  conservative  form,  include  the 
effects  of  buoyancy  through  density  variations  in  the  gravity  term  and 
contain  optional  resistance  factors  that  can  be  used  to  model  the  effects 
of  mechanical  structure  within  a  mesh  cell.  The  two  turbulence  model 
equations  are  described  in  Sec.  4.3. 

4.1.2  Donor  Cell  Differencing 

An  interpolated  donor  cell  treatment  is  applied  to  the  advective 
terms  in  all  the  transport  equations.  The  user  may  select,  through  a 
donor  cell  coefficient  (a),  center  differencing,  interpolated  donor  cell, 
or  full  donor  cell  differencing.  Center  differencing  provides  higher 
numerical  accuracy  at  the  expense  of  a  more  restrictive  stability 
criterion.  Full  donor  cell  differencing  gives  the  greatest  numerical 
stability,  but  is  less  accurate  than  center  differencing;  the  numerical 
inaccuracies  manifest  themselves  as  an  artificially  enhanced  viscosity 
that  makes  the  flow  more  diffusive  than  the  center  differencing  technique. 
Interpolated  donor  cell,  as  the  name  implies,  is  a  hybrid  combination  of 
the  two  extremes  just  discussed. 

In  this  work,  for  the  mesh  sizes,  time  steps,  and  flow  rates 
examined,  full  donor  cell  differencing  (a  =  a_  =  1.0)  is  always  used  to 
preserve  numerical  stability.  A  test  case,  using  the  coarse  mesh  described 
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in  the  next  section,  showed  that  interpolated  donor  cell  differencing 
(<*x  ”  az  “  *5)  results  in  essentially  identical  mean  velocity  predictions, 
and  slightly  higher  estimates  of  TKE  than  obtained  with  full  donor  cell. 
However,  the  proportional  increase  in  calculated  TKE  throughout  the  cell 
is  roughly  equal  to  the  artificial  increase  in  the  code-specified  inlet 
TKE  that  results  from  using  a  donor  cell  coefficient  not  equal  to  one. 
Accordingly,  the  use  of  full  donor  cell  differencing  is  not  judged  to 
cause  significant  numerical  inaccuracies.  Attempts  to  use  a  donor  cell 
coefficient  less  than  0.5  resulted  in  numerical  instabilities  for  feasibly- 
sized  time  steps. 

4.1.3  Solution  Procedure 

The  solution  procedure  applied  in  VARR-II  is  the  Simplified  Marker 
and  Cell  (SMAC)  method  developed  by  Amsden  and  Harlow  [401  and  described 
in  detail  in  their  report.  The  basic  procedure  starts  with  an  explicit 
computation  of  an  interim  velocity  field  that,  in  general,  does  not 
satisfy  the  continuity  equation.  Boundary  conditions  are  set,  and  then 
cell  pressures  are  adjusted  to  satisfy  continuity  at  each  mesh  cell. 

The  pressure  adjustments  update  the  entire  velocity  field  so  that  this 
step  must  be  performed  iteratively.  The  pressure  iteration  contains  a 
user-adjustable  relaxation  factor  (8q)  that,  when  selected  properly,  can 
improve  the  convergence  rate.  For  stability  reasons,  8Q  is  usually  bounded 

between  1.2  and  2.0.  In  this  work  the  value  B  =  1.85  was  used  for  all 

0 

computations. 


VARR-II  permits  the  use  of  either  cartesian  or  axisymmetric 
cylindrical  coordinates.  Since  the  experiment  was  performed  in  a 
rectangular  test  cell  (see  Sec.  3.1.3),  the  cartesian  system  is  selected. 
The  calculations  as  well  as  the  measurements  are  performed  on  one  half  of 
the  symmetric  plenum  model,  as  shown  in  Fig.  4.1. 

Two  mesh  sizes  are  used  in  the  computational  experiments:  a  10x12 
"coarse"  mesh  and  a  16x24  "fine"  mesh.  The  coarse  mesh  requires  far  less 
computer  time  and  is  used  for  the  extensive  parametric  analysis  that 
examines  the  effects  of  changing  boundary  and  inlet  conditions,  turbulence 
model  constants,  and  other  user-adjustable  features  in  VARR-II.  The  coarse 
mesh  is  depicted  in  Fig.  4.2.  Note  that  this  mesh  does '.not  represent  the 
the  experimental  test-cell  exactly,  especially  at  the  outlet  region. 
However,  comparison  of  the  computed  results  using  both  mesh  sizes  shows 
that  the  coarse  mesh  is  adequate  for  preliminary  estimates,  and  is 
significantly  less  expensive.  The  16x24  fine  mesh.  Fig.  4.1,  is  used  for 
the  final  parameter  studies  and  transient  runs. 

4.3  Modifications 

4.3.1  Turbulence  Model 

In  Sec.  2.1.2  the  rationale  for  using  the  more  complex  two-equation 
turbulence  models  (cf.  one-equation  or  Mixing  Length  models)  in  predictions 
of  recirculating  flows  is  addressed.  One  should  recall  that  the  two 
quantities  needed  to  determine  the  turbulent  viscosity  are  a  turbulent 


velocity  scale  and  a  turbulent  length  scale.  In  the  two-equation  model, 
the  velocity  scale  is  found  from  a  transport  equation  for  K,  the  turbulent 
kinetic  energy.  Attempts  to  write  a  second  transport  equation  directly  for 
the  length  scale,  i,  have  not  been  fruitful .  Instead,  workers  in  the  field 
have  developed  a  second  transport  equation  for  some  other  turbulence 
quantity  that  can  be  related  either  to  l  or  to  the  turbulent  viscosity 
directly. 

The  approach  adopted  in  VARR-II  (and  discussed  in  Sec.  2.1)  is  due  to 
Stuhmiller,  who  writes  the  second  equation  directly  in  terms  of  o,  the 
turbulent  viscosity,  yielding  the  K-a  model.  His  method  is  unique;  the 
more  widely  used  two-equation  model  is  some  form  of  the  K-e  model,  where 
the  second  equation  is  written  for  c,  the  turbulence  energy  dissipation 
rate,  and  a  is  ultimately  found  as  a  function  of  K  and  e.  It  should  be 
pointed  out  again  that  the  two  partial  differential  equations  that  comprise 
the  two-equation  turbulence  models  are  merely  heuristic  approximations  to 
the  unsolvable  exact  equations  that  truly  describe  the  physics  of  the  flow. 
In  this  sense,  the  "correctness"  of  a  particular  model  is  difficult  to 
define.  The  final  assessment  is  probably  best  made  by  asking  which  model, 
for  a  given  amount  of  computational  effort,  gives  results  that  are  closest 
to  valid  experimental  data. 

The  thermal -hydraul ics  code,  TEACH-T  [41]  ,  uses  the  K-e  model 
developed  by  the  staff  in  the  Mechanical  Engineering  Department  at  the 
Imperial  College  of  Science  and  Technology  in  England.  However,  since 
this  code  cannot  calculate  transient  flows,  the  K-e  model  in  TEACH-T 
is  transformed  into  a  mathematically  equivalent  K-a  model  and  included 
in  VARR-II  as  an  option.  In  addition  to  allowing  direct  comparison  of 


the  results  from  both  models,  this  exercise  gives  some  insight  into  the 
modelling  process  itself. 


4. 3. 1.1  K-e  to  K-a  Transformation 


The  two-equation  K-a  model  used  in  VARR-II  is  represented  by  two 
transport  equations  for  K  and  a: 


VARR-II  K-transDort  Eauation 


3t  ax 


_  *  o_2  flUi  +  13.  j  +  £  JL  I  rCT  i£_  I 

i  4k  L3xj  3xd  k  3xj  |_  **LI 

"  I?  ^"[j1  K  (<CI  ’ 4  a]  K  ' 


Here,  r,  a,  ,  and  a-j  are  the  model  free  parameters. 

The  two  equation  K-e  model  used  in  TEACH-T  is  represented  by  two 
transport  equations  for  K  and  e: 


TEACH-T  K- transport  Equation 


IK.  +  U  1K_  -  a 
at  uj  ax.  ° 

J 


°K  3Xj 
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TEACH-T  e -transport  Equation 


h  +  u  fulfil  + jlIjl  sc 

3t  Uj  ITT  K  |_^7  3X^  J  3Xj  3Xj  ae  3YT 


-  C2  ^  .  (4.4) 


Also,  since  a  as  a  function  of  time  and  space  is  the  quantity  of  primary 
interest,  this  model  requires  the  additional  relation: 

c/ 

a  -  -f-  .  (4.5) 

Here,  a^,  C-j ,  a£,  C2,  and  Cy  are  the  model  free  parameters. 

The  desired  transformation  is  obtained  by  using  Eq.  4.5  to  eliminate 
e  wherever  it  appears  in  Eqs.  4.3  and  4.4.  Equation  4.3  becomes 

-  /  3U,  3U .  \2  „  /  _  \  „2 


2 

,  3 

a  3K  j 

+  *T 

J 

aK 

1 

Equation  4.4  (after  some  involved  rearrangement)  becomes 


’ui  +  1 

2 

,  2a  3  , 

a  3K  \ 

3X .  IXT 

K  3Xi  1 

J 
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J  1  J 
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(4.7) 


Equations  4.6  and  4.7  are  simply  re-ordered  forms  of  Eqs.  4.3,  4.4, 
and  4.5.  In  theory,  with  the  exception  of  errors  introduced  by  the  finite 
difference  approximations,  use  of  Eqs.  4.6  and  4.7  in  a  fluid  dynamics 
code  should  produce  the  same  results  as  use  of  Eqs.  4.3,  4.4,  and  4.5.  In 
practice,  truncation  errors  and  differences  in  the  treatment  of  boundary 
conditions  can  be  expected  to  make  the  comparison  less  exact. 

It  is  interesting  to  see  that  two  independent  modellers,  applying 
physical  reasoning  and  judicious  approximation  along  two  entirely  different 
theoretical  pathways,  have  arrived  at  strikingly  similar  forms  for  their 
respective  models.  The  two  TKE  equations  (4.1  and  4.6)  are  term-by-term 
identical,  differing  only  in  the  values  of  the  model  free  parameters.  Even 
the  more  complex  a-transport  equations  (4.2  and  4.7)  share  identical  terms 
with  the  only  difference  being  an  extra  diffusion-like  term,  marked  with 
an  asterisk,  in  Eq.  4.7. 

Although  the  structural  similarity  of  the  two  models  is  made  apparent 
by  this  juxtaposition,  some  specific  differences  are  also  emphasized,  such 
as  the  enhanced  diffusion  effects  mentioned  above.  Another  example  is 
demonstrated  by  considering  the  decay  of  TKE  in  a  homogeneous  flow.  The 
two  "identical"  K-equations,  using  recomnended  values  of  the  irodel  constants, 
predict  turbulent  energy  decay  rates  that  differ  by  a  factor  of  two.  This 
discrepancy  is  a  reminder  of  the  arbitrariness  of  the  modelling  process  and 
the  variation  in  the  available  data  base  from  which  the  model  constants  are 
obtained. 

The  programming  changes  in  VARR-II  required  to  use  the  K-a  ver¬ 
sion  of  the  K-e  model  are  minimized  by  writing  a  generalized  pair  of  transport 
equations  that  represents  both  the  original  and  new  K-o  models.  Equations 


4.1  and  4.6  are  contained  In  this 


Generalized  K- transport  Equation: 
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Likewise,  Eqs.  4.2  and  4.7  are  contained  In  this 


Generalized  a- transport  Equation: 
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(4.9) 

Equations  4.8  and  4.9  represent  either  the  VARR-II  K-a  model  or  the 
TEACH-T  K-e  model,  depending  on  the  values  of  the  k^ .  The  k..  are  functions 
of  the  original  turbulence  model  constants  and  are  defined  in  Table  4.1. 


4.3.2  Inlet  Conditions 

VARR-II  provides  for  a  user-specified  inlet  flow  that  can  be  stated 
in  terms  of  velocity  or  volumetric  flowrate.  For  time-dependent  calcula¬ 
tions,  a  table  of  inlet  velocities  is  prescribed  and  the  code  performs  a 
linear  interpolation  for  time  values  not  in  the  table.  Although  this 
technique  permits  accurate  modelling  of  the  experimentally  determined 
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Inlet  velocities,  the  important  turbulence  quantities,  K  and  o,  are  not 
readily  controlled.  In  a  previous  work,  Chen  and  Golay  [42]  showed  that 
VARR-II  steady-state  predictions  are  sensitive  to  the  specified  levels  of 
inlet  K  and  a.  Accordingly,  the  code  has  been  modified  to  permit  the 
following  user-specified  inlet  turbulent  conditions: 

(1)  Normal  Treatment  -  nlet  and  are  calculated  by 

the  code  as  a  function  of  the  inlet  velocity,  mesh  size, 
and  TQJET  and  TSJET  (where  TQJET  and  TSJET  are  given, 
recommended  constants  based  on  experiment). 

(2)  Steady-State  Matching  Treatment  -  K,  -j  t  is  user- 
specified  at  each  inlet  cell;  ^  t  is  either  ratioed  to 
Kinlet  a  manner  consistent  with  the  normal  treatment, 
or  specified  as  a  constant  multiple  of  the  normal  inlet 
treatment  value. 

(3)  Transient  Matching  Treatment  -  Kiniet  is  user-specified 
at  each  inlet  cell  as  a  function  of  time  (limited  to  a 
maximum  of  three  inlet  cells);  a. ^  t  is  treated  as  in  (2) 
above. 

4.4  Boundary  Conditions 

At  a  solid  wall,  VARR-II  sets  perpendicular  velocities  to  zero  and 
provides  three  options  for  the  parallel  velocity:  free-slip,  no-slip, 
and  derived  boundary  conditions.  The  first  two  are  self-explanatory; 
the  derived  boundary  condition  analytically  prescribes  a  logarithmic 
velocity  profile  between  the  computational  cell  next  to  the  wall  and  the 


wall  itself.  The  purpose  of  this  boundary  condition  is  to  provide  the 
correct  shear  stress  at  the  wall  when  the  mesh:  size  is  too  large  to 
resolve  properly  the  velocity  distribution  near  the  wall. 
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4.5  S  urinary 

The  boundary  conditions,  special  numerical  features,  and  modifications 
have  been  extensively  exercised  in  this  work  while  constructing  the  steady- 
state  flow  field  from  which  the  computational  transient  is  begun.  Proper 
selection  of  these  parameters  and  options  is  important  because  they  often 
have  a  more  significant  impact  on  calculational  results  than  the  turbulence 
model  itself.  The  outcome  of  these  numerical  experiments  performed  on  the 
steady-state  flow  field  is  covered  in  Chapter  Five. 
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CHAPTER  V 

DISCUSSION  OF  EXPERIMENTAL  RESULTS  AND  COMPUTER  PREDICTIONS 
5.0  Introduction 

In  this  chapter  the  experimental  measurements  and  the  computational 
predictions  are  presented  and  discussed.  Since  the  primary  objective  of 
determining  transient  effects  is  most  conclusively  achieved  when  the 
initial  predicted  and  experimental  steady-state  flow  fields  agree  well, 
understanding  this  facet  of  the  problem  is  a  necessary  first  step.  It  is 
noted  in  Chapter  Three  and  Appendix  B  that  experimental  error  in  a 
previous  work  (Chen  [43]),  incurred  while  measuring  the  steady-state 
turbulence  field  in  the  same  apparatus,  may  be  responsible  for  some  of 
the  discrepancy  between  experiment  and  prediction  noted  in  that  study. 
Therefore,  new  steady-state  measurements  are  performed  in  this  work,  and 
the  data  reported  herein  supersede  the  data  of  Chen.  Moreover,  in  an 
effort  to  obtain  the  best  match  between  the  new  experimental  data  and 
the  computational  prediction  of  the  steady-state  full-flow  condition. from 
which  the  coast-down  transient  begins,  a  substantial  number  of  steady-state 
numerical  experiments  are  performed.  As  a  result,  a  significant  portion 
of  the  computational  effort  in  this  work  is  concerned  with  the  steady- 
state  prediction. 
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5.1  Steady-State  or  Full -Flow  Condition 

5.1.1  Experimental  Data 

The  measured  steady-state  values  of  U,  V,  u'2,  v'2,  k  (TKE),  and 
u'v'  are  listed  in  Appendix  C,  Table  C.l.  Estimated  experimental  errors 
are  typically  on  the  order  of  4%  for  the  mean  velocity,  8%  for  K  and  25% 
for  the  Reynolds  Stress.  Appendix  B  contains  the  details  of  the 
experimental  error  analysis.  All  steady-state  full -flow  conditions 
referred  to  in  this  analysis  are  for  an  inlet  Reynolds  number  of  70,000 
as  calculated  in  Appendix  A. 

The  measured  mean. velocity  vectors  are  plotted  in  Fig.  5.1  over 
one  half  of  the  symmetrical  tes ;  cell.  They  exhibit  the  expected 
recirculating  flow  pattern  over  the  -ain  body  of  the  flow.  In  Sec.  3.1.3, 
a  measured  defect  in  the  mass-flow  balance  occurring  at  the  top  of  the 
cell  is  described;  this  discrepancy  is  readily  observed  in  the  velocity 
plot.  However,  over  the  major  part  of  the  vortex  region,  a  good  mass-flow 
balance  obtains,  which  reflects  the  improved  flow  symmetry  in  the 
modified  test-cell. 

In  addition  to  being  listed  in  Table  C.l,  the  steady-state 
measured  K  and  Reynolds  Stress  are  directly  compared  to  calculated  values 
in  Figs.  5.2  and  5.3.  Since  the  scalar  quantity  K  is  accurately  measured 
and  exhibits  a  distinct  spatial  structure  that  should  be  reproduced  by 
an  accurate  turbulence  model,  it  is  considered  a  useful  measure  of 
turbulence  model  effectiveness.  Consequently,  spatial  distributions  of 
K,  in  addition  to  mean  velocity  and  Reynolds  Stress  comparisons,  are 
analyzed  in  the  following  section. 


5.1.2  Comparison  of  Experiment  with  Calculation  (Standard  K-g  and  K-e 
Models) 


5. 1.2.1.  Mean  Velocity  Plots 

Figures  5.4  and  5.5  are  the  calculated  mean  flow  vector  plots  for 
the  K-cr  and  K-e  models  respectively.  These  predictions  of  the  initial 
steady-state  flow  field  result  from  using  the  VARR-II  code  on  the  16  x  24 
fine  mesh  with  the  standard  K-a  turbulence  model  and  the  transformed  K-e 
model  described  in  Sec.  4.3.1.  In  addition,  the  following  code  variables 
are  used  and  referred  to  in  the  future  as  "standard"  settings: 

(1)  inlet  velocities--matched  to  experiment; 

(2)  inlet  turbulence  quantities  (K^n  and  ct^)— set  by  code; 

(3)  turbulence  model  free  parameters— recommended  values; 

(4)  boundary  conditions— free  slip. 

In  order  to  enhance  the  comparison  of  the  experimental  and  computed 
flow  fields,  the  normal  technique  of  printing  an  averaged  velocity  vector 
at  each  computational  mesh  cell  is  not  used.  Instead,  only  velocities  at 
those  cells  for  which  a  measured  velocity  exists  (approximately  every 
other  one)  are  depicted. 

The  K-cr  and  K-e  models  are  seen  to  produce  basically  identical 
mean  flow  fields.  Comparison  of  these  calculated  fields  with  the  experi¬ 
mental  data  (Fig.  5.1)  shows  generally  good  qualitative  and  quantitative 
agreement.  Discounting  the  top  row  for  3-D  effects  in  the  experiment, 
the  vortex  center  location  is  predicted  well.  The  calculated  vortex  is 
elongated,  as  in  the  experiment,  but  with  a  visibly  milder  vertical 
gradient  of  horizontal  velocity  near  the  top.  This  difference  may  be 
an  important  reason  for  observed  discrepancies  in  the  spatial  distribution 
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of  K  that  are  discussed  in  Sec.  5. 1.2. 2. 

The  calculated  upflow  between  the  chimney  entrance  and  the  top 
boundary  of  the  test  cell  overlies  almost  perfectly  the  experimental 
data,  suggesting  that  good  flow  symmetry  is  obtained  in  the  experiment. 

The  steepest  velocity  gradients  are  found,  in  both  the  experiment  and 
calculation,  on  the  right  side  of  the  vortex  adjacent  to  the  inlet  jet; 
however,  the  experimental  results  show  a  somewhat  sharper  gradient  that 
persists  over  the  top  of  the  vortex  which  is  not  predicted  by  the  code. 

The  gradient  discrepancies  are  noted  because  they  drive  production  of 
the  turbulence  parameters  K  and  a,  and  may  explain  the  differences  observed 
between  the  calculated  and  measured  values  of  these  quantities. 

The  calculated  velocity  field  also  shows  slightly  less  downflow 
than  is  observed  experimentally  on  the  left-most  line  of  vertical 
measurements.  This  is  not  surprising  since  the  free-slip  boundary 
condition  in  the  calculation  permits  a  higher  than  normal  average  fluid 
velocity  in  those  cells  adjacent  to  the  wall,  thereby  causing  a  slightly 
lower  than  normal  velocity  to  be  developed  in  the  next  column  of  interior 
cells. 

Overall,  especially  in  the  main  recirculating  zone  away  from 
the  region  of  three-dimensional  effects  on  the  data,  the  mean  velocity 
field  prediction  is  considered  good.  This  judgement  holds  for  both 
models  and,  significantly,  for  the  standard  code  variables  that  include 
the  somewhat  unphysical  free-slip  boundary  condition  and  the  unmatched 
inlet  turbulence  quantities. 
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5. 1.2. 2  Turbulence  Kinetic  Energy  (K) 

It  Is  pointed  out  in  Appendix  B  that  calculated  and  measured 
values  of  K  do  not  represent  the  exact  same  quantity,  and  this  must  be 
understood  before  the  predicted  and  experimental  values  are  compared. 
Even  though  VARR-II  is  a  two-dimensional  code,  it  is  clear  from  the 
theoretical  development  of  Stuhmiller's  turbulence  model  that  the  scalar 
quantity  K,  calculated  in  the  code,  is  defined  by 


i/  _  u'2  +  y'2  +  w'2 

*calc.  ~  2 

where  u',  v',  and  w*  are  the  fluctuating  components  of  velocity  in  all 
three  coordinate  directions.  The  two-channel  experimental  set-up  used  in 
this  work  permits  measurements  in  only  two  dimensions  so  that 


In  the  code,  consideration  of  the  third  dimension  momentum  equation  for  I? 

is  eliminated  by  assuming  that  all  derivatives  with  respect  to  this 

dimension  are  identically  zero.  This  assumption  represents  genuine 

two-dimensionality,  and  the  resulting  equations  apply  exactly.  Likewise, 

the  experiment  is  designed  so  that  W  and  derivatives  in  the  W  direction  are 

zero  or  as  small  as  possible.  However,  w'2  is  not  zero,  since  turbulence 

is  necessarily  three  dimensional.  Therefore,  one  must  account  for  the 

absence  of  w'2  in  the  reported  values  of  K 

r  meas. 

If  the  turbulence  is  isotropic,  then  =  u"1"^  =  v' 2.  Reference 
to  the  measured  data  in  Table  C.l  show  that  u ' 2  and  v ' 2  are  roughly 


■ 
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equal  throughout  the  test-cell.  Since  mean  velocities  and  gradients  in 
the  third  coordinate  direction  are  small  or  non-existent  on  the  plane  of 
the  two-dimensional  measurements,  one  would  not  expect  w'2  to  differ 
significantly  from  its  isotropic  turbulence  level.  In  that  case, 

i  o  ~  u  '2  +  y 1 2 

w  Z  =  - 1 — 

2 

Adding  this  additional  component  of  K  onto  the  measured  values  yields 


Kcalc.  *  '-5  <W> 


This  approximate  relationship  should  be  kept  in  mind  when  comparing 
measured  and  calculated  values  of  turbulence  kinetic  energy. 

Figure  5.2  compares  measured  and  calculated  K  throughout  the  cell 
for  both  models  using  the  standard  code  parameters.  Significant  differ¬ 
ences  are  observed,  even  between  the  two  calculations.  A  clearer  picture 
is  obtained  by  plotting  K  as  a  function  of  horizontal  position  at  a  number 
of  vertical  planes  in  the  cell.  This  is  done  at  three  representative 
locations  in  Fig.  5.6  for  the  data  and  both  turbulence  model  predictions. 

Throughout  the  main  recirculation  zone,  the  data  consistently 
show  two  characteristic  peaks  of  K  on  either  side  of  the  vortex.  The  peaks 
occur  in  the  regions  where  the  velocity  gradients  are  steepest.  This  is 
as  expected  since  K  can  only  be  produced  by  velocity  gradients.  Further, 
the  magnitudes  of  the  measured  peaks  follow  a  logical  pattern.  The  peak 
on  the  inlet  side  is  seen  to  increase  monotonically  as  one  follows  the 
flow  from  just  above  the  chimney  to  the  top  of  the  vortex.  Remembering 
that  the  inlet  K  level  is  dramatically  lowered  by  the  screens  in  the 


modified  inlet,  one  observes  the  inlet  flow  carrying  a  relatively  low  K 
and  entering  a  region  of  high  K  production.  As  the  level  of  K  is 
increased  by  the  strongly  positive  net  production  in  this  area  of  high 
velocity  shear,  the  fluid  proceeds  up  along  the  inlet  jet  constantly 
increasing  its  K  level  in  the  same  manner.  The  K  built  up  in  this  fashion 
appears  to  be  convected  along  a  streamline  (in  fact,  the  two  characteristic 
peaks  occur  on  exactly  the  same  streamline)  and  contributes  to  the 
relatively  larger  K  peak  observed  on  the  downflow  side  of  the  vortex. 

Here  the  velocity  distribution  shows  a  markedly  reduced  velocity  gradient, 
and  one  would  expect  the  production  of  K  to  fall  off.  The  data  are  con¬ 
sistent:  they  show  the  initially  high  K  level  at  the  top  of  the  vortex 
decreasing  monotonically  along  the  downflow  streamline,  indicating  that 
the  dissipation  rate  of  K  in  this  region  has  risen  to  the  point  where  net 
production  is  negative. 

The  calculated  and  measured  ((-distributions  are  compared  in  Fig.  5.6. 
The  location  of  the  K  peak  on  the  upflow  side  of  the  vortex  is  predicted 
reliably  by  both  models  while  the  magnitude  is  overestimated  somewhat. 

The  K-a  and  K-e  distributions  are  similar,  but  the  K-e  model  is  seen  to 
have  a  higher  overestimate  of  magnitude  and  an  apparent  greater  ability  to 
convect  K  over  the  top  of  the  vortex  to  the  downflow  side.  The  most  striking 
disagreement  between  calculation  and  measurement  is  the  absence  of  a  pre¬ 
dicted  K  peak  on  the  downflow  side  of  the  vortex.  A  perceptible  rise  in  K 
close  to  the  wall,  however,  is  observed  in  the  calculations,  and  in  fact 
decays  along  the  downflow  streamline  in  the  same  fashion  as  the  data. 

Possible  reasons  for  these  differences  are  discussed  in  Sec.  5.1.3  where 
attempts  are  made  to  improve  the  numerical  prediction. 
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Although  it  is  true  that  correct  prediction  of  the  Reynolds  Stress 
field  is  all  that  is  required  to  achieve  an  accurate  calculated  velocity 
field,  the  importance  of  correctly  predicting  K,  as  well  as  a  and  the 
Reynolds  Stress,  should  not  be  overlooked.  Since  the  calculated  magnitude 
and  spatial  distribution  of  K  play  an  Integral  role  in  the  resulting 
cr-field,  accurate  Reynolds  Stress  predictions  without  accurate  K  predictions 
are  fortuitous,  but  not  useful  determinants  of  turbulence  model  performance. 
If  one  hopes  to  develop  confidence  in  the  generality  of  a  particular 
turbulence  model,  and  improved  generality  is  the  rationale  for  these  more 
complex  transport-type  models,  then  the  ability  to  predict  K-distributions 
is  required.  Accordingly,  the  spatial  variation  of  K  across  the  main 
recirculating  zone  is  examined  frequently  in  this  work  as  a  useful 
measure  of  turbulence  model  performance. 

5. 1.2. 3  Reynolds  Stress 

Figure  5.3  is  a  comparison  of  steady-state  measured  and  calculated 
Reynolds  Stress  throughout  the  cell  for  both  models  using  the  standard 
code  parameters.  Physical  interpretation  of  the  Reynolds  Stress  distribu¬ 
tion  is  necessarily  less  direct  than  were  the  K-field  comparisons.  This 
is  so  because  of  the  larger  experimental  error  involved  in  measuring  the 
velocity  cross  correlation,  and  because  the  turbulent  viscosity,  a,  which 
is  transported  in  the  computation  and  reported  directly  by  the  code,  can 
only  be  inferred  from  experiment.  Nevertheless,  since  the  calculated 
and  measured  mean  velocity  fields  (and  hence  mean  velocity  gradients) 
are  essentially  congruent,  Reynolds  Stress  comparison  can  give  useful 
qualitative  Information  regarding  the  real  and  predicted  transport  of 
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the  Important  turbulence  quantity  a. 

Comparison  of  the  K-c  and  K-eftodel  estimates  of  Reynolds  Stress 
indicates  qualitatively  similar  spatial  distributions,  with  generally 
higher  magnitudes  in  the  K-e  case.  Since  calculated  velocity  gradients 
are  approximately  identical  In  both  cases,  one  expects  the  calculated  a 
to  be  the  main  contributor  to  the  differences,  and  this  is  in  fact  the  case 
over  most  of  the  flow  field.  Reference  to  the  generalized  K-transport 
equation  (Eq.  4.8)  used  by  both  models  shows  that,  all  other  factors 
being  equal,  a  higher  value  of  a  will  result  in  a  net  increase  in  the 
production  rate  of  K.  This  is  consistent  with  the  overall  higher  levels 
of  K  obtained  with  the  K-e  model  and  discussed  in  the  previous  section. 
Moreover,  these  findings  compare  well  with  the  results  of  TEACH-T  and 
VARR-II  calculations  performed  previously  on  the  same  geometry  [44]. 

Whether  used  In  its  original  form  in  the  TEACH-T  code  or  in  the 
K-cr-equivalent  form  in  the  modified  VARR-II,  the  K-e  model  consistently 
produces  turbulence  quantities  that  are  higher  than  those  predicted  by 
Stuhmi Tier's  K-a  model. 

Comparison  of  the  K-a  Reynolds  Stress  prediction  with  the  data 
are  generally  poor  over  most  of  the  main  recirculating  region  with  the 
notable  exception  of  the  downflow  side  of  the  vortex.  Here  data  and 
prediction  are  in  very  good  qualitative  as  well  as  quantitative  agreement. 
To  emphasize  the  importance  of  correctly  predicting  a  and  the  Reynolds 
Stress,  it  Is  noted  that,  with  the  exception  of  the  top  part  of  the 
test-cell  where  three-dimensional  flow  effects  are  important,  the  major 
mean  flow  discrepancies  always  correspond  to  points  of  significant 
disagreement  between  measured  and  calculated  Reynolds  Stress. 


Since  the  calculated  inlet  turbulence  quantities  that  feed  directly 
into  the  upflow  side  of  the  vortex  are  not  matched  to  the  experimental 
values,  this  discrepancy  may  cause  the  significant  Reynolds  Stress 
mis-match  observed  in  the  high-shear  region  on  the  upflow  side  of  the 
vortex.  Such  a  possibility  is  examined  in  the  following  section  where 
improvements  in  the  steady-state  prediction  are  sought. 

5.1.3  Parametric  Analysis  to  Improve  Steady-State  Prediction 

5. 1.3.1  Rationale,  Method,  and  Scope 

Before  the  transient  calculation  was  performed,  an  attempt  was 
made  to  improve  the  correspondence  between  the  numerical  prediction  and 
the  measurement  for  the  initial  steady-state  full-flow  condition.  The 
experimental  and  calculated  mean  velocities  are  shown  to  be  in  relatively 
close  agreement  in  Sec.  5. 1.2.1,  so  efforts  are  concentrated  on  improving 
the  K  and  Reynolds  Stress  spatial  distributions  where  significant 
differences  between  prediction  and  measurement  have  been  identified. 

The  user-adjustable  parameters  that  have  the  potential  to  change 
the  calculated  turbulence  field  for  a  given  steady-state  inlet  flow  rate 
are  listed  and  defined  below  as  the  following: 

(1)  a  *  a  —  donor  cell  coefficients; 

A  4 

(2)  boundary  conditions; 

(3)  inlet  turbulence  levels  —  the  values  of  K  and  a  ascribed 
to  the  inlet  flow  cells; 

(4)  turbulence  model  free  parameters: 


a  —  dissipation  coefficient  for  K. 
r  —  relative  diffusivity  coefficient  for  K, 
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ctQ  —  dissipation  coefficient  for  o, 

r-|  --  relative  diffusivity  coefficient  for  a. 

Based  on  physical  interpretation  of  the  data  and  comparison  with  the 
original  calculated  predictions,  variations  are  made  in  all  of  the  above 
parameters  in  a  systematic  way  to  determine  if  the  steady-state  numerical 
prediction  can  be  improved. 

Because  VARR-II  runs  are  expensive  in  computer  time,  the  faster- 
running  10  x  12,  cell  coarse-mesh,  described  and  depicted  in  Sec.  4.2, 
is  used  for  this  parametric  analysis.  The  relative  accuracy  of  the 
coarse  mesh  can  be  evaluated  by  comparing  the  results  of  a  fine  mesh 
calculation  to  a  coarse  mesh  calculation  when  all  other  parameters  are 
held  equal.  Figure  5.7  is  the  mean  velocity  field  for  a  coarse  mesh 
calculation  using  the  K-o  model  and  standard  settings  on  the  code  para¬ 
meters.  The  general  features  of  this  prediction  compare  well  with 
Fig.  5.4,  its  fine  mesh  analog.  More  importantly  for  our  purposes, 
spatial  distributions  of  K  and  a  are  essentially  equivalent  for  the  coarse 
and  the  fine  mesh  prediction  over  the  entire  recirculating  zone. 
Significant  differences  are  observed  at  the  level  of  the  plenum  outlet 
because  the  coarsh  mesh  outlet  is  twice  as  large  as  that  of  the  fine 
mesh.  Therefore,  predicted  turbulence  quantities  in  this  region  are  not 
used  in  the  evaluation. 

Approximately  one  hundred  different  combinations  of  the  above 
listed  variableswere  calculated  on  the  10  x  12  coarse  mesh  and  the  results 
compared  to  the  standard  code  prediction  and  the  data.  Table  5.1  depicts 
the  values  assumed  by  these  variables  during  the  course  of  the  study. 

In  order  to  further  limit  the  cost,  only  the  K-a  model  is  used  in  the 


parametric  analysis. 

The  main  features  of  the  experimental  data  that  the  standard  code 
settings  failed  to  predict,  as  discussed  throughout  Sec.  5.1,  are: 

(1)  magnitudes  of  K  and  Reynolds  Stress,  especially  along 
the  upflow  side  of  the  vortex; 

(2)  spatial  distribution  of  K  and  Reynolds  Stress  across  the 
recirculation  zone  (most  notably  the  absence  of  a  second 
K-peak  on  the  down-flow  side  of  the  vortex); 

(3)  the  decay  rate  of  K  along  the  upflow  side  of  the  vortex. 

New  predictions  are  judged  to  be  "improvements"  over  the  standard  code 
settings  when  they  display  a  better  overall  fit  to  the  data  with  respect 
to  the  above  criteria.  In  addition,  the  basically  good  agreement 
achieved  between  predicted  and  measured  mean  flow  speeds  must  be  main¬ 
tained. 

The  major  findings  of  this  exercise  are  discussed  individually  in 
the  sections  that  follow.  Although  they  are  treated  separately  in  the 
text,  various  combinations  of  the  parameter  changes  are  also  accomplished. 
Where  these  combined  effects  produce  unique  results,  they  are  identified. 

5. 1.3. 2  Donor  Call  Coefficient 

The  effect  of  varying  the  donor  cell  coefficient  is  described 
in  Sec.  4.1.2.  There  it  is  shown  that  full  donor  cell  differencing 
(“x  *  az  *  1)  des^rab1e  terms  of  numerical  stability,  and  that 
little  change  in  the  prediction  is  observed  for  values  of  ax  *  a2  as  low 
as  0.5.  In  all  the  following  analyses,  the  full  donor  cell  treatment  is 
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5. 1.3. 3  Boundary  Conditions 

The  alternative  boundary  conditions  available  to  the  VARR-II 
user  are  described  in  Sec.  4.4.  An  improvement  in  the  K-distribution, 
specifically  the  generation  of  a  second  K-peak  on  the  downflow  side  of 
the  vortex,  is  possible  with  a  change  of  the  boundary  condition.  The 
standard  free-slip  condition,  which  is  generally  adequate  for  turbulent 
flow  away  from  the  walls,  under-estimates  the  wall  shear  stress  and 
causes  a  corresponding  underestimate  in  the  production  of  K  near  the 
wall.  However,  if  additional  K  were  produced  along  the  top  wall  by  a 
change  in  the  boundary  condition,  it  might  be  convected  down  into  the 
region  of  K-deficiency  on  the  left  side  of  the  vortex. 

A  series  of  numerical  experiments  was  conducted  in  which  no-slip, 
as  well  as  partial-sl ip,  wall  boundary  conditions  were  applied.  The 
production  of  K  near  the  wall  was  greatly  enhanced  by  the  no-slip  option 
and  reduced  somewhat  in  the  partial-slip  case.  However,  the  substantial 
levels  of  additional  K  generated  by  the  revised  boundary  condition  remained 
near  the  wall  and  were  not  convected  or  diffused  into  the  interior  region 
to  produce  the  second  K-peak  observed  in  the  data.  The  fact  that  this 
second  K-peak  is  not  created  in  the  computation  by  diffusion  of  K  from 
near  the  wall  is  consistent  with  the  data.  There  it  is  seen  that  K  levels 
near  the  wall  are  lower  than  those  at  the  peak  on  the  down-flow  side  of  the 
vortex,  indicating  that  K  is  actually  diffusing  towards  the  wall.  On  the 
other  hand,  the  failure  of  the  turbulence  model  to  convect  K  produced  at 
the  top  wall  into  the  downflow  side  of  the  vortex  is  not  consistent  with 
the  data.  This  behavior  appears  to  be  the  result  of  higher  than  desired 
K-dissipation  rates  in  the  model,  and  is  addressed  further  in  Sec.  5. 1.3. 5. 
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Finally,  the  mean  velocity  field  was  significantly  distorted  in  the 
no-slip  case,  and  the  overall  magnitudes  of  K  throughout  the  cell  were 
increased  markedly  above  the  data  values. 

The  use  of  the  derived  boundary  condition,  or  wall  function  treatment, 
is  sometimes  recommended  for  these  cases  where  the  mesh  size  is  too  large 
to  describe  adequately  the  velocity  profile  near  the  wall.  But  this 
boundary  condition  exacts  a  significant  penalty  in  computer  run  time  [45] 
and  was  not  attempted.  The  results  of  the  no-slip  case,  however,  suggest 
that  the  derived  boundary  condition  could  not  have  offered  substantial 
improvement  since,  as  noted  above,  it  is  seen  in  the  no-slip  numerical 
experiments,  and  in  the  data,  that  diffusion  of  K  from  enhanced  production 
near  the  wall  cannot  be  the  source  of  the  second  K-peak.  As  a  result,  the 
free-slip  boundary  condition  is  seen  to  give  the  most  accurate  mean  flow 
prediction  for  the  lowest  computer  cost  and  is  apparently  not  responsible 
for  the  observed  maldistribution  of  the  turbulence  quantity,  K.  Therefore, 
the  free  slip  boundary  condition  is  used  for  the  transient  computation. 

5. 1.3. 4  Inlet  Turbulence  Levels 

In  addition  to  the  inlet  mass  flow  rate,  values  of  the  turbulence 
quantities  K  and  a  in  the  inlet  flow  must  be  specified  before  the  numerical 
solution  can  proceed.  With  the  standard  code  parameters,  these  values  are 
established  as  a  function  of  inlet  flow  velocity,  the  donor  cell  coefficient, 
and  two  input  variables,  TQJET  and  TSJET.  For  the  case  where  ax  3  a2  =  1  * 
the  relation  is 

K,„  *  TQJET  (UJ2, 


(5.1) 
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and 

ain  =  TSJET  (U.n)  (&X),  (5.2) 

where 

U-n  *  average  Inlet  velocity, 

5X  =  cell  width. 

TQJET  and  TSJET  are  obtained  from  experiments  in  turbulent  jets,  and 
recommended  values  are  given  in  the  VARR  User's  Manual. 

It  has  been  noted,  Chen  and  Golay  [46],  that  the  values  of  the  inlet 
turbulence  quantities  can  have  a  substantial  effect  on  the  predicted  mean 
flow  velocity  distribution  and  that  attention  to  these  quantities  may  be 
as  important  as  attention  to  the  turbulence  model  itself.  In  this  experi¬ 
ment  the  measured  values  of  inlet  K  are  approximately  two  orders  of 
maanitude  lower  than  the  values  set  by  the  code  for  the  same  inlet  flow 
rate;  this  small  measured  value  is  a  result  of  the  turbulence-suppressing 
effect  of  the  inlet  grids. 

It  is  somewhat  encouraging  to  note  that,  even  with  this  large  disparity 
of  K  at  the  inlet,  the  turbulence  model  brings  the  overall  K-levels  into 
rough  agreement  a  short  distance  away  from  the  Inlet.  However,  in  order 
to  achieve  this  agreement,  the  computed  levels  of  K  are  required  to  decay 
along  the  upflow  side  of  the  vortex  in  contradistinction  to  the  slight,  but 
monotonic  build-up  observed  in  the  data.  Therefore,  a  series  of  numerical 
experiments  was  performed  to  determine  the  effect  of  matching  the  inlet 
turbulence  levels  to  the  experiment. 

The  measured  inlet  K-levels  are  known  well  and,  with  the  code 
modifications  described  in  Sec.  4.3.2,  can  be  directly  applied  at  the 


94 


three  inlet  mesh  cells.  However,  inlet  levels  of  a  are  not  readily 
measured  in  this  region  where  mean  velocity  gradients  are  small  and  the 
nearly  Isotropic  turbulence  yields  a  Reynolds  Stress  that  is  close  to 
zero.  When  the  inlet  K  alone  is  matched  to  experiment,  and  inlet  a  is 
calculated  by  the  code  with  Ea.  5.2,  the  resulting  prediction  is 
severely  distorted  by  the  generation  of  excessively  high  values  of  o 
throughout  the  solution  space;  therefore,  some  method  of  "matching" 
is  required.  The  turbulent  jet  experiments  used  to  obtain  Eqs.  5.1  and 
5.2  indicate  that  the  levels  of  inlet  K  and  a  are  related,  and  matching 
only  to  the  experiment  destroys  this  relationship.  Therefore,  in 
cases  where  a.n  is  not  known,  and  K1>n  is  specified  from  experiment,  an 
estimated  value  for  ain  is  found  from  the  combination  of  Eqs.  5.1  and  5.2 
to  be 


When  this  value  of  ain  is  used,  the  mean  field  velocity  predictions  return 
to  their  normal  acceptable  fit  that  was  observed  in  the  unmatched  case. 
However,  there  is  no  improvement  in  the  distribution  of  turbulence  quantities. 
In  fact,  both  the  upflow  K  decay  rate  and  the  spatial  distribution  of  K 
across  the  cell  are  more  poorly  predicted  than  in  the  unmatched  case. 

This  is  a  surprising  result:  3  somewhat  better  overall  numerical 
prediction  is  obtained  in  the  unmatched  case  where  less  information  about 
the  actual  flow  is  used.  This  anomaly  is  thought  to  be  related  to  the 
numerical  method  and  the  turbulence  model  formulation.  Reference  to 
Eqs.  4.8  and  4.9,  the  generalized  transport  model  equations  for  K  and  a. 
show  their  dependence  on  the  ratio  of  a/K.  When  the  very  low  values  of 
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K.  and  a.  associated  with  this  experiment  are  used  as  matched  inlet 
in  in 

conditions,  large  spatial  gradients  in  both  K  and  a  are  established 
above  the  inlet  cells.  Since  the  full  donor  cell  treatment  computes 
derivatives  to  only  first  order  accuracy,  truncation  errors  in  the  largest 
gradient  approximations  could  be  causing  the  unrealistic  buildup  of  o 
that  is  observed. 

While  it  is  certainly  true  that  the  code  mean  flow  and  turbulence 
field  predictions  are  highly  sensitive  to  the  values  of  the  inlet  turbu¬ 
lence  quantities,  improved  predictions  may  not  always  result  from 
matching  these  quantities  to  known  experimental  values.  As  observed 
above,  the  matching  of  K  alone  produced  extremely  degraded  results. 
Therefore,  it  is  recommended  that  any  matching  of  inlet  parameters 
that  Is  attempted  should  roughly  maintain  the  (o/K)^n  ratio  that  Is 
established  by  Eq.  5.3  unless  specific  experimental  values  of  both  K  and 
a  are  known  with  certainty. 

Numerous  combinations  of  boundary  conditions  and  turbulence  model 
free  parameters  were  tried  with  the  matched  inlet  conditions.  In  all 
cases,  the  new  predictions,  considered  on  the  basis  of  the  evaluation 
criteria  listed  In  Sec.  5. 1.3.1,  are  judged  to  be  inferior  to  results 
obtained  with  the  standard  code  parameters  and  unmatched  inlet  conditions. 
As  a  result,  unmatched  inlet  turbulence  values  are  chosen  for  the 
transient  computation. 

5. 1.3. 5  Turbulence  Model  Free  Parameters 

Once  the  boundary  and  inlet  conditions  are  fixed,  the  remaining 
variables  that  can  be  used  to  Improve  the  quality  of  the  steady-state 
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numerical  prediction  are  the  turbulence  model  free  parameters  themselves. 

The  four  free  parameters  used  in  the  K- a  model  are  listed  and  defined  in 
Sec.  5. 1.3.1.  Guidance  as  to  the  specific  parameter  and  the  direction  of 
variation  that  is  most  likely  to  improve  the  prediction  can  be  sought  from 
the  analysis  of  the  experimental  data.  For  example,  the  experimental  data 
in  Fig.  5.6  show  clearly  that  the  magnitude  of  the  K-peak  on  the  upflow 
side  of  the  vortex  increases  as  one  proceeds  from  lower  to  higher  levels 
in  the  test-cell.  However,  both  numerical  predictions  exhibit  the  opposite 
behavior,  suggesting  the  need  for  an  adjustment  to  the  turbulence  model  to 
reduce  the  dissipation  rate  of  K.  Additional  analyses  of  this  nature 
were  made  based  on  other  differences  between  the  measured  and  calculated 
turbulence  fields,  and  the  suggested  changes  to  the  free  parameters  were 
tested.  A  partial  list  of  the  test  runs  accomplished  and  a  qualitative 
evaluation  of  their  outcome  is  given  in  Table  5.2. 

These  results  show  that  the  turbulence  field  and  mean  velocities  are 
highly  sensitive  to  relatively  small  changes  in  the  turbulence  parameters. 

In  addition,  the  complex  coupling  of  the  model  transport  equations  is 
reflected  in  the  somewhat  unpredictable  response  of  the  turbulence  variables 
to  perturbations  in  the  free  parameters.  For  example,  increasing  a,  the 
K  dissipation  coefficient,  results  in  higher  values  of  K  throughout  the 
cell;  decreasing  a  produces  the  opposite  behavior.  On  the  other  hand, 
increasing  aQ,  the  a  dissipation  coefficient,  decreases  both  a  and  K; 
decreasing  aQ  yields  the  opposite  behavior. 

Moderate  adjustments  to  the  relative  diffusivity  coefficients 
(r  and  r-j)  produce  relatively  small  changes,  in  the  expected  direction, 
in  the  spatial  distribution  of  the  scalars  K  and  a.  In  the  development  of 


the  K-o  turbulence  model,  rand  are  shown  to  be  connected  by  a  fixed 
ratio.  Significant  changes  in  the  two  coefficients  can  be  made  as  long 
as  this  original  ratio  is  maintained.  When  one  or  the  other  coefficient 
is  changed  enough  to  substantially  alter  their  ratio,  numerical 
instabilities  result. 

Although  adjustment  of  the  turbulence  parameters  can  slightly 
improve  the  agreement  between  certain  localized  features  of  the  experi¬ 
mental  data  and  the  computed  prediction,  such  changes  always  result  in 
deterioration  at  other  points  of  comparison.  No  combination  of  changed 
turbulence  parameters  has  been  identified  that  can  produce  a  demonstrably 
better  numerical  prediction  of  the  measured  turbulence  field  than  is 
obtained  using  the  standard  recommended  values.  This  statement  applies 
only  to  the  K-a  turbulence  model  because  funds  were  not  available  to 
exercise  the  K-e  model  in  the  same  way. 

An  order-of -magnitude  analysis  of  the  terms  in  the  K-tran sport 
turbulence  modal  equation  show  that  production  and  dissipation  dominate 
over  diffusion  and  advection  throughout  most  of  the  flow.  Further,  al¬ 
though  the  data  show  a  slight  but  persistent  build-up  of  K  with  height 
above  the  chimney  on  the  upflow  side  of  the  vortex,  the  model  consis¬ 
tently  shows  decay.  Even  when  inlet  K  is  matched  to  the  relatively  low 
experimental  values,  the  model  quickly  builds  up  a  high  value  of  K  and 
then  begins  to  decay  with  height.  This  behavior  may  also  be  related  to 
the  aforementioned  failure  of  the  model  to  convect  K  over  the  top  of  the 
vortex  and  into  the  downflow  region.  Although  it  appears  from  the  above 
that  K-dissipation  is  being  over-estimated  by  the  model,  no  tested  change 
to  the  K-dissipation  coefficient  was  able  to  correct  this  behavior.  This 
suggests  that  a  change  in  the  form  of  the  K-dissipation  term  may  be  needed 


to  improve  the  predictions. 


5.2  Transient  Comparison  and  Analysis 
5.2.1  General 

In  order  to  determine  with  certainty  whether  the  use  of  constant 
turbulence  model  free  parameters,  which  are  derived  with  steady-state 
assumptions  and  obtained  from  steady-state  experiments,  adversely  affects 
the  computational  prediction  of  flows  in  the  presence  of  a  momentum 
transient,  it  is  first  necessary  to  minimize  the  differences  between  the 
Initial  steady-state  prediction  and  the  corresponding  measured  steady 
field.  If  this  is  successfully  accomplished,  differences  observed  between 
computed  and  measured  transient  flow  behavior  can,  in  most  cases,  be 
reliably  ascribed  to  the  transient  effect.  The  two-equation  turbulence 
models  examined  in  this  work  provide  steady-state  mean  velocity  field 
predictions  that  generally  match  well  the  experimental  data  in  regions 
of  the  test-cell  where  two  dimensionality  Is  maintained.  Unfortunately, 
the  analysis  of  Secs.  5.1  and  5.2  have  shown  that  the  calculated  scalar 
turbulence  quantities  exhibit  significant  departures  from  the  measured 
values  at  a  number  of  locations.  Further,  the  numerical  experiments 
indicate  that  no  combination  of  revised  boundary  conditions,  inlet 
turbulence  levels,  or  turbulence  model  constants  can  genuinly  improve 
the  Initial  steady-state  match  between  data  and  computation. 

Concern  for  the  adequacy  of  the  steady-state  match  of  the  computed 
and  measured  turbulence  scalars  K  and  Reynolds  Stress  stems  from  two 
sources.  First,  there  is  the  notion  expressed  earlier  that  the  reliable 


prediction  of  such  accurately  measured  and  spatially  structured  turbu¬ 
lence  quantities  as  K  is  needed  before  genuine  confidence  in  the  gener¬ 
ality  of  the  model  can  be  developed.  Second,  any  transient  effects  that 
are  related  to  the  turbulence  model  constants  should  be  most  directly 
observable  in  those  dependent  variables  whose  calculation  is  most  closely 
linked,  to  the  values  of  the  constants.  The  turbulence  kinetic  energy, 

K,  is  the  only  directly  measured  quantity  that  is  dependent,  in  a  first- 
order  sense,  on  the  turbulence  model  free  parameters.  Here  it  is 
recalled  from  the  parametric  analysis  of  Sec.  5.1.3  that  relatively  small 
changes  in  the  values  of  the  individual  turbulence  model  constants  can 
markedly  affect  the  predicted  K  and  a  spatial  distributions  and  absolute 
magnitudes.  This  fact  suggests  that  the  potential  for  a  significant 
transient  effect,  that  is  associated  with  the  model  constants,  does 
exist.  Moreover,  if  steady-state  distributions  of  K  were  accurately 
predicted  by  the  model,  such  transient  effects  could  be  directly  eval¬ 
uated  by  observation  of  the  time-varying  K-field. 

The  predicted  and  measured  steady-state  fields  can  be  brought  into 
initial  agreement  by  a  straightforward  normalization  at  each  point  in  the 
cell,  and  the  computational  transient  calculated  and  compared  with  the 
data  in  this  fashion.  This  is  the  procedure  used  in  the  analysis  that 
follows.  However,  the  non-linear  nature  of  the  problem  makes  it  impos¬ 
sible  to  assure  that  this  initial  normalization  will  compensate  throughout 
the  entire  transient  for  the  apparently  fundamental  inability  of  the  model 
to  make  accurate  steady-state  predictions.  Therefore,  one  cannot  be 
certain  that  differences  observed  in  comparing  the  time-varying  value 
of  a  given  measured  and  calculated  quantity  are  the  sole  result  of 


transient  effects  associated  with  the  use  of  the  constant  turbulence 
model  parameters.  This  fact  is  noted  here  as  a  limitation  on  the 
analysis  of  the  data  that  follow. 


5.2.2  Description  of  Computational  Transient 


The  transient  calculation  was  performed  by  first  establishing  an 
initial  steady-state  flow  pattern  in  the  solution  space  using  measured 
values  of  the  inlet  flow  speeds  as  inlet  velocities  in  the  code.  Once 
the  steady  field  was  developed,  the  inlet  velocities  used  by  the  code 
were  matched  to  the  measured  inlet  velocities  as  they  decreased  during  the 
approximately  120  second  coastdown  transient.  The  general  shape  of  the 
coastdown  as  a  function  of  time  is  similar  to  that  shown  In  the  U  plot  in 
Appendix  F  for  the  position  I  ®  10,  J  8  2,  just  above  the  chimney.  The 
flowrate-versus-time  table  used  by  the  code  was  generated  from  the 
experimental  data  and  is  reproduced  in  Table  5.3. 

In  accordance  with  the  findings  of  the  previous  section,  the 
transient  calculation  is  performed  with  VARR-II  variables  set  as  follows: 

-  16  x  24  fine  mesh 

-  free  slip  boundary  conditions 

-  inlet  turbulence  parameters  set  by  code  (i.e.  not  matched  to 
experiment) 

-  K-o  turbulence  model  with  standard  free  parameter  values. 

To  save  funds,  the  calculation  was  stopped  at  eighty  seconds  into  the 
transient  since  all  observable  trends  in  the  calculation  were  firmly 
established  at  this  point  in  problem-time.  This  observation  is  readily 
verified  by  considering  the  plots  of  turbulence  field  parameters  versus 
time  in  Appendix  F. 


During  the  computed  transient,  all  of  the  turbulence  field  parameters 


were  printed  after  each  ten  seconds  of  problem-time.  This  interval 
corresponds  to  points  in  the  experiment  where  data  were  digitized  and 
recorded  for  the  ensemble  average.  Comparison  of  the  computed  and 
measured  quantities  during  the  transient  Is  presented  In  the  following 
sections. 

5.2.3  Mean  Velocity  Plots 

The  comparison  of  the  calculated  and  the  measured  mean  velocity 
fields  is  presented  In  Figs.  5.8  through  5.13.  Because  differences 
between  calculation  and  measurement  are  small  and  difficult  to  see,  only 
three  sets  of  mean  velocity  plots  are  included.  The  plots  correspond  to 
times  that  span  the  length  of  the  calculated  transient.  Note  also  that 
transient  measurements  are  made  at  only  42  of  the  79  1-0  coordinate 
positions  that  are  used  to  describe  the  steady-state  flow  pattern,  so  the 
vector  plots  in  these  cases  are  less  full. 

The  steady-state  mean  velocity  field  was  shown  to  be  predicted  well 
in  Sec.  5.1,  and  this  observation  generally  holds  for  the  transient  case 
also.  However,  when  the  transient  plots  are  carefully  overlaid  in  sequence, 
some  differences  in  the  computed  and  measured  time  behavior  can  be  discerned. 
The  effect  Is  slight,  but  a  definite  trend  is  observed  wherein  the  calculated 
result  moves  more  than  the  measured  amount  of  flow  to  the  outside  of  the 
vortex  as  the  transient  proceeds.  A  corresponding  deficit  In  down-flow 
velocity  Is  observed  to  be  growing  with  time  near  the  center  of  the  vortex. 
This  behavior  is  related  to  observed  trends  In  the  difference  between  the 
calculated  scalar  turbulence  quantities  and  their  measured  counterparts  in 
the  next  section. 
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5.2.4  Turbulence  Field  Parameters  Versus  Time 

Plots  of  the  turbulence  field  parameters,  U(t),  V(t),  K(t),  and 
-p  u'v'ft),  as  functions  of  time  during  the  transient,  are  provided  in 
Appendix  F  for  each  of  the  I-J  coordinate  locations  where  transient 
measurements  were  made.  On  a  single  page,  for  a  specific  I-J  coordinate 
position,  four  plots  compare  the  data  (+'s)  and  corresponding  calculated 
values  (solid  lines)  throughout  the  transient.  At  those  measurement  points 
where  a  specific  averaged  quantity  is  very  small  and  the  statistics  are 
poor,  a  large  normalization  factor  is  used  to  keep  the  data  scatter  from 
obscuring  adjacent  charts.  This  accounts  for  the  occasional  plot  in 
which  all  the  data  seem  to  fall  on  the  abscissa. 

The  plots  in  Appendix  F  show  that  in  most  cases,  the  normalized 
data  and  calculations  appear  to  follow  the  shape  of  the  general  inlet 
flow-versus-time  curve  defined  by  the  inlet  flow  rate  function  of 
Table  5.3  (the  scalar  quantities  seem  to  vary  as  the  square  of  the  flow 
rate  fractions  in  Table  5.3).  However,  a  close  inspection  of  all  the 
curves  reveals  significant  departures  by  the  data  and  the  calculation 
from  the  general  shape  of  the  inlet  flow  rate  curve.  The  magnitude  and 
sign  of  these  departures  is  seen  to  vary  considerably  throughout  the  flow 
region,  and  this  spatial  "scatter"  of  transient-related  differences 
virtually  eliminates  the  possibility  that  a  systematic  error  in  the 
measurement  technique  is  responsible  for  the  observed  departures. 

The  plots  in  Appendix  F  generally  show,  in  cases  where  the  data 
statistics  are  good,  that  the  discrepancies  between  measured  and 
calculated  transient  behavior  are  not  large.  For  the  flow  speed,  geometry, 
and  transient  rate  used  in  this  experiment,  the  existence  of  a  momentum 
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transient  does  not  cause  differences  between  the  time-dependent  measured 
and  calculated  normalized  turbulence  field  quantities  that  are  large 
compared  to  those  observed  Initially  between  the  two  unnormalized 
steady-state  fields. 

A  qualitative  connection  between  the  time-varying  mean  velocity 
differences  discussed  in  the  previous  section  and  the  plots  in  Appendix  F 
may  be  drawn  by  considering  the  spatial  distribution  of  the  time-dependent 
difference  between  measurement  and  calculation.  In  general,  when  the 
calculated  and  measured  values  of  a  particular  parameter  do  not  coincide 
during  the  transient,  they  tend  to  diverge  from  the  point  of  initial 
steady-state  agreement,  with  one  or  the  other  taking  on  higher' values 
through  time.  The  value  which  stays  high  in  these  cases  is  said  to  "lead"; 
that  which  takes  on  lower  values  is  said  to  "lag."  With  this  definition, 
the  spatial  distributions  of  lead/lag  characteristics  are  plotted  in 
Figs.  5.14  through  5.17. 

In  Figs.  5.14  and  5.15  the  previously-mentioned  tendency  to  over¬ 
estimate  the  amount  of  flow  along  the  outside  wall  during  the  transient 
is  seen  again.  Near  the  center  of  the  vortex  on  the  downflow  side, 
predicted  downflow  (U)  decreases  more  rapidly  than  the  data,  and  predicted 
crossflow  (7)  towards  the  wall  decreases  less  rapidly  than  the  data.  The 
net  result  of  these  two  effects  is  a  tendency  to  show  more  than  the 
measured  amount  of  fluid  along  the  outside  of  the  vortex,  which  trend 
continues  and  grows  as  the  transient  proceeds.  Calculated  upflow  along 
the  outside  of  the  chimney  is  also  seen  to  decrease  less  rapidly  during 
the  transient  than  the  data,  which  corresponds  again  to  the  notion  of 
predicting  more  than  measured  flow  at  the  outside  of  the  vortex. 


These  subtle  departures  from  the  data  of  the  behavior  of  the 
calculated  mean  velocities  is  likely  due  to  the  transient  Reynolds 
Stress  differences  mapped  in  Fig.  5.17.  Though  the  experimental  errors 
are  largest  for  the  Reynolds  Stress  measurements,  a  definite  pattern  of 
the  prediction  decreasing  in-time  less  rapidly  than  the  data  is  developed 
at  most  locations  In  the  vortex  center  and  near  the  cell  exit.  Differ¬ 
ences  in  the  transient  behavior  of  the  mean  flow  occurring  in  the  vortex 
center  are  more  difficult  to  detect  because  of  the  poor  statistics  of  the 
data  when  average  magnitudes  are  near  zero.  However,  when  the  velocity 
field  is  distorted  at  such  points,  the  effect  is  transmitted  throughout 
the  cell  and  observed  more  easily  in  the  regions  where  average  velocities 
are  higher. 

Anomalies  in  the  K-versus-time  comparison  are  generally  distributed 
around  the  outside  of  the  region  of  Reynold  Stress  discrepancy.  An  area 
of  strong  disagremeent  is  observed  in  the  downflow  region  just  inside 
the  vertical  wall.  Here  the  predicted  K  shows  a  more  rapid  than  measured 
rate  of  decrease  in  time  that  grows  during  the  transient.  Since  this  region 
receives  fluid  from  that  part  of  the  cell  where  three-dimensional  effects 
have  been  observed,  it  is  possible  that  this  discrepancy  is  caused  by 
these  effects. 

5.3  Transient  Effects 

The  above  comparison  of  data  and  calculation  shows  that  the 
transient  measurements  are  generally  well-predicted  after  the  initial 
steady-state  differences  are  normalized  out.  In  most  cases,  the 
calculated  curves  follow  the  data  exactly.  However,  at  a  number  of 
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locations  (shaded  areas  in  Figs.  5.14  -  5.17)  small  but  measureable 
differences  in  the  transient  behavior  of  specific  turbulence  quantities 
are  observed.  Such  small  differences  become  even  less  significant  when 
the  initial  normalization  factors  applied  to  the  data  and  calculation 
are  compared.  The  normalization  factors  account  for  the  difference 
between  the  measured  and  calculated  steady-state  values  of  each  quantity 
and  differ  in  some  cases  by  as  much  as  an  order  of  magnitude,  and  with 
more  frequency,  by  factors  of  two  to  five.  With  disparities  of  this 
order  between  data  and  the  turbulence  model  prediction  for  steady-state 
conditions,  the  10  -  20  percent  difference  in  the  transient  comparisons, 
as  is  observed  in  the  cases  of  worst  agreement,  cannot  be  reliably 
ascribed  to  transient  effects  alone. 

An  additional  reason  to  suspect  that  the  differences  seen  between 
data  and  calculation  in  the  plots  of  Appendix  F  do  not  indicate  a  strong 
transient  effect  is  their  spatial  irregularity.  If  the  model  closure 
assumptions  are  deficient  in  the  presence  of  a  momentum  transient, 
then  one  would  expect  this  deficiency  to  be  expressed  throughout  the 
entire  flow  area,  albeit  to  varying  degrees.  The  fact  that  most  of  the 
calculated  quantities  are  seen  to  follow  the  data  closely  during  the 
transient  suggests  that  such  a  general  transient  effect  is  absent. 

Further,  it  is  shown  in  the  steady-state  parametric  analysis  (Sec.  5.1.3) 
that  small  variations  in  the  turbulence  model  free  parameters-  have  a  signi¬ 
ficant  effect  on  the  entire  turbulence  field;  this  also  suggests  that  a 
significant  transient  effect  associated  with  the  turbulence  model  closure 
would  be  more  pervasive  than  is  observed  in  this  experiment. 

Chapter  1  questions  the  use  of  turbulence  models,  developed  from 
steady-state  closure  assumption  and  tested  in  essentially  steady 
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experiments,  in  the  analysis  of  transient  flows.  From  the  above  discussion, 
for  the  flow  speeds,  geometry,  and  transient  rate  used  In  this  experiment, 
it  is  concluded  that  the  steady-state  closure  is  adequate  under  the 
influence  of  a  momentum  transient. 


5.3.1  Characteristic  Time  Estimates 


The  question  of  the  applicability  of  the  above  result  to  flow 
transients  with  a  shutdown  rate  different  from  that  used  in  this 
experiment,  in  particular  faster  rates,  can  be  addressed  with  a  time-scale 
analysis.  This  exercise  also  provides  some  physical  insight  to  the  time- 
scale  over  which  each  term  in  the  turbulence  model  acts.  The  technique 
is  illustrated  for  the  turbulent  kinetic  energy  transport  equation 
(Eq.  4.8)  which  is  shown  to  be  shared  by  both  the  K-a  and  K-e  models  in 
Sec.  4. 3. 1.1. 


For  example,  to  estimate  the  characteristic  time  for  diffusion  of 
K,  select  the  diffusion  term  in  Eq.  4.8  and  let 


(5.4) 


where  a  is  assumed  approximately  constant.  Assume  an  exponential  form 
for  K(t) , 


K(t)  =  KQ  exp  [t/ip],  (5.5) 


where  s  characteristic  time  for  diffusion;  then  approximate 

JL  -  _1 
3X  Lc 


(5.6) 
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1  ft.  for  this  experiment.  Substituting 


where  Lc  *  characteristic  length  * 
Eqs  5.6  and  5.5  into  Eq.  5.4  gives 

v£  ■ 


A  similar  procedure  can  be  applied  to  each  term  in  Eq.  4.8  to  produce  the 
characteristic  times  for  advection  (t^),  production  (rp),  and  dissipation 
(t  ).  with  the  additional  approximations  (based  on  the  experimental 
design  and  a  typical  VARR-II  run) 


Uc  5  2  ft/sec 


a  =  10 


-3  fr 
sec 


^  =  10 


Table  5.4,  showing  the  form  and  estimated  value  for  each  of  the  character¬ 
istic  times  discussed  above,  can  be  produced.  Roughly  speaking,  2.3 
characteristic  times  must  pass  before  the  specific  term  of  interest  can 
'  change  by  an  order  of  magnitude;  therefore,  the  t's  are  a  measure  of 

;  relative  response  time. 

It  is  seen  from  Table  5.4  that  all  terms,  with  the  exception  of 
diffusion,  are  much  faster  than  the  characteristic  time  of  the  imposed 
momentum  transient,  Ty,  which  is  on  the  order  of  100  seconds.  The 
diffusion  time,  Tg,  is  of  the  same  order  as?-]-;  however,  the  magnitude 
of  the  diffusion  term  itself  is  generally  small  in  comparison  to  production 
i.  and  diffusion  (Sec.  5. 1.3. 5).  As  a  result,  three  of  the  four  terms  in  the 


K-equation  are  faster,  by  at  least  an  order  of  magnitude,  than  the 
characteristic  time  of  the  transient,  and  if  the  model  closure  assump¬ 
tions  are  valid,  can  be  expected  to  readily  follow  a  transient  of  the 
order  imposed  in  this  experiment.  The  slower  fourth  is  relatively  small 
in  magnitude,  so  if  the  assumed  form  and  turbulence  model  constant  in 
the  diffusion  term  are  not  valid  in  the  transient,  the  effect  is  not 
readily  detected. 

From  this  analysis,  one  can  estimate  the  fastest  transient  rates 
for  which  the  findings  of  this  chapter  are  applicable.  For  flows  where 
diffusion  is  not  important,  transient  rates  slower  by  an  order  of 
magnitude  than  the  longest  of  t.,  t  ,  or  t  should  produce  equivalent 
results. 

It  would  be  interesting  to  stress  the  turbulence  model  by  examin¬ 
ing  transient  effects  in  a  time-scale  regime  on  the  order  of  1  to  5 
seconds.  Unfortunately,  the  experimental  hardware  required  for  such  an 
investigation  (primarily  the  fast  A-to-D  and  a  much  larger  data  storage 
and  analysis  capability)  is  beyond  the  financial  scope  of  this  effort. 
However,  numerical  experiments  regarding  this  phenomena  are  feasible  and 
included  in  the  recommendations. 


Time  =  40  sec. 


ro  3 


Prediction 
Lags  Data 


Prediction 
Leads  Data 
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Table  5.3 

Percent  Inlet  Flow  vs.  Time 


TTME 

(sec.) 

INLET  VELOCITY 

(%  of  Initial  Full-Flow  Value) 

10 

99.4 

20 

97.6 

30 

94.3 

40 

89.5 

50 

83.1 

60 

75.6 

70 

67.7 

80 

60.0 

90 

53.2 

100 

47.1 

110 

41.8 

120 

38.1 

130 

35.6 

140 

34.7 
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Table  5.4  Results  of  Time  Scale  Analysis 


Term: 

Advection 

Production 

Diffusion 

Dissipation 

Characteristic 

L 

/  . 

2 

c 

/  Le 

a 

Time  (  t  ) 

Uc 

20  he) 

c 

r  o 

4  a  K 

Order  of 

Magnitude  for 

.5  sec. 

5  sec. 

670  sec. 

.5  sec. 

This  Experiment 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

Measurements  have  been  made  in  an  improved  experimental  apparatus 
of  the  mean  turbulence  field  parameters  in  a  recirculating  flow  under 
steady-state  and  transient  flow  conditions.  In  Chapter  3  a  unique 
transient  measurement  technique  is  developed,  and  the  errors  associated 
with  its  use  are  quantified  and  found  to  be  a  function  of  the  transient 
rate  and  shape. 

In  Chapter  4,  a  two-equation  K-e  turbulence  model  is  converted 
to  K-ct  form  for  use  in  the  VARR-II  fluid  dynamics  code,  and  it  is 
found  that  the  two  models  are  nearly  identical  in  mathematical  structure. 
Both  models  predict  well  the  qualitative  and  quantitative  mean  velocity 
fields  under  steady-state  and  transient  flow  conditions.  The  K-o 
model  is  found  to  predict  correctly  the  general  magnitudes  of  the  scalar 
turbulence  quantities,  K  and  Reynolds  Stress,  while  these  quantities 
are  consistently  overestimated  by  the  K-e  model.  With  respect  to  the 
spatial  distribution  of  these  scalars,  both  models  fail  to  predict 
the  measured  build-up  of  K  above  the  chimney  and  on  the  downflow  side 
of  the  vortex.  Since  a  local  balance  of  production  and  dissipation 
are  found  to  dominate  in  the  K-transport  equation,  attempts  were  made 
to  correct  the  predicted  spatial  distribution  by  changing  the  turbulence 
model  free  parameters,  in  particular,  the  parameter  associated  with  the 
dissipation  term.  The  fact  that  this  effort  was  unsuccessful  suggests 
that  the  fundamental  form  of  the  K-dissipation  term  needs  to  be  revised 
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before  improvements  to  the  steady-state  predictions  can  be  made. 

The  values  of  the  scalar  turbulence  quantities  at  the  inlet  to 
the  plenum  are  shown  to  influence  strongly  the  numerical  predictions. 
Matching  the  calculation  values  to  the  experimental  measurements  at  the 
inlet  produces  a  poorer  overall  prediction,  which  is  thought  to  result 
from  errors  in  the  numerical  simulation  that  are  amplified  by  the 
large  spatial  gradients  arising  in  the  matched  inlet  case. 

Comparison  of  transient  data  and  calculation  show  some  small 
differences  in  time-behavior  that  are  generally  less  than  10-20  percent 
at  points  of  worst  agreement.  However,  these  differences  are  not 
thought  to  be  attributable  solely  to  transient  effects  ori  the  turbulence 
model.  Transient  effects  which  occur  in  the  computation  as  a  result 
of  deficient  model  closure  assumptions  are  expected  to  manifest 
themselves,  to  first  order,  through  the  turbulence  model  constants. 

The  fact  that  minor  adjustments  to  the  turbulence  model  constants 
are  shown  to  produce  significant  variations  in  computed  results  suggests 
that  transient  effects  should  produce  a  noticeable  perturbation  through¬ 
out  the  solution  space;  no  such  effect  is  seen.  More  importantly,  the 
failure  of  the  turbulence  model  to  predict  accurately  the  magnitudes 
of  the  measured  steady-state  turbulence  scalars  results  in  start-of- 
transient  normalization  factors  that  vary  widely  throughout  the  plenum. 
These  normalization  factors  are  likely  to  be  functions  of  inlet  flow 
rate,  as  well  as  position,  and  therefore  can  cause  differences  in 
the  comparison  of  transient  data  with  calculation  that  are  not  true 
transient  effects.  Finally,  the  spatial  intermittency  of  transient 
discrepancies  observed  between  data  and  calculation  is  evidence  that 
the  pervasive  effects  one  would  expect  to  be  caused  by  the  turbulence 
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model's  failure  to  handle  transients  are  not  present.  Therefore,  the 
data  indicate  that  the  model's  ability  to  follow  a  momentum  transient 
of  the  order  used  in  the  experiment  is  essentially  verified,  and  that 
improvements  of  the  model's  steady-state  predictions  are  required 
before  more  comprehensive  assessments  can  be  made. 

The  time-scale  analysis  in  Chapter  5  indicates  that  the  con¬ 
clusions  arrived  at  regarding  transient  effects  in  this  experiment 
can  be  applied  to  other  transients  when  the  characteristic-times  of 
the  new  situation  meet  the  criteria  of  Sec.  5.3.1. 

6.2  Recommendations 

Two-equation  turbulence  models  are  thought  to  be  the  simplest 
type  that  can  successfully  predict  recirculating  flows.  Although 
this  experiment  shows  that  some  Important  features  of  the  data,  such  as 
the  mean  velocity  fields,  are  correctly  modelled,  and  that  the  general 
magnitudes  of  the  scalar  quantities  are  well -predicted  in  many  cases. 

It  is  clear  from  the  analysis  of  Sec.  5.1  that  further  effort  is 
required  to  improve  the  steady-state  predictions  of  recirculating 
turbulent  flow.  The  results  of  the  numerical  analysis  accomplished 
to  date  suggest  future  work  as  follows:  (1)  Examine  and  test  alternative 
forms  for  the  K-dissipation  term.  The  features  of  the  K-distribution 
that  are  not  predicted  well  by  the  current  model,  and  that  should  be 
used  as  a  guide  in  the  development  of  the  revised  dissipation  ternij 
are  discussed  in  detail  in  Sec.  5. 1.2. 2.  (2)  Use  a  three-dimensional 

code,  ideally  with  the  same  turbulence  model  used  here,  to  account 
for  possible  effects  from  the  top  part  of  the  cell.  (3)  Obtain 
through  numerical  experiments  a  better  understanding  of  the  problems 
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experienced  in  trying  to  match  the  inlet  turbulence  quantities 
to  the  experimental  data.  Use  of  a  finer  mesh  to  improve  the  numerical 
simulation  of  the  steep  gradients  that  occur  in  the  matched  inlet 
case  is  recommended  as  a  first  step. 

The  time-scale  analysis  showed  that  faster  transients  (with 
characteristic  times  on  the  order  of  1-5  seconds),  would  present  a 
more  severe  test  of  the  model  closure  assumptions.  A  fluid-flow 
experiment  similar  to  the  one  conducted  here  but  with  the  faster 
coastdown  rate  would  be  extremely  expensive;  however,  numerical 
experiments  examining  this  aspect  of  the  problem  would  be  worthwhile. 

The  observed  sensitivity  of  the  computational  predictions  to 
small  changes  in  the  turbulence  model  constants  is  somewhat  surprising. 
The  model  constants  are  generally  established  in  an  approximate  fashion 
from  many  data  sources,  and  one  would  not  expect  the  precision 
requirements  of  such  constants  to  be  high.  It  is  possible  that  the 
observed  sensitivity  is  problem-dependent,  and  it  would  be  useful  to 
know  if  the  same  behavior  is  seen  on  problems  of  larger  scale.  Of 
further  interest  along  the  same  lines,  the  sensitivity  to  changes 
in  the  model  free  parameters  of  quantities  that  are  of  primary 
interest  to  designers  (e.g.  temperature  and  mean  velocities)  should 
also  be  investigated. 
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APPENDIX  A 

DATA  REDUCTION 


Since  the  experiment  is  similar  to  that  performed  by  Chen  [  47], 
this  appendix  closely  follows  his  format.  Further,  an  error  in  his  data 
reduction  method  is  noted  and  quantified  in  the  first  section. 


The  theory  of  Laser  Doppler  Anemometry  shows  that,  with  the  reference 
beam  arrangement,  the  measured  Doppler  frequency  is  linearly  related  to 


the  fluid  velocity  by 


2sine 


(A.l) 


where 

U  =  the  component  of  local  flow  velocity  that  is  normal  to 
the  bisector  of  the  beam  intersection  angle,  and  lies 
in  the  plane  formed  by  the  beams; 
fp  3  Doppler  frequency  (or  shift); 

X  3  wavelength  of  laser; 

8  =  half-angle  of  the  beam  intersection. 

Care  must  be  taken  to  account  properly  for  the  index  of  refraction 
(n)  of  the  medium  in  which  the  measurements  are  made.  The  data  of  Chen 
have  been  called  into  question  in  Sec.  3.1.3  because  of  suspected  flow 
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asymmetry;  here  it  is  noted  that  additional  discrepancies  exist  because 
of  an  apparent  error  involving  the  index  of  refraction  of  water.  In 
Eq.  A.l,  X  and  sin  e  are  quantities  whose  values  depend  on  the  index  of 
refraction  of  the  medium  in  which  they  are  observed.  Hence, 


* water 


X 


air 

water 


and 


(sine) 


water 


<s1ne>ai, 

n  water 


(A. 2) 


Eqs.  A.l  and  A. 2  taken  together  show  that  the  index  of  refraction  effects 
cancel  and  need  not  be  considered.  Chen  corrected  the  angle  0  for  the 
index  of  refraction  of  water,  but  did  not  correct  x;  thereby  biasing  his 
reported  velocities  by  a  factor  of  1.33,  the  index  of  refraction  of  water. 
In  the  present  experiment, 

X  *  5145  x  10"  ®  meters  in  air; 

9 vert. 3  5*69°  in  a1r; 
eHor1z.=  5-58°  1n  a1r' 

Hence, 

8 vert,  [lie]-  <2'59>  fD  • 

and 

>Wtz.  [jlc]  *  <2'65>  fD  [MHzI  '  (A-3 


- - - •_  ■  ..... 

'  - 3 

-  **  '  ‘  . -n  fi.M  a .  i  >-■*•  ■  ■*• --  ' ■*-.  '  ---  ' 


Frequency  shifting  is  employed  to  remove  directional  ambiguity.  In 
such  a  case,  when  the  fluid  is  at  rest,  the  output  voltage  from  the  tracker, 
call  it  VQ,  corresponds  to  zero  velocity.  Since  the  tracker  output  voltage 
minus  the  VQ  offset  is  proportional  to  the  Doppler  frequency,  the  fluid 
velocity  is  (using  the  vertical  component  as  an  example) 


where 


U 


*  (2.59)  •  (V-V  )  •  K 


~MHz“ 
1  [volt_ 


(A. 4) 


K-|-  a  constant  that  depends  on  the  tracker  range  setting; 
Y  =  time  averaged  value  of  tracker  output  voltage. 


1 .2  Mean  Square  Amplitude 

The  output  from  the  tracker  is  a  fluctuating  voltage  that  can  be 
decomposed  into  its  mean  and  fluctuating  components  as 

V(t)  =  V  +  v' (t)  [volts]  .  (A. 5) 


The  two  voltage  components  on  the  RHS  of  Eq.  A. 5  are  proportional, 
respectively,  to  the  decomposed  velocity  components: 

U(t)  *  IT  +  u'(t)  .  (A. 6) 


Therefore,  since  the  proportionality  constants  are  the  same,  Eq.  A. 4  shows 
that  the  mean  square  amplitude  of  the  velocity  fluctuations  is 


Here,  v'^  3  j ■  J v'2(t)dt  is  the  output  of  a  mean  square  voltmeter. 


1.3  Reynolds  Stress 


The  velocity  correlation  is  expressed  as 


V2  *  [{2,59(Kl)(vl)]  •  1(2.65) ( Kg) (vg)  1 

*  6.86(K1K2)(v]v2)  ,  ( 

where  subscripts  one  (1)  and  two  (2)  represent  vertical  and  horizontal 
velocity  components  respectively.  The  quantity  v^  v2  is  the  output  of 
the  DISA  Turbulence  Processor,  an  analog  correlator. 

The  Reynolds  Stress,  in  the  two-dimensional  case,  is  defined  by 


R.S.  5  -  pu,'  u' 


1  2  * 


2.  Cell  Reynolds  Number 


The  Reynolds  Number  is  defined  by 


where: 

U  »  average  velocity  in  the  chimney; 

D  3  chimney  width  3  .188  ft. 

Therefore,  at  full  flow  (F  3  2.64  ft/sec.),  with  temperature  3  110°  F, 
Re  3  70,000.  During  a  transient,  this  value  decreases  monotonically  to 
Re  3  22,750  at  the  end  of  the  coast-down. 
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APPENDIX  B 
ERROR  ANALYSIS 

The  purpose  of  this  appendix  is  to  estimate  the  errors  in  the  veloci¬ 
ties,  mean  square  velocities,  turbulent  kinetic  energies,  and  double 
velocity  correlations  from  the  uncertainties  in  the  LDA  measurement  tech¬ 
nique  and  the  inaccuracies  of  the  electronic  instrumentation. 

1 .  Laser  Doppler  Anemometer  (LDA)  Errors 

Before  demonstrating  the  mathematical  procedure  by  which  the  overall 
experimental  errors  are  estimated,  one  must  evaluate  the  source  and 
magnitude  of  the  individual  error  contributions  associated  with  LDA  opera¬ 
tions.  In  principle,  the  output  from  a  perfect  (error  free)  LDA  is  a 
frequency  that  is  proportional,  with  known  proportionality  constant,  to  the 
fluid  velocity  being  measured.  However,  when  measuring  a  steady-state 
laminar  flow,  a  real  LDA  produces  an  output  that  is  not  the  single  constant 
frequency  that  one  expects,  but  rather  it  is  a  Gaussian-shaped  band  of 
frequencies  centered  about  the  "true"  value.  This  anomaly  arises  from 
factors  associated  with  the  complex  light-scattering  process  in  the  measur¬ 
ing  volume,  and  its  magnitude  can  be  calculated  with  reasonable  precision; 
see,  e.g.,  George  and  Lumley  [48], 

The  various  contributions  to  this  error  are  referred  to  as  Doppler, 
or  frequency,  "broadening"  terms.  The  broadening  errors  affect  turbulence 
measurements  also.  The  magnitude  of  the  turbulent  velocity  fluctuations 
is  derived  Trom  the  spectrum  of  frequency  fluctuations  about  the  mean  that 


are  processed  by  the  tracker,  and  contributions  to  that  spectrum  from 
Doppler  broadening  may  incorrectly  be  interpreted  as  local  fluid  velocity 
fluctuations. 

For  a  frequency  tracker,  such  as  used  in  the  current  experiment, 
the  tracked  frequency  is  proportional  to  the  output  voltage,  E(t),  and 
assuming  the  probability-density  function  of  E(t)  is  symmetrical  and 
Gaussian,  the  mean  square  turbulence  level  is  given  by  Mel  ling  and 
Whitelaw  [49]  as 


/e2 1 
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fe2 
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In  this  equation: 

e  3  RMS  value  of  the  voltage  fluctuations; 

E  =  mean  value  of  E(t); 
t  3  true,  or  correct,  value; 
o  3  observed  value; 
g  =  gradient  broadening; 
f  3  finite  transit  time  broadening; 
n  3  electronic  noise  broadening. 

Equation  B.l  shows  that  if  the  various  broadening  contributions  in  the 
brackets  can  be  evaluated,  the  necessary  corrections  can  be  made  in  a 
straightforward  manner. 

1.1  Gradient  Broadening 


If  a  velocity  gradient  exists  over  the  finite  size  of  the  measuring 
volume,  particles  crossing  the  control  volume  at  different  positions  will 
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have  different  velocities.  These  velocity  differences  cannot  be 
distinguished  from  velocity  variations  with  time,  and  have  to  be  considered 
as  a  broadening  contribution.  For  a  frequency  tracker,  this  term  is 
estimated  by 


(B.2) 


where  ax  is  the  dimension  of  the  scattering  volume  perpendicular  to  the 
direction  of  the  velocity  gradient  measurement  [Mel ling  &  Whitelaw].  This 
correction  is  negligible  for  the  measurements  in  this  report. 


Finite  Transit  Time  Broadeninc 


Due  to  the  finite  life-time  of  the  signal  associated  with  each 
scattering  particle  passing  through  the  measuring  volume,  the  Doppler 
spectrum  exhibits  a  finite  width  over  and  above  that  produced  by  the  fluid 
velocity  fluctuations  alone.  The  random  arrival  and  departure  of  particles 
to  and  from  the  scattering  volume  causes  fluctuations  in  phase,  and  hence 
frequency,  of  the  Doppler  signal.  The  phase  fluctuations  are  correlated 
only  over  the  transit  time  of  particles  across  the  measuring  volume,  which 
depends  on  the  particle  velocity  and  a  scattering  volume  dimension.  The 
broadening  contribution  that  arises  from  this  phenomenon  is  not  consist¬ 
ently  reported  throughout  the  literature;  therefore,  the  most  conservative 
expression  is  used  [Melling  and  Whitelaw]: 


Ng  is  the  number  of  fringes  in  the  measuring  volume.  For  this  experiment. 
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Ng  is  47,  which  gives 


9.6  x  10 


-3 


1.3  Electronic  Noise  Broadening 


RMS  voltage  fluctuations  on  the  LDA  output  may  arise  from  the  noise 
characteristics  of  the  frequency  tracker  circuitry.  The  result  is  an 
additional  broadening  term  that  Melling  and  Whitelaw  evaluated  for  a  DISA 
tracker  like  that  used  in  the  current  work.  They  found  a  worst  case 
value  of 


.002, 


which  will  be  used  for  both  trackers  in  this  analysis. 


(B  .4) 


Referring  to  Eq.  8.1,  the  effect  of  the  total  broadening  correction 

fe' 


is  seen  to  vary  with  the  value  of  t,  the  turbulence  intensity.  For  the 
value  of  the  broadening  contributions  estimated  in  this  section,  Eq.  B.l 
can  be  used  to  show  that  the  error  in  the  measured  value  of  the  mean  square 
fluctuation  is  always  less  than  1%  as  long  as  the  turbulence  intensity  is 
greater  than  3%.  Since  the  experimentally  determined  turbulence  intensity 
is  greater  than  3%  at  all  but  two  locations  near  the  test-cell  inlet,  a  1% 
error  is  assigned  to  the  measurements  of  mean  square  fluctuations  as  the 
result  of  Doppler  broadening. 


2.  General  Principle  of  Error  Calculation 

The  current  experimental  set-up  is  similar  to  that  used  by  Chen  [50]; 
therefore,  his  Error  Analysis  format  is  closely  followed. 

Since  the  errors  come  from  many  sources,  the  method  of  compounding 
errors  suggested  by  Wilson  [51]  is  used. 
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Suppose  that  the  final  result,  Y,  is  related  to  the  components, 
X.j>  by  the  relation 

Y  *  F(Xj,X2,...Xn) 

where  F  is  a  known  functional  form.  A  small  variation  in  the  Xi  will 
alter  Y  by  the  amount 


dY*ll1dXl+lf2d)<2+--+lf/Xn 


3.  Calculation  of  the  Error  in  Velocity 

The  velocity  is  calculated  by  the  expression 

•  <6-6> 

where 

f  *  the  Doppler  frequency 

x  *  the  wavelength  of  laser  light 

9  =  the  half  angle  of  the  beam  intersection 

Hence,  there  are  three  quantities,  f,  x,  and  9  which  will  make  a  contribu¬ 
tion  to  the  errors.  The  error  in  x  is  negligibly  small.  The  errors 
associated  with  9  and  f  are  discussed  below. 

3.1  Error  Due  to  the  Uncertainty  in  8 

A  nominal  value  for  the  angle  9  is  provided  by  the  optics  manufacturer 
as  a  function  of  lens  and  beam-separation  distance.  However,  since  the 
sine  of  a  small  angle  is  a  sensitive  function  of  9,  it  is  recommended  that 
9  be  measured  independently  for  the  particular  optics  arrangement  employed. 
In  this  case,  the  intersecting  beams  were  projected  onto  a  flat  surface 
located  approximately  five  meters  from  their  focal  point,  and  the  angle 
9  was  determined  by  trigonometry.  This  procedure  is  estimated  to  result  in 
a  1%  error  in  the  value  of  9. 

In  Eq.  B.6,  differentiate  V  with  respect  to  9: 

=  -f  xcos9/2sin29  (B.7) 


From  Eq.  B.5,  the  error  (normalized  by  the  velocity)  is 
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1  /3V^2/,„\2  _  ,  A0 

^  a9*  (A9)  -  (tIHe 


*  (^f)2  ,  (B.8) 

where  tan  0  «  0.  (In  this  experiment,  9vertical  =  5.69°  and  9horizontal  = 
5.58°.) 

3.2  Error  due  to  the  Uncertainty  in  f 

The  uncertainty  in  the  frequency  comes  from  the  following  sources: 

1 .  the  accuracy  of  the  frequency  tracker 

2.  the  accuracy  of  the  digital  voltmeter 

3.  the  accuracy  of  the  frequency  shifter 

4.  the  accuracy  of  the  A-to-D  converter 

5.  signal  broadening 

6.  transient  data  averaging  method 

Factor  one,  the  overall  accuracy  of  the  tracker  0s  specified  by  the 
manufacturer  is  1%  of  the  full  scale  deflection.  For  this  experiment, 
the  measured  frequency  was  close  to  50%  of  the  full  scale  reading,  hence 
the  accuracy  is  2.0%.  Factor  two,  the  accuracy  of  the  digital  voltmeter 
is  0.1%.  It  is  small  compared  to  other  errors,  and  thus  will  be 
neglected  in  the  calculation.  For  factor  three,  the  manufacturers  do  not 
specify  the  accuracy  of  their  frequency  shifter.  It  involves  both 
frequency  upshifting  and  downshifting  and  errors  might  occur  during 
these  processes.  The  extent  of  this  error  is  difficult  to  quantify; 
however,  measurements  of  the  frequency  shift  observed  with  no  flow 
indicate  that  it  is  generally  less  than  1%,  so  this  value  will  be 
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conservatively  assigned  to  factor  three. 

Factor  four,  the  accuracy  of  the  A-to-D  converter,  refers  only  to. 
the  quantization  scheme  used  since  linearity  errors  were  negligible  over 
the  range  of  interest,  and  voltage  offset  errors  were  subtracted  out  during 
data  reduction.  For  velocity  measurements  there  were  two  possible 
digitization  modes  corresponding  to  the  particular  tracker  range  in  use. 

The  worst-case  mode  will  be  analyzed  to  give  a  conservative  estimate.  In 
this  case,  ten  volts  were  quantized  over  255  levels,  i.e.. 


10  volts 
Z55 


.0392 


volts 

level 


Therefore,  a  sampled  voltage  could  be  assigned  to  a  level  that  is  in 

0392 

error  by,  at  most,  0.0392  volts.  The  accuracy  is  q-  ■>  0.5%. 

Factor  five,  signal  broadening,  is  discussed  in  Sec.  1  of  this 
Appendix.  The  broadening  correction  applies  mainly  to  the  mean  square 
value  of  the  frequency  fluctuations  (i.e.  TKE),  and  has  a  negligible 
influence  on  the  average  value  of  the  Doppler  frequency  itself  (i.e., 

IT  and  7). 

The  accuracy  of  the  transient  data  averaging  method  is  described  in 
detail  in  Sec.  3.4.3.  It  is  shown  there  that  the  experimental  variables 
were  selected  to  keep  each  contribution  to  the  error  below  1%.  Since  the 
errors  are  time  dependent  and  calculated  in  an  approximate  fashion,  a 
conservative  estimate  of  the  accuracy  assumes  that  all  sources  of  error 
add  throughout  the  transient.  The  accuracy  in  this  case  is  2%. 

In  Eq.  B.6  differentiate  V  with  respect  to  f: 


3  V  =  x 
37  ”  2  sin  0  ‘ 


(B  .9) 
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Hence  the  error  Is 

1  a  u  2  «  Af .  2 

z  (|J4  (Af.)2  =  E  (-J-)  (B.10) 

1  V2  3Ti  1  1  T 

where  i  =  all  six  factors  discussed  above. 

3.3  Other  Errors 

In  the  above  sections,  only  errors  associated  with  Eq.  B.6  were  dis¬ 
cussed.  The  other  possible  errors  are  listed  here: 

1.  Position  errors,  i.e.,  the  measuring  points  are  not 
exactly  the  desired  positions.  In  this  experiment 
the  test  section  is  relatively  large  compared  to  the 
measuring  volume,  hence  this  error  is  judged  to  be 
negligibly  small. 

2.  Errors  in  the  measurement  angle,  i.e.,  the  laser  beam 
is  not  perpendicular  to  the  test  section,  so  that  the 
measured  velocity  component  is  not  the  correct  one. 

A  1%  error  in  velocity  is  assigned  to  this  factor. 

3.4  Sample  Calculation  for  the  Error  in  the  Velocity 

The  total  error  for  the  velocity  measurement  can  be  calculated  by 


01) 


2 


combining  Eqs.  B.5,  B.8,  B.10,  and  other  errors  from  the  previous 
sections: 
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=  0.001125 


Hence,  the  error  in  the  velocity  measurement  is  3.4%. 


4.  Calculation  of  Error  in  the  Turbulent  Kinetic  Energy  (TKE) 

The  turbulent  kinetic  energy,  K,  is  calculated  from  the  expression 

K  «  \  (u '  2  +  v  '  2  )  , 

where  u'2  and  v'2  are  the  respective  mean  square  amplitudes  of  the  velocity 
fluctuations. 

Since  u ' 2  and  v'2  are  measured  independently,  a  conservative  estimate 
of  the  error  in  K  is  calculated  by  adding  the  errors  from  its  components,  i.e. 


Iak  I  .. 

Au '  2 

4. 

Av  '2 

1  K| 

v7! 

(B.ll) 


The  mean  square  amplitude  of  the  velocity  fluctations  is  found 
from  the  expression 


u 


'2  =  f'2 


r  X2  1 


- 

4sin  8 


(B.12) 


where  f'  is  the  amolitude  of  the  fluctuating  component  of  the  Doppler 

t 


frequency,  and  f 


,2’r  f  (f 


( t ) 3 ^dt ,  the  mean  square  amplitude.  The 


errors  in  u'2  come  from  errors  in  and  f'2. 


1  l 
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4. 1  Error  due  to  the  Uncertainty  in  e 

Using  Eq.  B.5,  the  error  normalized  by  u'^  is 


A0 

where  from  Sec.  3.2  in  this  Appendix. 

O 


4.2  Error  due  to  the  Uncertainty  in  f*2 


The  uncertainty  comes  from  the  following  sources: 

1 .  the  accuracy  of  the  frequency  tracker 

2.  the  accuracy  of  the  A-to-D  converter 

3.  the  accuracy  of  the  transient  data  averaging  method 

4.  the  accuracy  of  the  mean  square  voltmeter 

5.  signal  broadening 

Factor  one,  two,  and  three  are  addressed  in  Sec.  3  of  this  Appendix, 
the  values  are  2%,  0.5%,  and  2%  respectively.  For  factor  four,  voltmeter 
accuracy,  one  need  only  consider  the  accuracy  of  the  mean  square  output 
voltage  since  frequency-related  errors  are  a  component  of  factor  three. 
The  accuracy  for  factor  four,  then,  is  given  by  the  manufacturer  as  1%. 
Factor  five,  signal  broadening,  is  estimated  in  Sec.  1  to  be  1%. 

Hence,  the  normalized  error  from  uncertainty  in  f ' ^  is 


(B. 14) 
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4.3  Sample  Calculation  for  the  Error  in  f1? 

The  total  error  is  calculated  by  combining  Eqs.  B.13  and  B.14: 


*  4(.01)2  +  (.02)2  +  ( .005) 2  +  ( .02)2  +  (.01)2  +  ( .01 )2 

*  .001425 


Since  the  error  in  v'2  comes  from  the  same  sources. 


Using  these  values  and  Eq.  B.ll,  the  estimated  error  in  the  measured 
value  of  TKE  is 

=  7.6%  . 

■  2  +  ' 2 

Note  that  this  error  applies  to  the  measured  quantity  - - ^ - 

which  is  defined  here  as  K,  the  turbulent  kinetic  energy.  This  definition 
is  necessarily  artificial  since  it  omits  w'2,  the  third  component  of 
fluctuating  velocity,  which  could  not  be  measured  with  the  two-dimensional 
LOA  system  employed.  The  correct  representation  of  TKE  is 
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K  -  \  (u*2  +  v'2  +  w'2). 

The  absence  of  w‘2  is  not  treated  as  experimental  error  per  se;  its 
probable  magnitude  is  by  no  means  negligible,  and  its  effect  on  the 
comparisons  is  discussed  in  Sec.  5  of  the  main  report. 

5.  Calculation  of  the  Error  in  the  Velocity  Correlation 

The  velocity  correlation  is  represented  by  u'v'.  The  error  comes 
from  the  following  four  sources: 

1 .  Errors  in  u'  and  v' 

2.  The  accuracy  of  the  turbulence  processor 

3.  The  accuracy  of  the  A-to-D  converter 

4.  Transient  data  averaging  method 

In  the  first  factor,  since  the  errors  associated  with  u'  and  v’  are  not 
expected  to  be  correlated,  they  will  not  make  a  contribution  to  the  error 
in  the  velocity  correlation  measurement.  Factor  two,  the  accuracy  of  the 
turbulence  processor  which  does  the  multiplication,  is  specified  by  the 
manufacturer  as  2%  +  10  mV.  Since  a  constant  calibration  was  used,  the 
possible  +10  mV  electronic  drift  error  was  eliminated.  Factor  three,  the 
accuracy  of  the  A-to-D  converter,  is  negligibly  small  in  this  case. 

The  ten  second  time  constant,  t,  used  for  all  five  averaging  circuits 
was  least  satisfactory  for  the  velocity  correlation  measurements.  A 
precise  estimate  of  factor  four  is  not  possible  in  this  case  because  of  a 
lack  of  frequency  Information  concerning  this  signal.  An  estimate  of  the 
upper  bound  of  the  error  associated  with  this  value  can  be  obtained  from 
the  standard  deviation  of  the  u 1 v 1  ensembles  at  each  location.  Although 
these  are  generally  on  the  order  of  10%  of  the  mean,  there  are  numerous 
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locations  where  they  exceed  25%,  especially  where  the  mean  value  of 
the  correlation  is  close  to  zero.  For'  lack  of  a  better  measure,  the  total 
error  in  the  velocity  correlation  is  estimated  to  be  25%. 


i 


APPENDIX  C 


TABULATION  OF  DATA 

For  ease  of  interpretation  and  use,  the  data  in  this  appendix  are 
compiled  in  two  sections.  Section  one  data  are  steady-state  full-flow 
values  obtained  at  the  X-Y  (I-J)  coordinate  positions  listed  in  Table  C.l 
and  defined  over  the  test  cell  in  Fig.  C.l.  Section  two  data  are  a  series 
of  tables,  one  for  each  I-J  position  where  a  transient  was  measured, 
depicting  the  average  value  of  the  turbulence  parameters  as  a  function  of 
time  during  the  transient.  Also  included  for  each  steady-state  and 
transient  mean  value  are  the  estimated  standard  deviations  arising  from 
the  sampling  procedure  used  to  record  the  data. 

General  Statistics 

The  mean  values  and  estimated  standard  deviations  of  the  primary 
measured  quantities  were  found  from 


N 


Here,  in  the  steady-state  case,  N  refers  to  the  number  of  samples  taken 
while  the  flowrate  is  full  and  constant;  in  the  transient  case,  N  is  the 


total  number  of  transient  runs  employed  to  form  the  ensembl 


Also,  for  the  indirect  quantity 


the  estimated  standard  deviation  is  found  from 


where  and  S£  are  associated  with  the  two  components  of  K 
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Table  C.l  (cont.) 

Full  Flow  Steady-State  Case 
(Mean  Values) 


u'2  v'2  u'v' 

TKE 

u 

V 

[m/sec]2  [m/sec]2  [m/sec]2 

[m/sec] 

Position  [m/sec] 

[m/sec] 

x  103  x  IQ4  x  104 

x  10-3 

I  J 

12  2 
4 
6 
8 

10 

12 

14 

16 

14  2 

4 
6 
8 
10 
12 
14 
16 

16  2 
4 
6 
8 
10 
12 
14 
16 

18  2 
4 
6 
8 
10 
12 
14 
16 


.779 

.738 

.155 

.090 

.032 

-.251 

-.509 

-.712 

.757 

.730 

.164 

.077 

-.029 

-.265 

-.499 

-.671 

.740 

.729 

.156 

.055 

-.045 

-.297 

-.471 

-.674 

.664 

.670 

.235 

.053 

-.064 

-.304 

-.406 

-.617 


.0 

.035 

.008 

.008 

.025 

.015 

-.021 

-.028 

.016 

.067 

.024 

.028 

.038 

.023 

-.014 

-.022 

.014 

.009 

.029 

.030 

.038 

.031 

-.009 

-.011 

.032 

.153 

.103 

.063 

.058 

.065 

.032 

.025 


■  .344 

2.22 

2.33 

2.26 
2.03 

6.90 

7.24 
5.54 

.302 

1.90 

2.49 
1.69 

1.33 
5.83 

7.81 
4.96 

.333 

1.81 
3.77 

1.62 

1.50 

6.86 

7.63 
4.05 

.320 

1.35 
6.12 
1.74 
1.61 
7.08 

7.35 
5.18 


.381 

1.26 

2.56 

2.29 

1.71 
6.01 
8.76 

5.25 

.319 

1.27 

2.48 
1.59 
1.31 
6.20 
9.86 

5.49 

.280 

1.25 
3.19 

1.26 

1.26 
7.95 

11.3 

5.12 

.219 

.846 

6.67 

1.27 

1.35 

11.0 

12.9 

4.71 


-.196 

4.90 

3.25 

2.64 

3.51 

19.40 

25.4 

16.4 

.0 

4.61 

4.94 

1.65 

2.90 
18.0 
19.3 

15.1 

.0 

4.14 
9.54 

1.25 
2.75 

19.7 
13.9 

10.8 

.0 
1 .78 
28.0 
1.85 
3.30 

14.2 

4.15 
7.41 


.362 

1.74 

2.45 
2.28 
1.87 

6.45 
8.00 
5.39 

.310 

1.59 

2.48 

1.64 
1.32 
6.02 
8.84 
5.23 

.307 

1.53 

3.48 
1.44 
1.38 
7.41 

9.46 
4.58 

.270 

1.10 

6.40 

1.50 

1.48 
9.03 

10.1 

4.95 
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Table  C.1  (cont. ) 

Full  Flow  Steady-State  Case 
(Mean  Values) 


u'2 

•2 
v  c 

u  V* 

TKE 

U 

V 

2 

[m/sec] 

p 

[m/sec] 

2 

[m/sec] 

[m/sec] 

Position 

[m/sec] 

[m/sec] 

x  103 

x  104 

x  104 

x  103 

I  J 

20  2 

.560 

.045 

.346 

.167 

-.203 

.256 

4 

.603 

.204 

.703 

.452 

~0 

.577 

6 

.376 

.268 

6.91 

6.29 

25.6 

6.60 

8 

.072 

.134 

2.31 

3.34 

4.33 

2.82 

10 

-.081 

.095 

2.83 

3.38 

2.25 

3.11 

12 

-.223 

.107 

6.57 

15.8 

.0 

11.2 

14 

-.295 

.060 

6.29 

12.4 

3.75. 

9.34 

16 

-.560 

.051 

2.41 

3.63 

-0 

3.02 

22  2 

.449 

.048 

.459 

.138 

‘0 

.298 

4 

.478 

.233 

.665 

.  .220 

-.281 

.442 

6 

.445 

.442 

1.58 

.847 

-.905 

1.22 

8 

.212 

.393 

3.78 

5.66 

7.24 

4.72 

10 

.030 

.258 

3.14 

6.54 

-0 

4.84 

12 

-.084 

.199 

4.44 

11.1 

'0 

7.77 

14 

-.234 

.180 

4.73 

7.16 

-3.90 

5.94 

16 

-.444 

.144 

1.95 

3.83 

-3.81 

2.89 

Table  C.l 


Full  Flow  Steady-State  Case 
(Standard  Deviations) 


U 

V 

u'2 

2 

[m/sec] 

v*2 

[m/sec]3 

u  V 

[m/sec]3 

TKE 

[m/sec] 

Position 

[m/sec] 

[m/sec] 

x  103 

x  10J 

x  104 

x  10J 

I  J 

5  6 

.003 

.007 

.127 

.303 

1..59 

.232 

8 

.008 

.011 

.181 

.377 

1.59 

.295 

10 

.012 

.012 

.212 

.532 

2.22 

.399 

12 

.016 

.010 

.358 

.564 

4.40 

.473 

14 

.015 

.005 

.591 

.647 

4.65 

.617 

16 

.017 

.004 

.098 

.277 

1.05 

.208 

17 

.019 

.005 

.086 

.240 

1.26 

.180 

6  6 

.007 

.135 

.421 

1.08 

.313 

8 

.008 

.012 

.132 

.404 

1.46 

.300 

10 

.008 

.007 

.366 

.433 

1.78 

.401 

12 

.014 

.006 

.416 

.696 

3.22 

.573 

14 

.012 

.005 

.272 

.355 

.222 

.317 

16 

.015 

.005 

.600 

.257 

3.76 

.439 

17 

.009 

.005 

.625 

.225 

3.41 

.470 

7.5  2 

.002 

.003 

.008 

.006 

.057 

.006 

3 

~G 

.002 

.010 

.009 

.079 

.010 

4 

.002 

.002 

.075 

.028 

.398 

.054 

8  6 

.012 

.006 

.139 

.221 

1.53 

.171 

8 

.012 

.009 

.134 

.393 

2.25 

.293 

10 

.008 

.263 

.473 

2.47 

.371 

12 

.015 

.008 

.328 

.333 

1.84 

.314 

14 

.011 

.007 

.484 

.319 

4.29 

.410 

16 

.006 

.006 

.271 

.244 

1.91 

.258 

10  2 

.002 

.002 

.008 

.006 

.074 

.007 

4 

.003 

.002 

.039 

.019 

.209 

.030 

6 

.009 

.006 

.016 

.025 

1.40 

.211 

8 

.009 

-  - 

.111 

.315 

1.18 

.234 

10 

.005 

.12? 

.282 

1.64 

.217 

12 

.013 

.004 

.145 

.360 

2.47 

.274 

14 

.009 

.005 

.292 

.418 

2.71 

.360 

16 

.005 

.004 

.178 

.167 

1.27 

.173 

Table  C.l  (cont.) 


Full  Flow  Steady-State  Case 
(Standard  Deviations) 


u'2 

v'2 

u'v' 

TKE 

u 

7 

[m/sec] ^ 

2 

[m/sec] 

2 

[m/sec] 

[m/sec] 

Position 

[m/sec] 

[m/sec] 

x  103 

x  103 

x  104 

x  103 

I  J 

12  2 

.002 

.008 

.008 

.057 

.008 

4 

.002 

.002 

.063 

.039 

.259 

.053 

6 

.008 

.005 

.122 

.155 

1.01 

.129 

8 

.009 

.003 

.105 

.253 

.759 

.193 

10 

.008 

.005 

.098 

.197 

.761 

.155 

12 

.011 

.003 

.418 

.399 

2.38 

.408 

14 

.006 

.004 

.310 

.357 

2.79 

.331 

16 

.005 

.005 

.301 

.236 

2.21 

.270 

14  2 

.002 

.001 

.004 

.004 

.004 

4 

.002 

.005 

.055 

.063 

.471 

.058 

6 

.004 

.003 

.183 

.245 

1.27 

.216 

8 

.007 

.006 

.061 

.108 

.960 

.087 

10 

.009 

.003 

.056 

.099 

1.11 

.080 

12 

.009 

.004 

.245 

.290 

2.22 

.250 

14 

.004 

.005 

.275 

.387 

2.04 

.336 

16 

.006 

.006 

.233 

.236 

1.62 

.230 

16  2 

.002 

.001 

.004 

.007 

.006 

4 

.003 

.001 

.064 

.081 

.349 

.073 

6 

.006 

.004 

.187 

.182 

1.27 

.184 

8 

.008 

.004 

.093 

.105 

.552 

.099 

10 

.005 

.003 

.082 

.155 

.450 

.124 

12 

.008 

.006 

.287 

.385 

1.90 

.321 

14 

.005 

.005 

.173 

.367 

2.68 

.287 

16 

.006 

.006 

.172 

.306 

1.76 

.235 

18  2 

.002 

.001 

.007 

.004 

.006 

4 

.0 

.001 

.026 

.041 

.213 

.035 

6 

.008 

.003 

.363 

.363 

2.99 

.360 

8 

.007 

.003 

.099 

.114 

.594 

.107 

10 

.007 

.004 

.113 

.109 

.875 

.105 

12 

.005 

.007 

.221 

.690 

1.81 

.512 

14 

.008 

.007 

.222 

.688 

2.92 

.492 

16 

.005 

.004 

.341 

.340 

1.62 

.341 

164 


Table  C. 1  (cont. ) 

Full  Flow  Steady-State  Case 
(Standard  Deviations) 


'2 

u 

v'2 

u 1  v' 

TKE 

U 

V 

A 

2 

[m/sec] 

2 

[m/sec] 

[m/sec]2 

[m/sec. 

Position 

2 

[m/sec] L 

[m/sec] c 

x  103 

x  103 

x  104 

x  103 

I  J 

20  2 

.002 

.001 

.009 

.004 

.129 

.007 

4 

.002 

.002 

.020 

.020 

— 

.020 

5 

.005 

.003 

.336 

.250 

1.67 

.296 

6 

.006 

.008 

.100 

.279 

1.18 

.206 

10 

-0 

.003 

.140 

.210 

.941 

.178 

12 

.006 

.009 

.222 

.789 

— 

.579 

14 

.007 

.007 

.272 

.646 

2.08 

.495 

16 

.006 

.006 

.127 

.219 

— 

.174 

22  2 

.003 

.004 

.016 

.052 

.012 

4 

.042 

.001 

.027 

.011 

.120 

.021 

6 

.003 

.005 

.060 

.049 

.446 

.054 

8 

.005 

.002 

.189 

.198 

1.34 

.193 

10 

.006 

.005 

.106 

.425 

.294 

12 

.006 

.006 

.192 

.792 

— 

.533 

14 

.009 

.009 

.223 

.566 

1.87 

.430 

16 

.005 

.004 

.082 

.355 

1.29 

.258 
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Appendix  C  Section  2 
Tables  of  Transient  Data 


2**<S'M>  I  2**<8'W>  I  I  3**<S/'W>  I  CS'M)  I  <S/'M)  |C33S> 

3>IX  I  </*,">  |  <3**,A>  |  |  <ft>  |  <n>  |  3WH 

9  *r  *S  »I  SHOIXISOd 
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ABSTRACT 

This  report  discusses  the  realization  of  a  five-channel 
data  acquisition  system.  The  data  on  these  five  channels  are 
written  on  a  Tektronix  cassette  recorder  via  the  Z-80  system. 
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Features: 


SYSTEM  DBC  MARK  II 


(1)  Digitizing  1  to  5  channels  simultaneously 

(2)  Software  selection  of  1  to  5  channels 

(3)  Variable  sampling  time:  0.1  sec  to  25.5  sec 

(4)  Scattering  sampling  time  capability 

(5)  Variable  header  length 

(6)  System  can  be  expanded  very  easily:  PIO,  S10  &  RAM, 

EPROM  can  be  added  with  little  modifications  in  memory  (EPROM) 

(7)  System  is  a  talker.  It  can  easily  be  turned  into  both 
talker  and  listener;  hence  it  can  be  a  mini  4051  at  much 

lower  cost. 

(8)  External  programming  available.  There  are  128  programming 
steps  plus  two  manual  interrupts. 
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1.1  Introduction 

The  Nuclear  Engineering  Department  currently  operates  a 
Tektronix  4051  Graphic  Terminal.  With  the  recent 
acquisition  of  a  Tektronix  4924  Cassette  Drive,  it  is 
desirable  to  have  a  portable  system  that  can  communicate  with  the 
Cassette  Drive  since  the  4051  is  not  easily  moved  around. 

Since  the  main  use  of  the  cassette  drive  is  logging  data, 
the  desired  system  must  have  the  capabilities  of  accepting 
data  from  the  outside  world  and  the  capabilities  of  interfacing 
this  cassette  drive.  After  some  careful  considerations,  a 
microprocessor  is  best  fitted  for  such  a  system.  The  advantages 
of  using  a  microprocessor  to  implement  this  particular  system 
are: 

(1)  It  minimizes  hardwares. 

(2)  It  uses  program  to  control  all  the  activities.  A 

software  change  is  much  easier  and  cheaper  than 

hardware . 

(3)  It  consumes  less  power. 

(4)  Future  expansion  is  very  easy. 

There  are  disadvantages  of  using  microprocessor 
also.  Microprocessor  nowaday  are  still  very  slow 
and  are  quite  expensive.  However  speed  is  not  the 
main  concern  for  this  particular  application. 

1.2  Basic  System  Configuration  (Please  refer  to  the  System  Block 

Diagram) 

The  microprocessor  employed  for  this  application  is  the  Z80 
CPU.  The  system  also  consists  of  one  K  bytes  of  EPROM,  256  bytes 
of  static  RAM,  three  PIO’s  (Parallel  Input  Output  Port),  and  a 
CTC  (Counter  Time  Control) . 
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There  are  five  analog  inputs  which  will  accept  from  one 
to  five  external  voltages.  These  voltages  are  digitized 
sequentially  through  the  A/D  Circuitry  and  the  A/D  Timing 
Control  and  Input  Channel  Multiplexor.  The  digitized  data 
are  transferred  to  the  RAM  via  PIOA.  CTC  provides  the  timings 
for  the  A/D.  PIOB  is  used  for  interfacing  the  Tecktronix  4924 
Cassette  Drive.  PIOC  is  used  to  multiplex  the  10  digit  displays 
which  indicate  current  run  number,  number  of  sample  points  col¬ 
lected  and  the  number  of  sample  points  wanted.  — 

The  EPROM  is  the  brain  of  the  system.  It  contains  all  the 
operational  programs.  The  Z80  CPU  executes  the  instructions  supplied 
by  the  EPROM.  The  Data  Entry  Control  Logics  serves  two  important 
functions.  It  is  a  pseudo-DMA  (direct  memory  access)  control 
to  the  static  RAM.  It  is  also  maintenance  control.  This  is 
an  extremely  helpful  function  since  no  logic  analyzer  is 
available  to  debug  the  system.  There  are  three  buses  used  in 
the  system:  A  16-bit  Address  Bus,  an  8-bit  Data  Bus  and  a 
Control  Bus. 

1.3  System  Functions 

The  system  was  first  designed  to  perform  the  exact  functions 
as  the  Tektronix  4051.  However,  due  to  limited  space  of  the 
EPROM,  most  of  the  4051  statements  were  not  implemented.  The  system  is 
implemented  as  a  Talker  which  can  execute  the  FIND,  MARK  and 
WBYTE  Statements.  There  is  one  PIO  port  open  in  case  the 
system  is  needed  to  be  implemented  as  a  Listener.  Additional 
statements  can  also  be  implemented  when  more  EPROM  is  added. 

The  system  is  designed  under  the  IEEE  488  GPIB  standards. 

The  system  essentially  is  capable  of  executing  the  following 


statements: 
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FIND  @  2 : N 
MARK  @  2 :M, X 

WBYTE  @  34:63,95 

WBYTE  @34, 121: data 


WBYTE  @  121  .-data 


;N=nth  record,  OS  N<  9999 
;M*#  of  records,  0<M*255 
;X=#  of  bytes  per  record,  0<X*9999 
;this  statement  terminates  the 
handshake  operation  on  unit  2 
;this  statement  calls  for  a  byte 
>by  byte  transfer  operation.  Data 
;can  be  transferred  to  the  unit  as 
; soon  as  the  statement  is  executed 
;this  statement  put  the  system  in  a 
continuous  data  transfer  mode 
;without  terminating  and  restarting 
;  the  GPIB . 


The  system  is  designed  to  be  the  Talker  and  any  Listener 
handshake  with  the  system  is  designated  for  unit  2.  This  is  a 
fixed  unit  number.  This  number  can  be  changed  only  by  re¬ 
programming  the  EPROM. 

1.4  System  Implementation 

Since  no  TTY  is  connected  to  the  system,  the  only  way  to 
initiate  operations  are  by  depressing  switches. 

When  power  is  applied,  the  system  is  put  in  a  known  state 
through  the  Initialization  program.  All  the  PIOS  and  CTC  are 
provided  with  the  required  Control  Words  and  Interrupt  Vectors. 
Once  the  system  is  initialized  and  interrupt  mode  2  is  enabled, 
the  system  is  idle  and  loops  in  the  BCDPLAY  subroutine  which  will 
scan  all  ten  digits  sequentially. 

At  this  point,  user  can  put  the  system  in  Local  mode  by 
flipping  the  LOCAL/REMOTE  switch.  The  MEMORY/DATA  will  indicate 
whether  the  content  of  the  memory  is  being  scanned  or  data  is  to 
be  written  into  memory.  When  MEMORY/DATA  is  at  MEMORY  and 
READ/WRITE  is  at  READ,  the  content  of  a  specified  memory  location 
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is  scanned.  The  first  three  7- segment  displays  indicate  the 
memory  location  and  the  remaining  two  displays  indicate  its 
content.  The  address  as  well  as  content  are  hexidecimally 
represented.  When  MEMORY/DATA  is  at  DaPA  and  READ/WRITE  is 
at  WRITE ,  a  memory  write  is  to  be  performed.  The  basic 
procedures  for  doing  memory  read  or  memory  write  are: 

(1)  put  system  in  local 

(2)  key  in  desired  address 

(3)  select  read  or  write  operation. 

There  are  three  other  switches  which  are  associated  with 
this  function.  The  UPCOUNT  increments  the  displayed  address 
by  one.  The  DOWNCOUNT  decrements  the  memory  address  by  one. 

The  MEMORY/DATA  CLR  will  clear  the  Memory  Address  register 
if  it  is  a  memory  read  or  the  data  register  if  it  is  memory 
write.  The  actual  memory  write  occurs  when  the  second  half  of 
the  8-bit  byte  is  entered.  The  memory  address  will  automatically 
be  incremented  by  one  when  the  word  is  written. 

The  A/D  CLR  switch  clears the  analog  to  digital  circuit, 
unpicks  the  Solid  State  relay  and  clears  the  CTC.  The  SYS 
RESTART  switch  reinitializes  the  system.  There  are  five 
program  interrupt  switches.  They  are  WRITEHEADER,  FIND, 
USERINTERRUPT1 ,  START  and  STOP. 

The  WRITEHEADER  switch  interrupts  the  system  and  writes  the 
necessary  header  on  the  cassette  tape.  It  essentially  performs 
the  following  program  in  basic: 

N  =  1 

FIND  @  2:N 

MARK  @  2:1,X  ;X  is  pre-loaded  in  memory 

FIND  @  2 :N 

WBYE  @  34, 121:header  ;length  of  header  varies 

N  =  N  +  1 
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The  WRITEHEADER  switch  also  loads  the  CTC  with  the  correct 
time  constants  such  as  sampling  time  and  number  of  channels. 

Once  the  header  is  written,  the  system  is  ready  for  quantizing 
the  inputs.  The  START  switch  applies  an  external  trigger 
pulse  to  initiate  the  CTC  operation.  This  trigger  pulse  also 
triggers  a  50  microseconds  sample  and  hold  pulse.  All  five 
channels  are  sampled  at  the  same  instance.  The  trailing  edge  of 
and  H  pulse  triggers  a  3  microseconds  pulse.  This  pulse 
in  conjunction  with  the  System  Clock,  produces  the  start 
signal  required  to  start  the  A  to  D.  When  the  conversion  is 
done,  a  few  operations  are  performed  before  the  next  conversion 
can  begin.  The  next  analog  switch  is  turned  on,  supplying  the 
A  to  D  input  with  new  analog  signal  to  be  digitized.  A 
downcount  pulse  is  generated  and  applied  to  the  CTC  to  downcount 
the  number  of  channels  remaining  to  be  digitized.  The  previous 
digitized  signal  is  written  into  the  scratch  pad  memory  at 
location  zero.  The  scratch  pad  memory  is  then  upcounted  by  one. 

When  TO 2  of  CTC  comes  on,  it  indicates  that  the  channel 
being  quantized  is  the  last  channel.  As  soon  as  this  last 
channel  is  written  into  the  scratch  pad  memory,  a  DONE  signal 
is  generated.  This  DONE  signal  then  interrupts  the  system. 

The  system  services  this  interrupt  in  this  manner.  It  first 
checks  if  the  scattering  sample  operation  is  selected.  If  it 
has  been  selected,  then  the  system  checks  for  the  number  of 
STC  left  and  then  loads  the  necessary  sampling  time  for  the 
next  operation.  The  system  then  proceeds  to  read  in  the 
digitized  data  and  store  them  at  memory  location  starting  at 
4DA.  Once  all  the  data  are  written  into  the  RAM.  The  data  are 
transferred  to  the  GPIB  one  byte  at  a  time  until  all  digitized 
channels  are  written  on  the  Cassette  Drive.  A  dummy  byte  is 
also  written  to  provide  a  checking  mechanism  for  the  number  of 
samples  collected.  The  system  then  increments  the  sample 
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points  collected  register  by  one.  This  incremented  number  is 
compared  with  set  number  of  sample  points.  If  they  are  equal, 
the  system  will  stop  by  resetting  the  CTC  time  constant. 

The  STOP  switch  has  the  identical  function  as  mentioned  in 
the  last  sentence.  The  system  then  terminates  the  handshake 
operation  by  sending  the  EOI  (end  of  information) ,  UNTALK  and 
UNLISTEN.  If  the  comparison  is  not  equal,  the  system  will  remain 
in  the  wait  mode  for  more  data.  The  GPIB  is  not  terminated  in 
this  case. 

One  thing  should  be  clarified  here.  In  the  previous 
paragraph,  it  is  mentioned  that  data  are  first  transferred  to  the 
static  RAM,  then  to  the  GPIB.  This  may  seem  to  be  a  very  foolish 
idea.  There  are  two  reasons  for  doing  it  in  such  a  way.  From 
the  maintainence  standpoint,  it  is  easier  to  check  the  memory 
locations  to  see  if  the  A  to  D  are  working  correctly  than  to 
carry  the  Cassette  Drive  to  the  4051  to  verify  the  results. 

The  second  reason  is  that  transferring  a  group  of  data  at  a  time 
saves  more  programming  steps,  hence  more  EPROM  space  available 
for  other  programs.  The  fact  that  the  minimum  sample  time  is 
0.1  second  provides  ample  time  for  time-sharing  the  CPU. 
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Chapter  2 


SYSTEM  OPERATING  PROCEDURES 


COMMENTS : 

To  enter  information,  one  needs  to  do: 


(1)  MEM/DATA  at  MEM 

(2)  READ/WRITE  at  READ 

(3)  Press  MEM/DATA  CLR 

(4)  Key  in  Address 

(5)  MEM/DATA  at  DATA 

(6)  READ/WRITE  at  WRITE 

(7)  Press  MEM/DATACLR 
(81  Key  in  Data 

(9)  MEM/DATA  at  MEM 
(10)  READ/WRITE  READ 


(A)  If  the  tape  is  a  brand  new  one,  perform  Step  1  to  7.  If 
tape  is  a  used  tape,  go  to  Part  B. 

(1)  turn  power  on  both  the  system  and  the  Cassette  drive 

(2)  press  SYSRESTART 

(3)  press  A/DCLR 

(4)  switch  to  "LOCAL" 

(5)  Enter 


Address 

Content 

4F*5 

00 

4F6 

00 

4E6 

00 

4E7 

10 

4E8 

00 

4E9 

01 

403 

BS 

D2 

02 

406 

CD 

5E 

03 

409 

CD 

7C 

03 

40C 

CD 

AF 

02 
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(6)  press  "USERINTERRUPT1”  once.  Observe  that  tape  will 
be  marked  and  then  stopped 

(7)  go  to  Part  B  Step  12. 

(B)  If  tape  is  used,  perform  the  following  steps: 


B.l  For  the  very  first  run: 


(8) 

(9) 

(10) 
(11) 
(12) 


4F5 

00 

Program  FIND 

4F6 

01 

4D1 

#  STC 

4D2 

MSD  TC 

4D3 

LSD  TC 

4D4 

MSD  i TC 

Duration  of  i.TC 

4D5 

LSD  &  TC 

4EO 

#  ch.  in 

use 

Starts  00  to  04? 

4E4 

MSD  #  of 

samples 

4E5 

LSD  wanted 

4E6 

MSD  #  of 

bytes  in 

BYTE  #  in  MARK  N, I 

4E7 

LSD  the  next  run 

4E8 

MSD  00 

Record  #  in  MARK 

4E9 

LSD  01 

N,X 

turn  power  on  both  the  system  and  the  cassette  drive 
press  SYSRESTART 
press  A/D  CLR 
switch  to  "LOCAL" 

enter  Address  Content  Comments 


For  the  location  of  Dummy  data,  it  depends  on  how  many 
channels  are  in  use. 

For  1  channel 


2  channels 


4  DA  =  Data 
4DB  =  Dummy  data 
4DA, 4DB  *  Data 
4 DC  =  Dummy  data 


•t 


3  channels 


4DA, 4DB , 4DC  -  Data 
4DD  =  Dummy  data 

4  channels  4DA,4DB,4DC,4DD  -  Data 

4DE  =  Dummy  data 

5  channels  4DA, 4DB, 4DC, 4DD, 4DE  =  Data 

4DF  =  Dummy  data 

(13)  enter  header  starting  at  4A0  for  a  maximum  of 
48  bytes 

(14)  switch  to  "REMOTE" 

(15)  press  "WRHEADER"  to  start  writing  the  header 

(16)  when  cassette  drive  stops,  press  "START"  to  start 
the  current  run 

(17)  go  to  Part  C  Step  24  for  next  run  operating  procedures 

B.2  If  the  operation  is  not  the  first  run,  the  last  run  #  must 
be  known.  Then  proceed  as  following: 

(18)  turn  power  on  both  the  system  and  the  cassette  drive 

(19)  press  SYSRESTART 

(20)  Press  A/D  CLR 

(21)  Switch  to  "LOCAL" 

(22)  enter  all  required  parameters  as  described  in  Step  12 
except  locations  4F5  and  4F6  and  enter 

4F5  0  last  run  #  plus  one 

4F6  ^ 

4E1  last  run  # 

(23)  go  to  Step  13. 

(C)  For  the  next  run  operation,  perform  the  following: 

(24)  DO  NOT  press  SYSRESTART 

(25)  press  A/D  CLR 

(26)  if  no  changes  at  all,  go  to  Step  15.  If  there  are 
changes ,  do 

(27)  switch  to  "LOCAL" 

(28)  make  any  necessary  changes,  go  to  Step  14. 
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(D) 


If  system  hangs. 

do 

(29) 

switch  to  " 

LOCAL" 

(30) 

enter 

Address 

Content 

400 

FB 

401 

ED4D 

(31) 

switch  to  " 

REMOTE" 

(32) 

press  SYSRESTART 

(33) 

go  to  Part 

A  Step  2 

(34) 

if  Steps  29 

to  32  fail, 

drop  power  and  go  to 

Part  A  Step 

1 

3 . 1  Introduct ion 

Since  no  logic  analyzer  is  available  Co  debug  and  troubleshoot  the 
system,  it  is  necessary  to  have  a  maintenance  panel  in  the  system.  One 
third  of  the  entire  system  real-estate  was  devoted  to  the  maintenance  panel. 
There  are  many  built-in  maintenance  procedures  which  can  be  tested  in  local 
or  remote  mode.  It  is  hoped  that  this  panel  can  be  used  to  help  users  to 
narrow  down  any  specific  problems,  whether  they  are  hardware  or  software 
related. 

3.2  Maintenance  Panel  Design 

The  main  concerns  in  designing  the  panel  are: 

(1)  What  is  content  of  memory? 

(2)  What  is  content  of  address  line? 

(3)  Can  we  read/write  to  any  writable  memory? 

(4)  Can  we  program  all  the  PIO's  and  CTC? 

(5)  Can  we  test  all  the  programs  locally? 

(6)  Can  we  bind  individual  programs  to  form  a  master  program? 

(7)  Can  we  talk  to  the  outside  world? 

To  satisfy  all  these  concerns,  three  7-segment  LED's  are  used  to  display 
the  address  and  two  7-segment  LEDs  are  used  to  indicate  the  content  of  memory. 

A  switch  is  designated  to  control  remote  or  local  mode.  There  is  single-pulse 
switch  to  single  step  the  Z80.  The  memory  Write/memory  Read  switch  provides 
the  read  or  write  operation  to  the  static  RAM.  Memory  Address/Data  controls 
what  is  being  entered  through  the  Hex  keyboard.  There  are  reset  switches 
to  clear  the  internal  register  or  restart  the  system.  The  combination  of  all 
these  switches  provide  all  the  necessary  capabilities  to  test  the  system  locally 
Once  the  fundamental  system  is  operating,  simulated  programs  can  then 
be  Initiated  and  other  peripherals  can  be  tested.  This  is  the  concept  for 
testing  all  the  PIO's  and  CTC. 
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The  USERINTERRUPT1  Switch  Is  used  to  test  PIOs  and  CTC.  A  test  program 
for  a  particular  PIO  or  CTC  must  be  generated  and  written  into  the  memory 
starts  at  memory  location  A03.  The  length  of  the  program  cannot  exceed  128 
words.  The  procedures  for  running  this  test  will  be  described  in  the  next 
section. 


3.3  Maintenance  and  Testing  Procedures 

3.3.1  EPROM  Content  Check 

(1)  Power  up  system. 

(2)  LOCAL/REMOTE  at  LOCAL 

(3)  MEM/DATA  at  MEM 

(4)  MW  ON/OFF  at  OFF 

(5)  IOW  ON/OFF  at  OFF 

(6)  SYS  CLOCK/SINGLE  PULSE  at  SYSCLOCK 

(7)  Press  MEM  ADDR/DATA  REG  CLR 

The  first  three  digits  will  be  the  current  EPROM  address  and  the  next 
two  digits  is  the  content  of  this  address.  Press  UPCOUNT  or  DOWNCOUNT  to 
check  the  content  of  the  EPROM.  1 

3.3.2  RAM  Read  and  Write  Test 
(A)  Write  Test 

(1)  Power  up  system 

(2)  LOCAL/REMOTE  at  LOCAL 

(3)  IOW  ON/OFF  at  OFF 

(4)  SYS  CLOCK/SINGLE  PULSE  at  SYSCLOCK 

(5)  MW  ON/OFF  at  OFF 

(6)  MEM/DATA  at  MEM 

(7)  Press  MEMADDR/DATA  REG  CLR 

(8)  Enter  the  desired  address  to  be  written,  address  can  be  anyone 
400  to  4FF. 

(9)  MEM/DATA  at  DATA 

(10)  MW  ON/OFF  at  ON 

(11)  Press  MEMADDR/DATA  REG  CLR 

(12)  Enter  desired  data.  Each  datum  is  8  bits  long  and  is  divided 
into  two  4-bit  hexidecimal  code.  The  higher  order  4-bit  code  is 
enter  first  then  the  lower  order  4-bit  code.  The  memory  address 
will  automatically  increment  when  the  second  4-bit  is  entered. 


After  all  the  data  at  the  desired  addresses  are  write  we  can  go  back 
to  read  what  we  have  just  written. 

(13)  MW  ON/OFF  at  OFF 

(14)  MEM/DATA  at  MEM 

(15)  Press  MEM  ADDR/DATA  CLR 

(16)  Enter  the  desired  address  and  check  its  corresponding  content. 

3.3.3  Use  of  on-line  program  to  test  PIO's. 

If  3.3.1  to  3.3.2  are  working  correctly,  PIOs  and  CTC  can  then  be 
checked  on-line.  The  USERINTERRUPT  switch  is  intended  for  this  operation. 

The  procedure  is  same  as  memory  write.  The  format  is  as  follows: 

RAM  ADDR  CODE  PROGRAM 

403  CDD202  CALL  PUSHPAR 

406  Codes  for 

the  cor¬ 
responding 
program 

End  of  program 
C3AF02  JP  POPPAR 

3. 3. 3.1  Testing  the  Display  LED 

This  program  can  select  and  display  each  LED  individually. 


RAM  ADDR 

CODE 

PROGRAM 

403 

CDD202 

CALL  PUSHPAR 

406 

3E  XY 

LD  A,  XY 

408 

D3  11 

OUT  (11),  A 

40A 

06  02 

LD  B,  02 

40C 

11  FFFF 

LD,  DE,  -1 

40F 

2139DE 

Loopl: 

LD  HL,  14814 

412 

19 

Loop2 : 

ADD  HL,  DE 

413 

38FD 

JR  C,  L00P2 

415 

10F8 

DJNZ,  L00P1 

417 


C3AF0  2 


JP  POPPAR 
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x  ■  digit  select,  0  <_  x  <_  9 
y  ■  content  to  be  displayed  0  <_  y  <_  9 

Select  REMOTE 

Once  the  program  is  entered,  depressing  the  USER  INTERRUPT 
switch  once  should  display  the  selected  digit  with  the  desired  content 
for  approximately  2  seconds. 

3. 3. 3. 2  Test  MARK  statement  or  mark  the  tape  locally 


(1) 

Enter 

RAMADDR 

CONTENT 

4E6 

#  of  bytes  to 

be  marked 

4E7 

4E8 

#  records  to  be  marked 

4E9 

0  <  #  <  255 

RAMADDR 

CODE 

PROGRAM 

403 

CDD202 

CALL  PUSHPAR 

406 

CD7C03 

CALL  MARK  N,X 

409 

C3AF02 

JP  POPPAR 

(2) 

Select  REMOTE 

(3) 

Press  USERINTERRUPT 

3. 3. 3. 3  Test  FIND  statement  or  finding  the  tape  locally 
(1)  Enter 

RAM  ADDR  CONTENT 

4F7  Record  No.  for  user's  programmable 


4F8 


"FIND”  switch 


(2)  Select  REMOTE 

(3)  Press  FIND  RECORD  switch 


3. 3. 3. 4.  Test  FIND  statement  using  USER  INTERRUPT  switch 

(1)  Enter 


RAM  ADDR 

4F5 

4F6 


CONTENT 

Record  t  to  be  found 


RAM  ADD  CODE  PROGRAM 

403  CDD202  CALL  PUSHPAR 

406  CD5E03  CALL  FINDN 

409  C3AF02  CALL  POPPAR 

(2)  Select  REMOTE 

(3)  Press  USERINTERRUPT  switch 

3. 3. 3. 5  Test  WBYTE  statement  or  writing  data  using  USERINTERRUPT  switch 
(1)  Enter 


RAM  ADDR 

CONTENT 

4A0 

length  of 

data 

in  Binary 

4A1  to  X 

data,  x  £ 

what 

4A0  specified 

4F4 

79  ; 

this 

is  the  op  code 

RAM  ADDR 

CODE 

PROGRAM 

403 

CDD202 

CALL  PUSHPAR 

406 

21A004 

LD  HL,  A004H 

409 

3AA004 

LD  A,  (04AOH) 

40C 

47 

LD  B,  A 

40D 

3E79 

LD  A,  79 

40F 

32F404 

LD(04F4H) ,  A 

412 

CD6001 

CALL  HANSHAKE  1 

415 

CD3F03 

CALL  WBYTE 

418 

CD9601 

CALL  UNTKLN 

41B 

C3AF02 

CALL  POPPAR 

(2) 

Select 

REMOTE 

(3) 

Press 

USERINTERRUPT  1. 

3. 3. 3. 6  Conclusion 

Using  the  programs  described  above,  the  USERINTERRUPT1  switch 
can  perform  any  combination#  of  statements  like  FIND,  MARK  and  WYBTE. 
With  a  little  practice  and  imagination,  the  DBC  can  be  programmed  to 
perform  numerous  functions  like  generating  test  patterns  to  providing 
various  timing  signals.  I  shall  leave  this  for  my  users. 


Byron  C .  Cheung 
August  6,  1979  at  MIT 
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APPENDIX  E 

CROSS-CELL  VELOCITY  MEASUREMENTS 

In  order  to  quantify  the  three-dimensional  effects  occuring  at  the 
top  of  the  test -cell,  detailed  measurements  were  made  between  the  Plexi¬ 
glas  face  plates  on  a  horizontal  plane  located  approximately  at  the 
1*24  level  in  Fig.  C.l.  These  measurements  were  made  by  R.  W.  Sawdye  and 
are  recorded  in  Tables  E.l  and  E.2. 

The  data  show  a  definite  three-dimensional  structure  that  includes 
substantial  down-flow  occurring  near  the  face  plates  across  the  entire 
cell,  which  explains  the  failed  centerplane  mass  balance  that  is  described 
in  Chapter  Three.  It  appears  that  the  inlet  jet  impinges  on  the  cross¬ 
piece  that  forms  the  top  boundary  of  the  test  cell  and  rolls  outward  and 
down,  in  a  fairly  even  fashion,  along  both  face-plates.  Further,  the 
maximum  vertical  velocities  no  longer  occur  at  the  cell  centerplane,  but 
are  displaced  a  few  millimeters  to  one  side.  The  horizontal  velocities, 
on  the  other  hand,  generally  maintain  their  peaks  on  the  centerplane 
between  the  two  face  plates. 


I 
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Table  E.l 


Vertical  Velocities 


Measurements  by:  R.  W.  Sawdye 
FFTF  Test  Cell 

Flow  Rate  =  37.8  GPM  (steady-state) 

Temp.  =  110°F 

Vertical  Position,  I  *  24 


Off-Center  Distance  (mm) 


15 

10 

5 

0  5 

10 

15 

Approx. 

on 

Chimney 

Center 

/ 

/ 

-.361 

-.328 

.439 

.494  .577 

.441 

-.206 

/ 

-.336 

-.351 

.431 

.494  .604 

.408 

-.243 

/ 

/ 

-.347 

-.351 

.452 

.484  .573 

.368 

-.279 

/ 

/ 

-.323 

.338 

.455 

.474  .520 

.348 

-.311 

') 

-.332 

.305 

.468 

.461  .483 

.323 

-.322 

Plexiglas 

,  -.340 

.273 

.466 

.423  .437 

.276 

-.336 

Face 

-.329 

.111 

.429 

.366  .386 

.215 

-.344 

Plate 

-.322 

.169 

.334 

.303  .314 

.146 

-.343 

-.315 

.092 

.241 

.234  .232 

.078 

-.334 

-.287 

.028 

.157 

.158  .149 

.017 

-.308 

-.272 

-.029 

.082 

.092  .072 

-.037 

-.287 

-.262 

-.076 

.0i4 

-.025  -.007 

-.089 

-.279 

A 

-.293 

-.167 

-.064 

-.054  -.104 

-.171 

-.300 

-.332 

-.289 

-.219 

-.202  -.229 

-.273 

-.336 

-.316 

-.318 

-.322 

-.342  -.347 

-.336 

-.338 

Ann  my 

Cell  Side-Piece  Location 

/////'  /' 

/  /  y  /  f 

s  /  // s 

Table  E.l  Measured  Vertical  Velocities  (m/sec.)  on  a 

Horizontal  Plane  Between  Plexiglas  Face  Plates. 
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Table  E.2 

Horizontal  Velocities 


Measurements  by:  R.  W.  Sawdye 
Conditions:  Same  as  Table  E.l 


Off-Center  Distance  (mm) 


Approx . 
Chimney 

15 

10 

5 

0 

1 

1 

5 

10 

15 

Center  ~ 

y\ 

.019 

.021 

.034 

.022 

.013 

.032 

.064 

y 

.072 

.090 

.112 

.112 

.102 

.124 

.124 

'< 

X 

X 

.148 

.148 

.196 

.201 

.187 

.188 

.171 

y 

y 

X 

X 

.228 

.224 

.279 

.287 

.280 

.264 

.199 

y 

X 

X 

.275 

.287 

.365 

.374 

.361 

.306 

.239 

y 

Plexiglas 

.306 

.346 

.447 

.447 

.430 

.329 

.259 

y 

Face 

.343 

.358 

.508 

.501 

.479 

.352 

.259 

PI  ate 

.337 

.363 

.522 

.522 

.487 

.323 

.257 

•<| 

.313 

.306 

.468 

.496 

.456 

.270 

.239 

y 

.283 

.259 

.407 

.446 

.389 

.226 

.252 

y 

y 

.263 

.180 

.331 

.383 

.310 

.192 

.270 

y 

X 

X 

.263 

.150 

.241 

.321 

.245 

.178 

.292 

y 

y 

.275 

.162 

.184 

.243 

.203 

.194 

.313 

y' 

.253 

.191 

.180 

.190 

.192 

.211 

.300 

— 

.148 

.184 

.160 

.152 

.161 

.180 

.218 

Cell 

Side-Piece 

Location 

//////// 

rrrm 

-rrrn 

Table  E.l  Measured  Horizontal  Velocities  (m/sec.)  on  a 

Horizontal  Plane  Between  Plexiglas  Face  Plates. 
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APPENDIX  F 

PLOTS  OF  TURBULENCE  FIELD  PARAMETERS  VS.  TIME 

The  following  39  pages  of  plots,  one  page  for  each  I-J 

position  marked  with  a  (+)  in  Figs.  5.14  -  17,  portray  the 

turbulence  field  data  (+’s)  and  computed  transient  behavior  (solid 

lines)  as  a  function  of  time  during  the  transient.  Each  parameter, 

« 

both  data  and  calculation,  are  normalized  to  one  (1.0)  at  their 
individual  maximum  values  (usually  occurring  at  the  time  origin). 

These  normalization  constants  are  printed  above  each  graph  with  the 
convention  that  the  maximum  data  value  lies  over  the  maximum  calculated 
value.  The  labeling  scheme  is  demonstrated  on  the  first  page 
(I  8  5,  J  =  6). 

Note  that  turbulent  kinetic  energy  (K)  is  designated  (0)  on 
the  plots  in  Appendix  F.  Thus  QMAX  is  the  normalization  factor  for 
turbulent  kinetic  energy  as  a  function  of  time. 

Typical  values  of  the  standard  deviations  obtained  from  the 
ensembles  are  depicted  on  the  plots  of  page  252.  Precise  values  of 
the  standard  deviations  at  each  point  in  time  for  each  position  are 
given  in  Section  2  of  Appendix  C. 
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APPENDIX  G 


NOTES  ON  LDA  OPERATIONS 

A  number  of  general  statements  in  Chapter  Three  regarding  LDA 
theory  are  addressed  here  with  emphasis  on  the  specific  technique 
and  equipment  used  in  this  experiment.  In  particular,  this  appendix 
expands  on  the  following  three  sections: 

3. 2. 3. 2  Light  Scattering  Particles 

3. 2. 3. 3  Optical  Arrangement 

3. 2. 3. 4  Point  Positioning 


Light  Scattering  Particles 

In  cases  where  demineralized  water  with  artificial  seeding  is 
used,  the  scattering  particle  size  can  be  controlled  and  optimized. 
When  using  the  fringe  mode,  the  particle  diameter  should  be  smaller 
than  the  fringe  spacing  to  obtain  maximum  signal  modulation.  Gener¬ 
ally,  particle  diameters  on  the  order  of  ten  to  thirty  percent  of  the 
fringe  spacing  are  recommended.  In  the  reference  beam  mode,  the 
virtual  fringe  spacing  may  be  used  for  the  particle  diameter  calcula¬ 
tion.  The  half-micron  latex  particles  that  were  used  successfully  in 
a  previous  experiment  [Chen]  were  used  again  in  this  experiment,  and 
no  further  particle-size  optimization  v/as  attempted. 

The  scattering  particle  concentration  in  the  reference  beam 
mode  can  be  much  higher  than  that  used  in  the  fringe  mode.  In  the 


latter  case,  one  particle  per  measuring  volume  is  the  ideal  homogen¬ 
eous  concentration.  In  this  reference  beam  experiment,  a  concentra¬ 
tion  of  approximately  ten  to  one  hundred  particles  per  measuring 
volume  was  found  to  give  an  easily  tracked  and  continuous  Doppler 
signal.  The  so-called  "Doppler  ambiguity,"  introduced  by  the  simul¬ 
taneous  presence  of  many  particles  in  the  measuring  volume,  is 
discussed  and  accounted  for  in  Appendix  B. 

Optical  Arrangement 


In  order  to  measure  two  orthogonal  velocity  components  simultan¬ 
eously,  it  is  necessary  to  discriminate  between  the  Doppler  signals 
that  correspond  to  each  direction.  This  is  usually  accomplished 
by  using  either  different  polarization  or  different  wavelengths 
(colors)  for  the  two  velocity  directions  to  be  measured.  In  the 
dual-beam  or  fringe  mode  of  LDA  operation,  polarization  rotators  must 
be  added  to  the  optics  package  or  else  a  two-color  laser  modification 
(as  well  as  new  optics)  are  required.  The  expense  of  such  charges  is 
avoided  in  this  experimental  set-up  by  using  the  reference  beam  mode. 

In  such  a  case,  the  existing  optics  package  and  inexpensive  polarized 
lens  covers  on  each  of  the  two  photomultiplier  tubes  may  be  used  to 
provide  a  dual  channel  measurement  capability. 

It  is  noted  in  Section  3.2. 3.3  that  vibration  and  alignment 
problems  are  reduced  when  the  reference  beams  are  directed  to  pass 
through  the  measuring  volume.  Although  it  is  not  necessary  for  the 
reference  beams  to  pass  through  the  measuring  volume,  the  existing 
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optics  package  is  designed  to  operate  in  such  a  fashion.  To  split  off 
the  reference  beams,  guide  them  around  the  test-cell,  and  redirect 
them  accurately  down  the  bore  of  the  photo-multiplier  tubes  would  be 
an  unnecessary  complication.  Such  a  procedure  would  require  numerous 
mirrors  and  lenses  to  be  anchored  on  different  parts  of  the  experimental 
apparatus.  These  components  would  require  frequent  tedious  adjustments 
and  would  be  subject  to  different  sources  of  vibration.  The  existing 
set-up  anchors  the  entire  optics  package  at  a  single  point  which 
eliminates  the  problem  of  differential  vibrations  and  greatly  simpli¬ 
fies  alignment. 

Point  Positioning 

The  converted  milling-machine  used  for  point  positioning  and 
described  in  Sec.  3. 2.3.4  is  highly  accurate  when  used  properly. 
However,  a  possible  measurement  error  can  be  introduced  when 
reversing  the  direction  of  travel  because  of  play  in  the  threaded 
translation  mechanism.  This  phenomenon,  known  as  "backlash"  to 
experienced  machinists,  must  be  compensated  for  in  order  to  maintain 
positioning  accuracy. 

For  a  particular  translation  axis,  the  backlash  is  removed  as 
follows: 

(1)  Before  reversing  direction,  one  must  note  the  dial 
indicator  reading  and  then  turn  the  crank  slowly 
in  the  new  direction  until  the  play  is  taken  up  and 
the  threads  just  begin  to  bite  (generally  on  the 


order  of  10  hash-marks  on  the  indicator). 

(2)  Next,  one  loosens  the  knurled  knob  that  locks  the 
indicator  dial  and  rotates  the  dial  to  a  new  baseline 
on  the  indicator  ("0"  is  usually  most  convenient). 

(3)  Finally,  one  tightens  the  knurled  knob  and  continues 
cranking  in  the  new  direction  until  he  has  translated 
the  desired  distance. 

With  practice,  this  technique  can  provide  consistent  positioning 
accuracy  on  the  order  of  0.5  millimeters. 


